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The  new  directed  energy  and  advanced  kinetic  energy  weapons  have 
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1.  INTRODUCTION 


Recent  developments  In  directed  energy  and  In  advanced  kinetic 
energy  weapons  have  placed  a  renewed  emphasis  on  the  importance  of 
dynamic  fracture  processes  in  evaluating  the  vulnerability  of  airborne 
and  spaceborne  structures.  These  new  weapon  systems  can  subject  targets 
to  high  amplitude  load  pulses  of  very  high  rate  and  very  short 
duration.  Under  such  conditions,  the  concepts  of  classical  9tatic 
fracture  mechanics  can  no  longer  characterize  the  conditions  for 
initiation  and  propagation  of  preexisting  flaws  and  new  methods  must  be 
applied.  As  shown  in  a  previous  AFOSR-sponsored  program,^  the  highly 
transient  nature  of  the  stress  and  strain  fields  surrounding  the  crack 
tip  must  be  taken  into  account,  and  the  effect  of  extremely  high  loading 
rates  on  the  material  deformation  and  microdamage  response  must  be 
evaluated. 

This  report  summarizes  the  results  of  a  five-year  research  program 
to  improve  our  understanding  of  dynamic  fracture  behavior  by 
investigating  both  the  crack  tip  loading  conditions  and  the  material 
fracture  properties  for  high  loading  rates.  Section  2  reviews  the 
objectives  and  highlights  accomplishments  of  the  program.  The  results 
are  summarized  in  Section  3  and  discussed  in  detail  in  the  appendices. 
Section  4  lists  the  project  personnel  and  describes  special  activities 
during  the  program.  A  list  of  publications  and  presentations  resulting 
from  this  and  the  previous  programs  is  also  presented  in  Section  4. 


2.  RESEARCH  OBJECTIVES  AND  ACCOMPLISHMENTS 


The  objectives  of  this  research  were  to  establish  the  conditions 
controlling  dynamic  crack  initiation  and  extension  and  to  generate 
reliable  experimental  data  to  validate  proposed  dynamic  fracture 
theories*  To  accomplish  these  objectives,  the  following  tasks  were 
performed  during  the  five-year  program: 

Task  1:  Develop  a  simple  test  procedure  for  obtaining  reliable 

initiation  fracture  toughness  values  at  strain  rates 
representative  of  impact  loading. 

Task  2:  Using  the  test  procedure  developed  in  Task  1,  generate  the 

necessary  data  to  verify  and  extend  time  modified  theories 
for  dynamic  crack  instability. 

Task  3:  Measure  critical  conditions  and  establish  criteria  for 

crack  instability  under  dynamic  mixed-mode  I  and  II 
loading . 

Task  4:  Obtain  values  of  the  propagation  toughness  K™  and 

establish  the  relationship  between  and  tne  dynamic 
initiation  toughness  KI(j, 

The  successful  completion  of  these  tasks  has  resulted  in  the 
following  significant  accomplishments: 

•  A  simplified  impact  fracture  test,  the  one-point-bend  test 

(1PBT),  was  developed  and  fully  characterized.  The  new 
test  permits  reliable  dynamic  fracture  toughness 
measurements  at  rates  in  excess  of  10^  MPa  s-^.  This 

testing  procedure  is  now  documented  in  the  ASM  Metals 
Handbook,  Ninth  Edition,  Volume  8  on  Mechanical  Testing. 

•  The  1PBT  procedure  was  modified  to  permit  dynamic  mixed¬ 
mode  I  and  II  fracture  experiments  with  controlled  ratios 
of  mode  I  and  mode  II  loading. 

•  We  established  experimentally  that,  for  high  strength 
steels  and  similar  alloys,  no  time-modified  fracture 
criterion  is  necessary  when  the  loading  time  to  crack 
initiation  is  greater  than  20  p.s.  In  these  situations,  the 
classical  fracture  criterion  is  adequate  to  predict  crack 


initiation,  provided  the  correct  dynamically  applied  stress 
intensity  factor  is  equated  to  the  dynamic  fracture 
toughness. 

The  dynamic  mixed-mode  I  and  II  failure  curve  was 
determined  for  a  high  strength  steel  for  mode  II/raode  I 
ratios  beween  0  and  0.5.  Mixed-mode  fracture  in  this  range 
appears  essentially  controlled  by  the  mode  I  dynamic 
fracture  toughness.  A  lower  bound  mixed-mode  failure 
envelop  was  also  established  for  higher  mixed-mode  ratios. 

For  specimen  thicknesses  up  to  12.7  mm,  we  demonstrated 
that  the  propagation  toughness  of  high  strength  steel  can 
be  over  twice  the  value  of  the  dynamic  initiation 
toughness.  This  large  increase  in  toughness  was  attributed 
in  part  to  a  resistance  curve  effect  associated  with  the 
formation  of  small  shear  lips  and  with  the  bowing  of  the 
crack  front  during  extension. 


!>.V> 


m  ■'m 

.Vj 

/V.V.Y: 

r'M-'K- 


\  «■.  •-  <r . 

•  V  • 
\V.Y 


~7 

A  v. 


VliMVwV AIVWIY  V."%i  lTw  UV  iT'm  L’V  V7V  mT<  VXWl 


3.  SUMMARY  OF  RESULTS 


The  results  of  the  research  are  summarized  below  according  to  the 
four  main  tasks  of  our  work.  We  first  report  on  the  development  of  new 
experimental  techniques  necessary  to  investigate  fracture  at  high  loading 
rates.  Then  we  review  how  these  techniques  were  used  to  investigate 
dynamic  crack  initiation  under  pure  mode  1  loading,  mixed-mode  I  and  II 
loading,  and  dynamic  crack  propagation.  Detailed  discussions  of  each 
topic  are  presented  in  Appendices  A  through  H. 
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DEVELOPMENT  OF  NEW  DYNAMIC  FRACTURE  TESTING  TECHNIQUES. 

One  main  achievement  of  the  research  program  was  the  development  and 
analysis  of  a  novel  impact  fracture  test,  the  one-point-bend  test  (1PBT), 
for  measuring  reliable  dynamic  initiation  fracture  toughness.  The 
concept  of  the  LPBT  was  used  further  to  develop  an  impact  test  capable  of 
producing  crack  Initiation  under  controlled  conditions  cf  mixed  mode  I 
(tensile  opening)  and  mode  II  (in-plane  shear)  loading  and  to  investigate 
dynamic  crack  propagation. 

The  1PBT  for  Mode  I  Dynamic  Crack  Initiation 

A  testing  method  that  produces  an  adjustable  and  smoothly  varying 
load  history  is  desirable  to  obtain  reliable  measurements  of  the  dynamic 
initiation  fracture  toughness.  Unfortunately,  conventional  tests  such  as 
the  Charpy  impact  test  or  the  dynamic  tear  test  do  not  produce  a  smoothly 
varying  stress  Intensity  history.  Thus,  interpretation  of  the  test 
results  is  not  straightforward.  To  circumvent  these  difficulties  we 
investigated  a  test  configuration  in  which  a  simple  edge-cracked  specimen 
Is  loaded  in  bending  by  impacting  It  at  the  midsection,  but  without 
supporting  it  at  the  outer  edges. 
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The  test  arrangement  and  the  loading  mechanism  for  this  one-point- 
bend  test  are  shown  in  Figure  1.  The  loading  at  the  crack  tip  is 
achieved  strictly  by  inertial  effects.  The  resulting  stress  intensity 
history,  shown  in  Figure  2,  has  approximately  a  sinusoidal  time 
dependence,  with  the  period  of  oscillation  given  by  the  first  natural 
frequency  of  vibration  of  the  specimen.  The  maximum  stress  intensity 
amplitude  K^ax  is  controlled  in  part  by  the  specimen  geometry  and 
density,  but  mainly  by  the  impact  velocity.  The  9tress  intensity  history 
can  easily  be  measured  during  the  test  by  placing  a  strain  gage  a  short 
distance  from  the  crack  tip. 

We  investigated,  analytically  and  experimentally,  how  the  stress 

Max 

intensity  duration  TQ  and  amplitude  vary  with  impact  conditions, 

specimen  geometry,  dimensions  and  material  properties  (Appendix  A  and 
Appendix  B) .  We  showed  that  TQ  and  K^ax  can  be  expressed  as : 

TQ  «=  Wv/p7e  h  (L/W,  a/W,  v)  (1) 

^  *  Vimp  VL/W*  a^w»  v>  (2) 


where  W,  L,  and  a  are  the  width,  length,  and  crack  depth  of  the  specimen, 
respectively;  p,  E,  and  v  are  the  mass  density,  Young's  modulus,  and 
Poisson's  ratio,  respectively;  Vimp  is  the  impact  velocity;  and  h^  and  h£ 
are  nondlmenslonal  functions.  We  also  investigated  the  reliability  of 
the  strain  gage  technique  for  measuring  the  stress  intensity  history  and 
established  procedures  for  determining  the  point  of  crack  initiation 
(Appendix  C) • 

The  1PBT  can  be  performed  on  a  pendulum  impact  tester,  on  a  drop- 
weight  tester,  or  on  any  other  apparatus  capable  of  impact  velocities 
greater  than  3  ra  s-^.  Its  primary  advantage  is  that  the  measured  stress 
intensity  history  varies  smoothly  and  incorporates  dynamic  effects 
completely.  Therefore,  in  contrast  to  other  dynamic  fracture  tests,  no 
limits  need  to  be  imposed  on  the  impact  velocity  and  the  test  duration  to 
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fracture.  Another  advantage  of  the  1PBT  is  that  the  test  duration  TQ  and 

* 

the  stress  intensification  rate  Kj  can  be  varied  over  an  order  of 
magnitude  range  (80  to  800  ps  and  3  x  10^  to  3  x  10^  MPa  m^^  s-^)  by 
changing  the  specimen  geometry  and  the  impact  velocity. 

The  1PBT  is  now  a  well-established  technique  for  simply  and  reliably 
measuring  dynamic  initiation  fracture  toughness  values.  Its  use  has  also 
been  extended  recently  to  nonmetallic  materials,  and  it  has  been 
successful  in  discriminating  between  two  commercial  alumina  ceramics  in 
terms  of  their  dynamic  fracture  toughness. 

The  Modified  1PBT  for  Dynamic  Mixed-Mode  I  and  II  Crack  Initiation 

To  investigate  dynamic  crack  initiation  under  mixed-mode  loading 
conditions,  it  is  necessary  to  control  the  mixed-mode  ratio  Kjj/Kj  that 
is  achieved  during  the  experiment.  Controlling  this  ratio  can  be 
particularly  difficult  because,  due  to  structural  and  stress  wave 
dynamics,  the  mixed-mode  ratio  generally  varies  continuously  in  a  dynamic 
mixed-mode  experiment.  To  overcome  this  difficulty,  we  chose  the  double¬ 
crack  1PBT  configuration  shown  in  Figure  3(a)  (Appendix  F).  This 
configuration  presents  the  advantage  that  the  specimen  dynamics  are 
relatively  simple  and  easy  to  understand:  inertia  will  causa  shear 
stresses  on  the  ligaments  ahead  of  the  off-centered  cracks  and  bending  as 
in  the  single  crack  one-point-bend  test.  By  controlling  the  shear 
stresses,  we  can  adjust  the  mode  II  stress  intensity,  and  by  controlling 
the  bending,  we  can  adjust  the  mode  I  stress  intensity.  The  straight 
crack  configuration  also  makes  it  easier  to  fatigue  precrack  the 
specimens  than  it  is  with  the  slanted  crack  configurations  used  in  other 
Investigations . 

To  achieve  contol  of  the  mixed-mode  ratio  kxj/kI*  We  modified  the 
configuration  in  Figure  3(a)  by  statically  prestressing  the  specimen 
using  the  arrangement  shown  in  Figure  3(b) .  Prestressing  introduces  a 
mode  I  loading  component  during  the  early  times  of  the  experiments  when 
the  mode  II  loading  due  to  the  inertia  of  the  specimen  is  maximum.  The 
inertia  of  the  prestressing  fixture  itself  permits  maintaining  this  mode 


II  loading  component  at  later  times  when  the  dynamic  mode  I  loading 
component  due  to  specimen  inertial  bending  becomes  significant.  Stress 
intensity  histories  obtained  with  this  arrangement  are  shown  in  Figure  4 
for  an  experiment  in  which  the  specimen  fractured  at  a  mixed-mode  ratio 
of  0.51. 

We  investigated  the  dependence  of  the  stress  intensity  histories  In 
the  double-crack  one-point-bend  specimen  on  the  crack  length,  the  crack 
separation,  and  the  impact  conditions.  Three  strain  gages  were 
positioned  at  different  angles  near  the  crack  tip  to  measure  the  stress 
intensity  histories  (Appendix  C) .  It  was  found  that  beyond  a  limit 
Impact  velocity  for  which  vibrational  model  transition  occurs,  the  stress 
intensity  amplitudes  are  proportional  to  the  impact  velocity  and  increase 
with  increasing  initial  crack  length.  The  distance  between  the  two 
cracks  has  little  effect  on  the  mode  II  stress  intensity  and  is  chosen  to 
facilitate  fatigue-precracking. 

With  the  arrangement  shown  in  Figure  3(b) ,  we  were  able  to  produce 
controlled  stress  intensity  histories  with  mixed-mode  ratios  varying  from 
0  to  0.5,  with  maximum  Kjj  amplitudes  in  excess  of  50  MPa  m1^2* 

Therefore,  a  simple  and  reliable  technique  is  now  available  for  dynamic 
mixed-mode  fracture  investigations.  Further  studies  are  required  to 
fully  understand  the  doub!  e-crack  1PBT  specimen  dynamics.  In  particular 
the  effect  of  the  impact  ■ ilocity  on  the  specimen  response  needs  to  be 
examined.  Moreover,  adjustment  of  the  prestres.  ~ng  support  mass  and 
position,  as  well  as  modifications  of  the  impact  har  er  geometry,  will 
result  in  increased  mixed-mode  ratios  and  improved  control  of  the 
experiment. 


The  one-point-bend  test  Is  also  attractive  for  crack  propagation 
investigations.  The  simplicity  of  inertial  crack  tip  loading  and  the 
good  control  over  the  loading  amplitude  and  duration  permit  the  crack 
extension  distance  and  the  crack  velocity  to  be  adjusted  precisely. 
Further,  the  simplicity  of  the  boundary  conditions  in  the  1PBT  makes 
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numerical  simulations  of  the  experiment  particularly  easy.  Finally,  the 
one-point-bend  test  permits  the  investigation  of  propagating  cracks 
initiated  under  impact  conditions  in  contrast  to  most  other  test 
procedures,  which  use  quasi-static  initiation  conditions. 

During  this  program,  we  demonstrated  that  very  small,  controlled 
crack  jumps  can  be  produced  with  the  1PBT  (Appendix  G) .  With  this 
technique  we  were  able  to  investigate  the  very  first  increments  of  crack 
extension  to  relate  dynamic  initiation  toughness  and  propagation 
toughness.  We  developed  an  empirical  model  of  the  1PBT  to  obtain  first 
order  estimates  of  the  stress  intensity  at  the  tip  of  the  extending  crack 
(Appendix  D).  The  model  is  based  on  the  stress  intensity  history 
measurements  for  specimens  with  stationary  cracks  of  various  lengths.  It 
assumes  that  the  1PBT  is  displacement  controlled,  with  the  first  flexural 
mode  of  the  specimen  controlling  the  stress  intensity  at  the  crack  tip. 
The  model  accounts  for  the  effect  of  crack  extension  and  for  the  effect 
of  material  inertia  on  the  stress  intensity  at  the  crack  tip.  The  model 
is  implemented  on  a  microcomputer  and  calculates  the  stress  intensity  and 
history  using  the  propagation  toughness  or  the  crack  velocity  as  input 
variables. 

The  model  was  used  to  evaluate  the  first  increments  of  crack 
extension  during  the  1PBT  experiments.  We  also  used  it  to  analyze  the 
strain  history  at  the  original  crack  tip  after  crack  initiation  (Appendix 
C,  Appendix  G  ),  and  we  demonstrated  that  the  strain  history  can  yield 
crack  velocity  and  propagation  toughness  information. 

DYNAMIC  MODE  I  CRACK  INITIATION  FOR  SHORT  LOADING  TIMES 

In  a  previous  AFOSR-sponsored  research  program,^  it  was 
experimentally  established  that,  for  very  short  duration  pulse  loads 
(less  than  2. 2-^8  duration),  the  classical  quasi-static  linear  elastic 
fracture  mechanics  criterion  for  crack  instability  must  be  modified  to 
account  for  microstructural  rate  effects.  Kalthoff  and  Shockey^ 
postulated  that  cracks  would  become  unstable  only  after  the  dynamic 
fracture  toughness  has  been  exceeded  for  a  minimum  time  (minimum  time 


theory).  Based  on  fracture  tests  with  pulse  loads  of  several  durations 
(18  to  80  ps),  Homma  et  al.^  have  estimated  that,  for  4340  steel  in  the 
HRC  50  condition,  the  minimum  time  for  fracture  should  be  about  7  ps. 

This  value  is  somewhat  uncertain,  however,  because  the  fracture  toughness 
reported  by  Homma  et  al.  is  only  about  half  the  value  reported  in  the 
literature  for  similar  loading  rates. 

More  recently,  Ravi-Chandar  and  Knauss^  in  experiments  with 
Homalite,  and  Kalthoff^  in  experiments  with  high  strength  steel,  found 
that  when  the  time  to  monotonically  load  a  crack,  to  fracture  is  gradually 
decreased  below  about  30  ps,  the  stress  intensity  value  at  which  the 
crack  becomes  unstable  increases  rapidly.  For  a  loading  time  to  fracture 
of  5  to  10  ps,  the  stress  intensity  factor  at  instability  is  about  twice 
the  value  when  the  loading  time  is  50  ns  or  more.  Kalthoff  rationalized 
this  behavior  by  suggesting  that  a  so-called  incubation  time  is  necessary 
to  nucleate,  grow,  and  coalesce  microdefects  ahead  of  the  macrocrack. 
Conceptually,  the  incubation  time  is  similar  to  the  minimum  time 
postulated  by  Kalthoff  and  Shockey. 

We  attempted  to  establish  whether  fracture  incubation  time  has  a 
significant  influence  on  crack  instability  conditions  in  high  strength 
steel  at  conventional  impact  loading  rates,  that  is,  for  fracture  times 
down  to  20  pb  (Appendix  E).  For  this  investigation,  we  used  the  1PBT, 
varying  the  impact  velocity.  In  tests  at  low  stress  intensity  rate,  the 
specimens  were  impacted  at  a  velocity  sufficient  to  produce  a  maximum 
stress  intensity  slightly  above  the  anticipated  critical  value.  The 
loading  times  to  fracture  were  about  100  ps,  well  above  the  anticipated 
incubation  time.  Thus  these  tests  provided  a  measure  of  the  dynamic 
fracture  toughness,  KI(j. 

Tests  at  high  loading  rates  were  designed  to  determine  a  possible 
incubation  time.  In  these  tests  the  specimens  were  grossly  overloaded  by 
increasing  the  impact  velocity  to  achieve  times  to  fracture  of  20  to  25 
ps.  Thus,  if  fracture  incubation  time  is  important,  and  of  the  order  of 
several  microseconds,  we  expected  fracture  to  occur  at  significantly 
higher  stress  intensity  values  than  the  value  measured  in  the  low  rate 


tests.  The  difference  between  the  time  at  which  K-^  was  reached  and  the 
time  of  instability  would  provide  an  estimate  of  the  incubation  time. 
Intermediate  rate  tests  provided  additional  data  to  monitor  the  influence 
of  the  loading  rate  on  the  value  of  the  stress  intensity  factor  at  crack 
initiation. 

The  results  of  the  dynamic  fracture  teats  are  summarized  in  Figure 
5.  Within  the  range  of  loading  rates  studied,  the  stress  intensity 
factor  at  initiation  is  58.5  MPa  and  essentially  constant.  This 

value  is  about  92%  of  the  static  toughness  and  represents  the  dynamic 
initiation  fracture  toughness  KId*  The  results  therefore  demonstrate 
that  a  fracture  incubation  time  is  unimportant  for  these  loading  rates. 
Further,  our  results  imply  that,  if  a  fracture  incubation  time  exists,  it 
is  much  smaller  than  20  to  25  Us. 

We  conclude  that,  for  fracture  times  in  excess  of  20  ns,  crack 
initiation  predictions  for  4340  steel  in  the  HRC  50  condition  (and 
probably  for  similar  high  strength  alloys)  do  not  require  a  time  modified 
fracture  criterion.  Rather,  our  work  indicates  that  the  initiation 
conditions  can  be  determined  solely  by  using  the  classical  concepts  of 
fracture  mechanics,  provided  the  correct  dynamic  stress  intensity  history 
is  used. 

DYNAMIC  MIXED-MODE  CRACK  INITIATION 

In  most  practical  situations,  cracks  are  not  only  loaded  in  the 
tensile  opening  mode,  i.e.,  mode  I  loading,  but  also  experience  In-plane 
shear,  i.e.,  mode  II  loading.  Crack  initiation  conditions  under  mixed¬ 
mode  loading  no  longer  depend  on  a  single  parameter,  but  rather  must  be 
expressed  in  terms  of  a  failure  curve  in  the  Kj-KTj  plane,  where  Kj  and 
Kjj  represent  the  mode  I  and  mode  II  stress  intensity  factors, 
respectively.  Several  attempts  have  been  made  to  establish  the  mixed¬ 
mode  failure  locus  under  quasi-static  loading  conditions.  These 
investigations  suggest  that  a  failure  curve  based  on  a  critical  maximum 
tangential  stress  at  a  critical  distance  ahead  of  the  crack  tip  may  be 
adequate-  However,  they  also  suggest  that  the  failure  locus  depends 


on  the  specific  material  and  on  the  raicromode  of  fracture.  In  contrast 
to  the  static  mixed-mode  case,  very  few  data  are  available  for  dynamic 
mixed-mode  crack  initiation,  and  the  lack  of  these  data  prevents  the 
development  and  evaluation  of  reliable  mixed-mode  initiation  criteria. 

To  remedy  this  lack  of  information,  we  performed  mixed-mode  Initiation 
experiments  using  the  impact  test  technique  to  establish  the  failure 
curve  for  4340  steel  in  the  HRC  50  condition  (Appendix  F) . 

The  results  of  these  experiments  are  summarized  in  Figure  6,  which 
plots  the  mode  1  stress  intensity  factor  normalized  by  the  dynamic  mode  I 
initiation  toughness,  versus  the  mode  II  stress  intensity  factor, 

also  normalized  by  In  the  dynamic  experiments  the  mixed-mode  ratio 

does  not  stay  constant  but  varies  continuously  with  time.  The  open 
circles  in  Figure  6  therefore  represent  mixed-mode  states  that  have  been 
traversed  during  the  tests,  but  did  not  cause  initiation;  the  solid 
circles  represent  the  mixed-mode  states  that  resulted  in  crack 
instability.  The  lines  connecting  the  symbols  represent  typical  loading 
paths  followed  during  the  experiments. 

We  draw  several  conclusions  from  the  results  in  Figure  6.  First, 
for  Kxi^KI  between  0  atld  0.5,  crack  initiation  appears  essentially 
controlled  by  the.  mode  I  stress  intensity  factor.  Initiation  occurs  when 
KI  M  KId*  (The  two  isolated  points  on  the  far  right  of  Figure  6  are 
believed  to  be  shifted  artificially  to  higher  values  of  Kj  at  initiation 
because  of  the  difficulty  of  detecting  the  precise  point  of  crack 
initiation  for  higher  mixed-mode  ratios,  see  Appendix  F.) 

Further,  based  on  the  points  of  the  Kj-Kj-j-  plane  traversed  during 

the  fracture  experiments,  but  for  which  no  failure  occurred,  a  tentative 

lower  bound  mixed-mode  failure  curve  can  be  drawn  for  4340  steel  HRC 

50.  We  estimate  that  the  pure  mode  II  dynamic  toughness  Kj-jj  is  at  least 

equal  to  the  pure  mode  1  initiation  toughness  K^;  i.e.,  >_  58  MPa 

1  /2 

m  '  .  Assuming  that  Kj-j^  has  the  same  low  rate  sensitivity  as  Kj^,  and 
taking  for  the  static  mode  II  fracture  toughness  the  value  reported  by 
Shah,^  88  MPa  we  estimate  that  Kjjj  should  be  about  85  MPa  m^^. 

This  value  is  consistent  with  the  proposed  lower  bound  failure  curve. 
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Failure  Curve 


FIGURE  6  RESULTS  OF  DYNAMIC  MIXED-MODE  FRACTURE  EXPERIMENTS  ON  4340  STEEL 
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These  results  are  shown  in  Figure  7,  which  plots  as  a  function 
of  the  crack  velocity,  as  is  customary  in  crack  propagation 
investigations,  to  show  the  strong  crack  velocity  dependence  of  the 
propagation  toughness*  Also  shown  in  Figure  7  are  results  for  the  same 

0  Q 

material  published  in  the  literature.  *  At  the  crack  velocities 
investigated,  our  results  fall  significantly  above  the  published 
results.  However,  Figure  7  is  somewhat  misleading  because  part  of  the 
elevation  of  the  resistance  to  crack  extension  above  the  initiation 
toughness  is  due  to  a  resistance  curve  effect,  whereby  the  propagation 
toughness  increases  above  the  initiation  toughness  as  a  function  of  crack 
extension  distance.  The  mechanism  for  this  elevation  is  the  formation  of 
small  shear  lips  and  the  bowing  of  the  crack  front,  which  occur  during 
the  first  millimeter  or  so  of  crack  extension  in  the  relatively  thin 
specimens  used  for  our  investigation  (9.8  and  12.7  mm).  Moreover,  we 
believe  that  the  velocity  dependence  of  the  propagation  toughness  is  also 
influenced  by  the  formation  of  shear  lips  a  id  by  crack  tunneling,  but 
this  crack  velocity  dependence  could  not  be  assessed  independently  from 
the  resistance  curve  effect. 

Our  findings  have  four  practical  implications.  First,  they  indicate 
that  cracks  initiated  by  pulse  loads  in  structures  less  than  25  mm  thick 
may  not  necessarily  be  unstable,  because  of  the  resistance  curve  effect, 
which  requires  an  increasing  applied  stress  intensity  factor  to  sustain 
crack  extension  at  an  aopreciable  crack  velocity.  The  second  implication 
is  that  the  thickness  requirements  to  obtain  thickness  independent 
propagation  toughness  values  are  much  more  stringent  than  the  thickness 
requirements  for  the  measurement  of  dynamic  initiation  toughness.  The 
third  implication  is  that,  because  of  the  three-dimensional  nature  of  the 
crack  propagation  process,  caution  should  be  applied  in  interpreting 
experimental  data  obtained  from  surface  measurements.  Finally,  a 
comparison  of  propagation  toughness  data  obtained  with  cracks  initiated 
from  blunt  notches  or  from  sharp  cracks  suggests  that  the  initiation 
conditions  may  affect  the  extension  behavior  over  distances  that  are 
significant  compared  with  the  plastic  zone  size  at  the  initial  notch  or 
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REFERENCE:  Giovanola,  J.  H..  “Investigation  and  Application  of  the  One -Point  Bend 
Impact  Test,’"'  Fracture  Mechanics  Seventeenth  Volume,  ASTM  STP  90S,  1.  H  Under¬ 
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ABSTRACT:  This  paper  investigates  a  dynamic  fracture  test  that  produces  smoothly 
varying  stress  intensity  histories  of  controlled  amplitude  and  duration.  A  test  configura¬ 
tion  has  been  adopted  in  which  a  simple  edge-cracked  specimen  is  loaded  in  bending  by 
impacting  it  al  the  midsection  without  supporting  it  at  the  outer  edges.  The  stress  inten 
sity  history  is  measured  by  a  strain  gage  near  the  crack  tip  and  has  approximately  a  sinus¬ 
oidal  time  dependence 

The  dependence  of  stress  intensity  amplitude  and  duration  on  the  impact  velocity  and 
on  the  specimen  material  and  dimensions  has  been  determined  both  analytically  and  ex¬ 
perimentally.  The  limit  on  the  specimen  aspect  ratio  to  obtain  a  smooth  stress  intensity 
history  controlled  only  by  the  first  mode  of  vibration  of  the  specimen  has  also  been  estab¬ 
lished.  Stress  intensity  pulses  of  durations  varying  between  90  and  500  ps  can  be  readily 
achieved,  and  the  maximum  amplitude  can  be  independently  adjusted  over  an  order  of 
magnitude  (20  to  200  MPa  m1  }). 

The  new  experimental  procedure  has  been  used  to  measure  the  dynamic  fracture  tough¬ 
ness  of  4340  steel  (HRC  50)  at  three  loading  rates.  Comparison  of  the  dynamic  results  with 
the  value  of  the  static  toughness  indicates  little  strain -rate  sensitivity  of  the  fracture 
toughness  of  4340  steel  in  the  range  of  loading  rales  investigated  (quasi-static  to  3  x  10* 
MPa  m‘  Js  '). 

KEY  WORDS:  dynamic  fracture,  stress  intensity  history,  impact  bend  test,  unsupported 
bend  specimen,  dynamic  fracture  toughness.  4340  steel,  raie  sensitivity 


Nomenclature 

a  Crack  length 
a  Crack  velocity 
B  Specimen  thickness 

CcR  Elastic  compliance  of  bend  specimen  due  to  crack 
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Elastic  compliance  of  uncracked  bend  specimen 
Specimen  deflection 

Specimen  deflection  along  center  line  due  to  crack 
Uncracked  specimen  deflection  along  center  line 
Young's  modulus 

Maximum  kinetic  energy  of  specimen 
Maximum  strain  energy  of  specimen 
Dynamic  stress  intensity  factor 
Static  fracture  toughness 
Dynamic  initiation  fracture  toughness 
Maximum  stress  intensity  amplitude 
Specimen  mass 
Hammer  mass 
Hammer  load 
Maximum  hammer  load 

Maximum  effective  dynamic  load  on  the  specimen 

Specimen  length 

Time 

Time  at  which  maximum  stress  intensity  amplitude  is  reached 

Duration  of  stress  intensity  history 

First  natural  period  of  oscillation  of  specimen 

First  natural  period  of  oscillation  of  the  uncracked  specimen 

Hammer  velocity  at  fMAX 

Hammer  impact  velocity 

Velocity  of  the  specimen  center  of  mass  at  rMA* 

Specimen  wndth 

Coordinate  along  specimen  axis 
Mass  density  of  specimen  material 
Poisson’s  ratio 


Introduction 

A  dynamic  fracture  test  having  an  adjustable  and  smoothly  varying  load 
history  is  desirable  for  unambiguous  determination  of  dynamic  fracture  be¬ 
havior.  Unfortunately,  conventional  tests,  such  as  the  Charpy  impact  test  or 
the  dynamic  tear  test,  do  not  produce  a  smoothly  varying  stress  intensity  his¬ 
tory.  Thus  interpretation  of  the  results  is  not  straightforward.  This  paper  dis¬ 
cusses  a  modification  of  the  conventional  impact  bend  test,  first  proposed  and 
investigated  by  Kalthoff  and  co-workcrs  | /-.?],  which  produces  a  smoothly 
varying  stress  intensity  history.  The  amplitude  and  duration  of  the  stress  in¬ 
tensity  history  are  easily  varied  by  varying  the  impact  velocity  and  the  speci¬ 
men  dimensions. 

Wc  have  investigated  the  influence  of  specimen  dimensions  and  impact  ve¬ 
locity  on  the  stress  intensity  history  both  analytically  and  experimentally.  We 
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have  then  used  the  new  test  procedure  to  measure  the  dynamic  fracture 
toughness  of  4340  steel  at  three  strain  rates. 

Background 

Dynamic  effects  in  impact-loaded,  three-point-bend  specimens  (3PBS) 
have  been  discussed  in  several  recent  publications.  In  particular,  Kalthoff 
and  co-workers  have  shown  that  the  stress  intensity  history  results  from  the 
superposition  of  the  specimen  free  oscillations  on  the  fundamental  oscillation 
of  the  hammer-specimen-support  system  [1,2, 4-6 1.  Bbhme  (6)  has  proposed 
a  semi-analytical  description  of  the  behavior  of  the  hammer-specimen- 
support  system  to  predict  the  influence  of  specimen  dimensions,  impact  ve¬ 
locity,  and  material  properties  on  the  stress  intensity  history.  He  has  further 
demonstrated  that  the  dynamic  effects  due  to  the  specimen  oscillations  can  be 
estimated  for  specimens  of  arbitrary  size  and  material  from  a  set  of  normal¬ 
ized  master  curves.  Each  master  curve  depends  only  on  the  ratios  of  crack 
length  to  width,  a/W,  and  specimen  length  to  width,  S/W. 

Bohme  and  Kalthoff  [/]  have  also  observed  that  soon  after  impact  the  3PB 
specimen  momentarily  loses  contact  with  the  supports.  Therefore,  during  the 
initial  phase  of  the  test,  the  supports  do  not  influence  the  loading,  and  the 
stress  intensity  history  obtained  for  specimens  impacted  with  or  without  sup¬ 
ports  is  identical.  The  crack  tip  loading  is  then  strictly  achieved  by  inertia. 
Only  at  later  times,  when  the  specimen  makes  contact  with  the  supports 
again,  does  its  behavior  differ  from  that  of  the  unsupported  specimen;  the 
stress  intensity  history  for  the  supported  specimen  continues  to  rise,  whereas 
that  for  the  unsupported  specimen  decreases  to  zero. 

On  the  basis  of  these  observations,  Kalthoff  and  co-workers  [ 2,3}  have 
demonstrated  that  the  bend  specimen  impacted  without  supports— called 
one  point -bend  test  (1PBT)— provides  a  convenient  configuration  to  obtain 
controlled  and  smoothly  varying  stress  intensity  histories. 

The  test  arrangement  and  loading  mechanism  for  the  1PBT  are  illustrated 
in  Fig.  1.  Figure  2  shows  the  resulting  stress  intensity  history  in  a  specimen 
with  a  blunt  notch,  as  measured  with  a  strain  gage  placed  near  the  notch  tip. 
This  stress  intensity  history  is  roughly  sinusoidal  and  is  characterized  by  the 
half  period  7o  and  the  maximum  amplitude  AT^ax.  Figure  2  also  shows  the 
corresponding  hammer  load.  On  impact,  the  center  portion  of  the  specimen 
is  accelerated  away  from  the  hammer,  whereas  the  end  portions  of  the  speci¬ 
men  lag  behind  because  of  inertia.  This  causes  the  specimen  to  bend  and  to 
load  the  crack  tip.  The  specimen  deflection  also  progressively  unloads  the 
hammer  tip  until  the  maximum  deflection  (and  hence  the  maximum  stress 
intensity)  is  reached.  As  the  specimen  bends  back,  the  stress  intensity  de¬ 
creases,  whereas  the  hammer  load  increases  again  because  the  relative  dis¬ 
placement  of  the  hammer  and  of  the  specimen’s  center  portion  are  now  in 
opposite  directions  (/]. 
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(si  Before  Impact  (b)  After  Impact 

FIG.  1 — Test  arrangement  and  hading  mechanism  for  the  IPBT. 


TIME,  i* 

FIG.  2 — Typical  stress  intensity  and  load  histories  obtained  in  IPBT. 


According  to  this  interpretation  of  the  stress  intensity  records,  the  first  nat¬ 
ural  period  of  oscillation  of  the  specimen  controls  the  duration  of  the  approxi¬ 
mately  sinusoidal  stress  intensity  history.  Thus  the  test  duration  can  be  con¬ 
trolled  by  changing  the  dimensions  of  the  specimen  to  change  the  natural 
frequency.  Furthei  -nore,  the  results  of  Bohme  [6]  suggest  that  the  maximum 
stress  intensity  amplitude  is  proportional  to  the  impact  velocity. 

The  purpose  of  the  present  investigation  was  to  verify  the  foregoing  predic¬ 
tions,  to  establish  the  range  of  load  durations  that  could  be  achieved  with  the 
IPBT,  to  determine  more  precisely  the  role  of  specimen  dimensions  and  ma¬ 
terial  properties  on  the  stress  intensity  history,  and  to  apply  the  IPBT  to  dy¬ 
namic  fracture  toughness  testing. 
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Investigation  of  the  Dynamic  Behavior  of  the  1PBT 


Analysis 

Here  we  assume  that  the  stress  intensity  history  is  sinusoidal  with  a  period 
27o  equal  to  the  first  natural  period  of  the  specimen  and  with  an  amplitude 
^•max  tjlen  determine  how  these  two  parameters  depend  on  the  specimen 
material  and  dimensions  and  on  the  impact  velocity. 

Details  of  this  approximate  analysis  are  presented  in  the  Appendix  together 
with  references  to  earlier  analyses.  We  find  that  To  and  A'|’1ax  are  given  by 


To  =  Ty 


*■(£•  -W'  yh'vyliJ,2(ww' 


(i) 


and 


*rX  =  V,mp  yfW  s^E  hj^.  ~  v)  (2) 

where  Tur)  is  the  first  natural  period  of  oscillation  of  the  uncracked  beam  with 
free  ends;  A, ,  A2,  and  A,  are  dimensionless  functions  of  the  specimen  dimen¬ 
sions  and  of  Poisson's  ratio;  and  the  other  symbols  are  as  defined  in  the  No¬ 
menclature  section  of  this  paper.  In  Eq  2,  it  has  been  assumed  that  the  ham¬ 
mer  mass,  M ,  is  much  larger  than  the  specimen  mass.  m. 

Equations  1  and  2  show  that  for  geometrically  similar  specimens  (and  ne¬ 
glecting  a  mild  dependence  on  Poisson's  ratio),  To  scales  with  the  specimen 
size  and  with  the  square  root  of  the  ratio  of  mass  density  to  elastic  modulus; 
K is  proportional  to  the  impact  velocity  and  scales  with  the  square  root  of 
specimen  size,  mass  density,  and  elastic  modulus.  Both  To  and  are  in¬ 
dependent  of  specimen  thickness  and  increase  with  crack  length  for  given  val¬ 
ues  of  S,  W,  and  Viinp.  Equations  9  and  18  in  the  Appendix  can  be  used  to 
obtain  estimates  of  T0  and  /(|“1AX.  The  estimates  of  /f|'1AX  will  only  yield  ap¬ 
proximate  trends,  however,  because  two  coefficients  are  not  known  precisely 
in  Eq  18. 


Experimental  Verification 

To  verify  the  predictions  of  Eqs  I  and  2,  we  have  tested  several  steel  and 
aluminum  specimens.  Table  1  gives  a  matrix  of  the  specimens  tested  and  their 
dimensions  and  aspect  ratios.  The  specimens  had  a  blunt  notch  uf  0.5  mm  (20 
mil)  tip  radius  to  prevent  fracture  during  the  tests.  We  measured  the  stress 
intensity  history  with  a  single  element  strain  gage  placed  a  small  distance— 3 
to  5  mm  (0.1  to  0,2  in.)— away  from  and  directly  above  the  notch  tip.  The 
gage  was  oriented  to  measure  strains  perpendicular  to  the  crack.  The  stress 
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mental  results.  They  agree  well  with  the  analytical  predictions  (Eqs  1  and  2) 
and  demonstrate  that  the  interpretation  of  the  specimen  behavior  in  the  one- 


point-bend  test  and  the  assumptions  made  in  the  analysis  are  correct. 

Figure  4  represents  the  stress  intensity  history  for  seveial  steel  and  alumi¬ 
num  specimens  impacted  at  5.3  m/s.  First  we  observe  that  for  specimens  with 
high  S/W  values,  higher  modes  of  oscillation  also  significantly  affect  the 
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FIG.  4 — I continued }. 
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stress  intensity  history.  Beyond  a  value  of  S/W  =  7,  the  stress  intensity  his¬ 
tory  displays  large  oscillations  superimposed  on  the  fundamental  oscillation 
(Fig.  4b).  This  behavior  complicates  the  interpretation  of  fracture  tests.  Thus 
specimens  with  S/W  values  less  than  seven  are  preferable.  (We  have  not  in¬ 
vestigated  whether  a/W  ratios  higher  than  0.3  would  suppress  the  higher  os¬ 
cillations  in  geometries  with  S/W  >  7).  By  taking  into  account  this  limitation 
on  S/  W,  test  durations  T0  in  the  range  of  90  to  500  /rs  can  be  readily  achieved 
with  1PB  specimens  of  practical  size. 
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I e I  Specimen  5 
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FIG.  4 — (continued). 


Test  durations  measured  in  Fig.  4  and  predictions  from  Eq  9  in  the  Ap¬ 
pendix  are  compared  in  Table  2.  Estimates  are  given  only  for  specimens  for 
which  the  compliance  can  be  evaluated  with  reasonable  confidence  from 
handbook  tabulations.  The  agreement  between  predictions  and  measure¬ 
ments  is  sufficiently  good  to  confirm  that  the  first  mode  of  oscillation  of  the 
specimen  largely  controls  the  test  duration  for  specimens  with  S/W  <  7.  The 
largest  error  is  on  the  order  of  30% .  The  greatest  discrepancies  arise  for  those 
stress  intensity  histories  that  show  marked  inflections,  and  they  may  be  at¬ 
tributable  to  the  influence  of  higher  mode  oscillations. 
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TABLE  2 — Comparison  of  predicted  and  measured  1PBT  durations. 


Specimen 

Size 

Number 

Predicted 
Half  Period 

(r,/2). 

Measured 
Duration  1PBT 
(T0). 

To  -  1/2  T, 

S/ IV 

To 

1 

300 

284 

-6.0 

7.0 

2 

1123 

8S4 

—  31 .5 

14.0 

3a 

261 

304 

14.3 

3.8 

3b 

294 

300 

2.0 

3.8 

3c 

34S 

340 

-1.5 

3.8 

4 

508 

434 

-17.1 

S.6 

Figure  4  further  demonstrates  that  the  stress  intensity  histories  for  steel 
and  aluminum  specimens  are  remarkably  similar  except  for  the  amplitude. 
Again,  this  finding  is  expected  from  Eqs  1  and  2.  Because  the  ratio  p/E  is 
about  the  same  for  steel  and  aluminum,  specimens  with  the  same  dimensions 
will  have  the  same  fundamental  period.  Because  higher  oscillations  are  multi¬ 
ples  of  the  fundamental,  we  anticipate  that  steel  and  aluminum  specimens 
will  have  a  stress  intensity  history  with  the  same  global  shape.  According  to 
Eq  2,  K?ax  values  for  steel  and  aluminum  specimens  of  identical  geometry 
should  be  approximately  in  the  ratio  of  3  to  1 .  The  measured  ratios  for  speci¬ 
men  sizes  1  to  5  in  Table  1  are  3.2,  3.1 , 2.6,  3.1 ,  and  2.9  to  1 ,  respectively. 

Tests  on  the  same  specimen  at  several  impact  velocities  show  that  i> 
proportional  to  Vin,p,  as  predicted  by  Eq  2.  This  is  illustrated  in  Fig.  5,  which 
shows  the  stress  intensity  histories  for  steel  specimen  4  tested  at  3.5  and  5,.' 
m/s,  plotted  after  scaling  with  the  impact  velocity. 

The  effect  of  geometrical  scaling  is  represented  in  Fig.  6.  The  stress  inten 
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FIG.  5 — Normalized  stress  intensity  histories  Jor  specimen  4  impacted  at  3.5  and  5.3  m/s. 
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TIME/WtDTH,  ms/m 

FIG.  0 — Normalized  stress  intensity  histories  /or  two  geometrically  similar  specimens  of  differ¬ 
ent  sizes. 


sity  divided  by  the  square  root  of  the  specimen  width  has  been  plotted  as  a 
function  of  time  divided  by  the  specimen  width  for  two  similar  specimens  (5 
and  6).  The  agreement  of  the  two  curves  is  relatively  good  considering  that  the 
scaling  of  the  in-plane  dimensions,  in  particular  of  the  crack  length,  was  not 
exact  (see  Table  1). 

Figure  7  shows  the  influence  of  increasing  the  crack  length  on  the  stress 
intensity  history  when  all  other  dimensions  are  kept  constant.  Both  70/2  and 
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FIG.  7 — Stress  intensity  histories  Jar  specimen  3  with  three  difjerent  crack  lengths 
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K JHAX  increase  with  crack  length  as  predicted  by  the  analysis.  Note,  however, 
that  the  initial  portion  of  the  stress  intensity  histories  is  the  same  for  all  three 
crack  lengths.  This  result  occurs  because  in  the  early  part  of  the  loading,  the 
dynamic  stress  intensity  factor  depends  on  time  and  not  on  crack  length.  Only 
after  elastic  waves  have  traversed  the  crack  length  several  times  does  the  stress 
intensity  become  proportional  to  the  square  root  of  crack  length  1/2,7.?]. 

Finally,  multiple  tests  on  the  same  specimen  at  the  same  impact  velocity 
have  demonstrated  that  the  stress  intensity  histories  are  quite  reproducible 
for  large  specimens  (less  than  6%  deviation  on  For  small  specimens, 

the  scatter  increases  somewhat  (15%  deviation  on  A'^AX),  presumably  be¬ 
cause  small  specimens  are  more  sensitive  to  misalignment  errors. 


Dynamic  Fracture  Tests 

As  part  of  an  investigation  of  dynamic  crack  instability  criteria,  we  have 
used  the  IPBT  to  determine  the  dynamic  fracture  toughness  of  commercial 
aircraft-quality  4340  steel  at  three  different  strain  rates.  Specimens  corre¬ 
sponding  to  sizes  3  and  4  in  Table  1  were  cut  from  63.5  by  12.7  mm  (2.5  by  0.5 
in.)  bars,  and  quenched  and  tempered  to  a  Rockwell  C  hardness  of  50.  The 
specimens  were  notched,  fatigue  precracked  in  accordance  with  ASTM  Stan¬ 
dard  Test  Method  for  Plane-Strain  Fracture  Toughness  of  Metallic  Materials 
(E  399),  and  instrumented  with  a  strain  gage. 

Six  tests  were  performed  at  each  of  three  impact  velocities — 3.5,  5.3.  and 
50  m/s  (1 1.5,  17.4,  and  164  ft/s) — to  achieve  stress  intensity  rates  of  3.3  X 
10\5  X  10\  and  3  X  10h  MPa  in1  2  s 1  (3.0  X  IO\  4.6  X  10s,  and  2.7  X  10" 
ksi  in.|,J  s"1).  The  experiments  at  impact  velocities  of  3.5  and  5.3  m/s  were 
performed  on  the  modified  pendulum  tester  with  specimens  of  size  4  with 
a/W  —  0.3  and  0.38,  respectively.  To  achieve  the  50  m/s  impact  velocity,  we 
resorted  to  a  compressed  air  launcher  to  accelerate  a  steel  hammer  plate.  The 
test  arrangement  remained  otherwise  the  same  as  for  the  experiments  on  the 
pendulum  tester.  The  high-impact-velocity  specimens  were  of  size  3.  For 
comparison,  we  also  measured  the  static  fracture  toughness  in  triplicate  fol¬ 
lowing  the  standard  procedure  in  ASTM  E  399. 

For  the  dynamic  fracture  tests,  we  obtained  the  stress  intensity  from  the 
strain  measurements  using  the  static  calibration  method.  The  dynamic  initia¬ 
tion  fracture  toughness.  A|d,  was  determined  from  the  experimental  stress 
intensity  history  as  the  stress  intensity  value  at  the  point  where  crack  exten¬ 
sion  is  first  noticeable  on  the  record.  Crack  extension  manifests  itself  by  either 
a  clear  drop  in  the  measured  strain  or  by  a  gradual  deviation  from  an  essen¬ 
tially  linearly  rising  stress  intensity  curve. 

Figures  8  to  10  show  typical  stress  intensity  histories  measured  at  the  three 
loading  rates,  together  with  the  point  of  crack  initiation.  The  curves  in  Figs.  8 
to  10  represent  the  true  stress  intensity  history  only  up  to  the  point  of  crack 
initiation,  because  once  the  crack  propagates  the  stress  intensity  is  no  longer 
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directly  proportional  to  the  recorded  strain.  Figure  8  shows  that  for  the  tests 
at  3.3  X  lO*  MPa  m1  2  s_l,  the  crack  started  to  propagate  and  then  arrested 
after  extending  2  to  5  mm  (0,1  to  0.2  in.).  At  the  higher  rates,  the  specimens 
fractured  completely. 

At  low  impact  velocities,  crack  initiation  causes  a  clear  drop  in  the  mea¬ 
sured  strain,  and  the  point  at  which  K\a  is  evaluated  can  therefore  be  easily 
determined  (Fig.  8).  At  the  highest  impact  velocity,  the  strain  record  has  a 
smoother  appearance;  thus  the  point  of  crack  initiation  was  chosen  as  the 
point  of  separation  of  the  recorded  strain  history  from  a  straight  line  fitted  to 
the  main  portion  of  this  curve  (Fig.  9).  For  the  intermediate  impact  velocity, 
some  test  records  showed  a  significant  departure  from  linearity  before  a  sharp 
drop  in  strain;  others  only  showed  a  rapid  drop  in  strain.  Therefore,  for  the 
intermediate  impact  velocity,  either  one  or  the  other  definition  of  the  point  of 
crack  initiation  was  used,  as  appropriate  for  the  particular  test  record 
(Fig.  10). 

The  results  of  the  static  and  dynamic  fracture  experiments  are  summarized 
in  Fig.  11.  Although  the  scatter  in  data  increases  with  loading  rate,  the  aver¬ 
age  value  of  the  dynamic  fracture  toughness  K id,  58.5  MPa  m'  2  (53.5  ksi 
in.1'2),  is  essentially  the  same  for  the  three  rates,  and  it  is  only  slightly  lower 
than  the  static  K\c  value.  63.7  MPa  m1  2  (58  ksi  in.1  2). 

Discussion 

Results  of  the  Fracture  Tests 

The  dynamic  fracture  toughness  measurements  reported  here  demonstrate 
that  in  the  relatively  narrow  range  of  strain  rates  3.3  X  10'' to  3.0  X  lO^MPa 


FIG.  8— Stress  intensity  and  hammer  loud  histories  for  fracture  test  at  low  stress  intensity  rate. 
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TIME.ms 


FIG.  9— Stress  intensity  history  Jar  fracture  test  at  high  stress  intensity  rate. 


TIME,  us 

FIG  10— Stress  intensity  and  hammer  load  histories  for  fracture  test  at  intermediate  stress 
intensity  rate. 

m1'5  s'1  (3  X  10s  to  2.7  x  i0h  ksi  in.1'2  s-1),  K,d  for  4340  stee!  HRC  50  is  not 
sensitive  to  changes  in  loading  rate.  Furthermore,  comparison  of  the  dynamic 
fracture  toughness  value  58.5  MPa  m1  2  (53.3  ksi  in.1  2)  with  the  static  frac¬ 
ture  toughness  value  63.7  MPa  nt*'2  (58.0  ksi  in.1  2)  shows  that  the  overall 
strain-rate  sensitivity  of  the  fracture  toughness  is  only  mild.  If  it  is  assumed 
that  the  strain-rate  sensitivity  of  the  fracture  toughness  is  a  consequence  of 
the  strain-rate  sensitivity  of  the  flow  stress,  then  our  results  are  consistent 
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FIG.  II—  Summary  of  fracture  test  results  for  four  stress  intensity  rates. 


with  dynamic  flow  stress  data  for  4340  steel  |/4,/5).  These  data,  obtained  for 
lower  hardnesses  (HRC  40  and  HRC  30),  indicate  little  sirain-rate  sensitivity 
for  the  flow  stress  over  a  range  of  strain  rates  from  10' 4  s'1  to  104  s'1. 

However,  our  results  contradict  the  dynamic  fracture  data  of  Hontnia  et  al 
[16,17]  for  nominally  the  same  material  and  for  comparable  loading  rates. 
Homma  et  al  [16,17]  report  a  K\A  value  of  31.7  MPa  m1'2,  more  than  30% 
lower  than  the  static  toughness.  This  discrepancy  may  reflect  differences  in 
microfailure  modes.  Homma  et  al  1/6)  observed  a  transition  from  fully  dimple 
fracture  at  low  strain  rates  to  a  mixture  of  dimples  and  cleavage  facets  at  high 
strain  rates.  Scanning  electron  microscope  (SEM)  observations  did  not  reveal 
a  change  in  microfracture  mode  in  our  experiments.  Under  both  static  and 
dynamic  conditions,  fracture  occurred  by  what  appears  a  mixture  of  dimple, 
intergranular,  and  transgranular  failures. 

The  difference  between  the  results  of  Homma  et  al  and  the  present  results 
may  also  be  due  to  a  difference  in  the  definition  of  the  point  of  crack  initia¬ 
tion.  Homma  et  al  measure  Kti  after  20  to  50  pm  of  crack  extension,  whereas 
the  K ^  value  we  report  probably  corresponds  to  200  to  500  pm  of  crack  exten¬ 
sion. 


Significance  of  the  1PBT 

The  present  investigation  has  confirmed  the  suggestion  by  Kalthoff  et  al 
12,3]  that  the  1PBT  provides  an  attractive  test  configuration  for  dynamic  frac- 
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ture  experiments.  It  produces  s  smoothly  varying  stress  intensity  history  that 
can  be  easily  characterized  in  a  calibration  experiment.  Furthermore,  the 
loading  rate  and  the  pulse  amplitude  and  duration  can  be  readily  adjusted 
over  an  order  of  magnitude  range  for  each  parameter.  The  analytical  and 
experimental  results  given  in  this  paper  provide  the  information  necessary  to 
choose  the  appropriate  combination  of  specimen  dimensions  and  impact  ve¬ 
locity  to  achieve  desired  test  conditions. 

The  results  for  the  IPBT  expressed  in  Eqs  1  and  2  also  supplement  the 
results  obtained  for  the  impacted  3PBT  reported  by  Bbhme  [6|.  They  provide 
an  independent  analysis  of  the  free  specimen  oscillations  that  were  shown  to 
play  an  important  role  in  the  impacted  3PBT,  and  they  yield  the  same  conclu¬ 
sion  regarding  the  representation  of  these  dynamic  effects  by  master  curves 
that  depend  only  on  specimen  aspect  ratios. 

The  maximum  stress  intensity  amplitude  that  can  be  obtained  in  a  IPBT 
using  a  conventional  pendulum  impact  tester  is  rather  low,  particularly  in 
experiments  with  lighter  alloys  such  as  aluminum  or  titanium  alloys.  For 
these  materials,  devices  capable  of  higher  impact  velocities  are  required.  This 
problem  may  also  be  circumvented  by  attaching  ballast  plates  to  the  specimen 
extremities.  The  high-strain-rate  fracture  experiments  also  showed  that  the 
point  of  crack  initiation  is  not  always  clearly  defined  on  the  records  of  strain 
as  a  function  of  time.  Further  study  of  this  problem  is  desirable  to  obtain  an 
unambiguous  definition  of  the  point  of  dynamic  crack  initiation. 

An  interesting  aspect  of  the  IPBT  was  brought  to  light  by  the  fracture  ex¬ 
periments  in  which  crack  arrest  occurred.  They  indicated  that  the  one-point- 
bend  test  could  provide  a  convenient  experimental  procedure  to  study  not 
only  dynamic  crack  initiation  but  also  crack  propagation  and  arrest.  The  ad¬ 
vantage  of  the  one-point-bend  test  for  crack  propagation  investigations  would 
be  that  because  of  the  inertial  loading,  the  driving  force  on  the  running  crack 
tip  and  the  crack  velocity  could  be  easily  controlled  by  changing  the  impact 
velocity.  Further,  because  of  its  very  simple  boundary  conditions,  the  one- 
point-bend  test  should  be  much  easier  to  simulate  in  numerical  calculations 
than  the  conventional  wedge-loaded  compact  tension  crack  arrest  test  config¬ 
uration. 

Summary  and  Conclusions 

The  stress  intensity  history  obtained  by  impacting  an  unsupported, 
cracked,  bend  specimen  has  been  investigated.  A  combination  of  experi¬ 
ments  and  analysis  was  used  to  establish  the  influence  of  specimen  material 
and  dimensions  and  impact  velocity  on  the  maximum  stress  intensity  ampli¬ 
tude  and  the  test  duration.  Stress  intensity  pulses  of  durations  varying  be¬ 
tween  90  and  500  /is  can  be  readily  achieved.  For  steel,  maximum  amplitudes 
exceeding  100  MPa  m|,J  have  been  produced  with  a  pendulum  impact  tester. 

The  applicability  of  the  test  procedure  has  been  illustrated  by  measuring 
the  dynamic  fracture  toughness  of  4340  steel  (HRC  50)  at  three  different 
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strain  rates.  The  tests  indicated  little  strain-rate  sensitivity  for  the  fracture 
toughness  of  this  material. 

Finally,  it  is  suggested  that  the  orie-point-bend  test  could  also  conveniently 
serve  to  investigate  dynamic  crack  propagation  and  arrest. 
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APPENDIX 


Over  the  years  several  authors  have  presented  analyses  of  the  impacted-bend  speci¬ 
men  to  help  in  evaluating  dynamic  fracture  experiments  \t>,18-22\-  Here  wc  present 
simple  derivations  to  obtain  the  test  duration  To  and  the  maximum  stress  intensity 
amplitude  for  the  1PBT.  The  resulting  expressions  are  to  be  considered  more  as 
a  dimensional  analysis  than  as  precise  quantitative  predictions. 

Eitlmation  of  the  Duration  T, 

To  evaluate  the  test  duration  To,  we  assume  that  it  is  equal  to  one  half  the  funda¬ 
mental  period  of  oscillation  T\  of  the  cracked  beam.  Ireland  has  previously  proposed 
an  empirical  formula  to  estimate  the  period  of  inertial  oscillation  in  the  impacted 
3PBT  [2/].  We  estimate  T\  using  Rayleigh's  method.  To  do  so.  we  assume  for  the 
deflection  d  of  the  beam  the  following  function  of  position  x  along  the  axis: 


d(x)  =  d{ 


•(7*) 


0  <  x  <  — 
2 


The  first  term  on  the  right  of  Eq  3  represents  the  center  line  deflection  of  the  un¬ 
cracked  beam,  and  the  second  term  represents  an  estimate  of  the  deflection  due  to  the 
crack. 

Next  we  express  d^  and  in  terms  of  the  respective  elastic  compliances  CUN  and 
Ccr,  and  of  an  effective  dynamic  load  amplitude  Pmax  •  where  P^x  can  be  regarded 
as  the  maximum  hammer  load  corrected  to  account  for  the  inertia  forces  in  the  sense 
of  the  analysis  reported  in  Refs  19  and  20.  Thus 
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The  beam  deflection  at  any  point  x  and  time  r  is  given  by 
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Using  Eq  5,  we  calculate  the  maximum  kinetic  energy  of  the  specimen 
££!Tx  =  CJUN  |  +  ^  CUNCCB  +  {  C’n 

The  maximum  strain  energy  is  simply 

fST*  =  \  (^maxACun  +  CCR) 

Equating  kinetic  energy  and  strain  energy  yields 


Ts  =  2* 


IpBW- 


(^f 


un  ^  <CunCCr  Wr 


(Ccr  +  Cun) 


(6) 


(7) 


(8) 


Substituting 


S' 


Cmn  —  - r  1  +3(1  +  e)l— - 

UN  4 EBW3  \  S 


and  rearranging  yields 


where  Q  -  Ccr/Cun  and  Tm  =  (4  S2/»  W)sfip/E  is  the  fundamental  period  of  the 

uncracked  beam,  c 

The  dependence  of  Tt  on  specimen  dimensions  and  material  properties  given  by  bq 
9  is  the  same  as  the  dependence  anticipated  from  Ireland's  formula  1 21]. 


Eatlmatlon  of  the  Maximum  Amplitude  /(*** 

To  estimate  the  maximum  stress  intensity  amplitude,  we  express  it  as 


^MAX 


3  S  Vira 

7  ~dw2 


F{a/W)P 


FFF 

MAX 


(10) 


where  F(a/W)  is  the  handbook-tabulated  stress  intensity  calibration  appropriate  for 
the  given  S/W  ratio  |231-  We  then  proceed  to  express  FMAX  using  Newton  s  law  and 
conservation  of  momentum  and  energy. 

At  every  instant  r,  the  force  Ph  on  the  hammer  is  given  by 

Ph«)  =  c4)  U1) 

where  P***  is  defined  in  Fig.  2  andgfl/Ti)  is  a  nondimensional  function  that  can  be 
obtained  experimentally. 
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Integration  of  Newton's  law  for  the  hammer  from  the  initial  time  of  impact  t  -  0  to 
the  time  tMAX  at  which  the  specimen  deflection,  and  hence  the  stress  intensity,  is  maxi- 

mum  gives 


(o  P^Xg[jr)dl  =  TtPZAX\ G  (^)  -  G(0> 


(12) 


=  7’|/>JJAXu  =  M(V,mp  -  V„) 


where  M  is  the  hammer  mass.  Vllr,p  is  the  initial  impact  velocity,  VH  is  the  residual 
hammer  velocity  at  time  rMAX.  and  a  is  a  coefficient  resulting  from  the  integration  of 
g{t/T, ).  Considering  the  approximately  sinusoidal  shape  of  the  hammer  load  in  Fig.  2 
(period  »  T0  *  T,/2).  we  see  that  a  is  on  the  order  ot  =  0.15. 

From  Eq  12  we  obtain 


V„  =  V,„ 


-j"  u  M  A  X 

aT,PH 


M 


(15) 


Application  of  momentum  conservation  to  the  hammer-specimen  system,  with  sub 
stitution  of  Eq  13  for  VH,  yields 


Vs  = 


aT,P^x 


114) 


where  m  is  the  specimen  mass  and  Vs  is  the  velocity  of  the  center  of  mass  of  the  speci¬ 
men.  which  at  iMAX  is  also  the  uniform  velocity  of  the  whole  specimen. 

Finally,  the  energy  conservation  equation  can  be  written 


y  vLf M  =  y  VJ„M  +  y  mV*  +  -^  >2«\,n  +  C< „)  (15) 


where  the  last  term  of  the  right-hand  side  represents  (he  strain  energy.  P^\x-  Ci  s. 
and  CCR  have  been  defined  above. 


We  now  write  f^wx  as  a  tnultiple  of  P\[ 


pti-f  -  _  pM" 

^  ‘  h 


(16) 


The  results  of  Ref  20  indicate  that  I  /A  should  be  between  0  and  I 
By  substitution  of  Eqs  13,  14,  and  16  in  Eq  15  we  obtain 
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Finally,  combining  Eqs  9.  10,  and  17.  we  obtain 

*rX  =  V,mpV™  s[\2pE  F(a/w)^  j  Jl  +  3(1  +  v)(~)  LW) 


Equations  9  and  18  are  the  expanded  versions  of  Eqs  1  and  2.  The  assumption  To  = 
T \/2  has  been  made  in  writing  Eq  I. 
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One-Point  Bend  Test 

By  Jacques  Henri  Joseph  Giovanola 
Research  Engineer 
SRI  International 


INVESTIGATIONS  of  the  dynamic  be¬ 
havior  of  the  impacted  three-point  bend  spec¬ 
imen  have  shown  that  immediately  after 
impact  by  the  hammer,  the  ends  of  the  spec¬ 
imen  lose  contact  with  the  supports  (Ref  71). 
Consequently,  the  initial  part  of  the  response 
curve  discussed  in  the  preceding  section  is 
strictly  the  result  of  inertial  loading,  without 
contribution  from  the  reactions  at  the  sup¬ 
ports.  Therefore,  if  high  impact  velocities  are 
applied  to  achieve  high  loading  rates  and 
short  times-to-fracture.  it  becomes  unneces¬ 
sary  to  support  the  bend  specimen  during 
impact,  and  a  so-called  one-point  bend  test 
can  be  performed  (Ref  72-74). 

The  one-point  bend  test  allows  reliable 
toughness  data  to  be  obtained  at  loading  rates 
one  to  two  orders  of  magnitude  higher  than 
the  maximum  loading  rate  considered  in 
ASTM  E  399  for  fracture  tests  (Ref  l).  Sev¬ 
eral  tests  performed  with  different  specimen 
sizes  and  impact  velocities  on  high-strength 
4340  steel  have  demonstrated  the  usefulness 
of  the  one-point  bend  test  and  have  provided 
dynamic  fracture  toughness  values  in  good 
agreement  with  values  obtained  with  more 
complicated  tests. 

Test  Principle 

The  one-point  bend  test  uses  a  single-edge 
cracked  specimen  and  the  same  testing  ar¬ 
rangement  as  a  conventional  three-point 
bend  test,  except  that  the  end  supports  are 
removed.  The  procedure  is  illustrated  in  Fig. 
13.  When  the  hammer  strikes  the  specimen, 
the  center  portion  of  the  specimen  is  acceler¬ 
ated  away  from  the  hammer;  the  end  portions 


of  the  specimen  lag  behind  by  an  amount  d0- 
because  of  inertia  This  causes  the  specimen 
to  bend  and  to  load  the  crack  tip. 

A  typical  stress-intensity  history  resulting 
from  such  inertial  loading  is  shown  in  Fig. 
14.  The  impact  conditions  for  the  test  were 
such  that  no  crack  extension  occurred.  The 
stress-intensity  history  is  approximately  si¬ 
nusoidal;  it  increases  monotonically  up  to  the 
maximum  amplitude  and  then  decreases 
monotonically.  Later  oscillations  in  the 
stress-intensity  history  are  of  much  smaller 
amplitude.  For  comparison,  the  hammer  load 
is  also  shown  in  Fig  14.  It  is  apparent  that 
there  is  no  direct  relationship  between  the 
hammer  load  history  and  the  crack  tip  stress- 
intensity  history  in  the  one-point  bend  test. 

The  stress-intensity  history  obtained  in  the 
one-point  bend  test  depends  only  on  the  im¬ 
pact  velocity  and  on  the  specimen  geometry, 
size,  and  material  It  is  independent  of  the 
fracture  properties  of  the  material  tested. 
Therefore,  the  concepts  of  response  curve 
and  time-to- fracture  described  earlier  in  this 
article  can  also  be  applied  to  the  onc-point 
bend  test  when  used  for  dynamic  fracture 
toughness  measurements  in  the  small-scale 
yielding  regime. 

As  described  previously,  dynamic  fracture 
measurements  involve  two  steps:  determina¬ 
tion  of  the  response  curve,  i.e.,  the  stress- 
intensity  history,  for  a  given  specimen  geom¬ 
etry  and  class  of  material  and  measurement 
of  the  time-to-fracture  with  a  fatigue  pre¬ 
cracked  specimen.  The  dynamic  fracture 
toughness  is  then  evaluated  using  Eq  23. 

Alternatively,  the  stress-intensity  history 
and  the  stress-intensity  level  at  which  crack 
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initiation  occurs  can  be  measured  for  each 
specimen  by  measuring  the  crack  tip  strain 
intensification  with  a  calibrated  strain  gage 
placed  near  the  crack  tip.  The  measured 
strain  histories  can  be  related  to  the  stress- 
intensity  histories  by  two  methods.  In  the 
first  method,  the  strain  gage  signal  it  cali¬ 
brated  against  a  known  applied  stress  inten¬ 
sity  in  a  static  three-point  bend  experiment.  It 
is  then  assumed  that  the  static  calibration  also 
applies  in  the  dynamic  case.  In  the  second 


method,  the  strain  history  «,(»)  is  related 
directly  to  the  stress-intensity  history  by 
means  of  the  elastic  singularity  solution.  The 
stress  intensity  is  then  given  by: 


KM  * 


r,  .  e 

30  ( 
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36\  I 

i  +  si.ir 

sm~-  v 

1  -sinr 
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L  2 

2  \ 
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2/J 

(Eq  25) 


where  rand  6  are  the  polar  coordinates  of  the 
center  of  the  strain  gage  with  ongin  at  the 
crack  tip,  and  E  and  v  are  Young’s  modulus 
and  Poisson's  ratio,  respectively  Both  meth¬ 
ods  work  equally  well  if  the  strain  gage  is 
positioned  properly  This  aspect  is  discussed 
in  a  subsequent  section. 

Dependence  of  Stress-Intensity  His¬ 
tory  on  Test  Parameters.  As  illustrated  in 
Fig.  14,  the  stress-intensity  history  is  charac¬ 
terized  by  the  duration  of  loading,  T0.  and 
the  maximum  stress-intensity  amplitude, 
K?“.  It  has  been  shown  (Ref  74)  that  T„  and 
#f7“*  can  be  expressed  as: 

(Eq  27) 

where  W,  L,  and  a  are  the  width,  length,  and 
crack  depth  of  the  specimen,  respectively:  p, 
E,  and  v  are  the  mass  density.  Young’s  mod¬ 
ulus,  and  Poisson’s  ratio,  respectively;  Vlm(, 
is  the  impact  velocity;  and  6,  and  h2  are  non- 
dimensional  functions. 

Equations  26  and  27  have  been  verified 
experimentally  (Ref  74)  and  provide  the  scal¬ 
ing  rules  for  the  one-point  bend  test  stress 
intensity  history  Given  the  stress-intensity 
history  for  one  material,  one  specimen  size 
and  geometry,  and  one  impact  velocity,  Eq 
26  and  27  allow  construction  of  the  suess- 
imensity  history  for  any  material,  impact  ve¬ 
locity.  and  specimen  size  as  long  as  geomet¬ 
rical  similarity  of  the  specimen  is  preserved 
In  particular,  for  given  specimen  dimensions 
and  material,  the  one-point  bend  test  stress- 
intensity  history  is  proportional  to  the  impact 
velocity  and  is  independent  of  specimen 
thickness.  For  a  specimen  and  impact  veloc¬ 
ity,  it  can  further  be  shown  that  both  *7“* 
and  T„  increase  with  an  increase  in  crack 
length. 

The  experimental  one-point  bend  lest 
stress-intensity  histories  for  several  specimen 
geometries,  normalized  according  to  Eq  26 
and  27,  are  shown  in  Fig.  15  and  16.  These 
curves  were  obtained  from  steel  specimens 
with  blunt  notches.  Figure  15  illustrates  the 
effect  of  changing  the  length-over-width  ra¬ 
tio  Figure  16  presents  the  normalized  stress- 
intensity  curves  for  a  specimen  with  the  stan¬ 
dard  ASTM  E  399  (Ref  I)  proportions  (U 
W  =  4.2)  and  for  four  ratios  of  crack  length 
to  specimen  width  (alW  =  0.3, 0.4, 0.5.  and 
0.6). 

Advantages  and 
Limitations 

The  primary  advantage  of  the  one-point 
bend  lest  is  that  the  measured  stress-intensity 
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Fig.  IS  Normalized 
•tress-intensity  histories  for  two 
specimens  with  different  LIW 
ratios 
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Fig.  16  Normalized 
stress-intensity  histories  for 
ASTM  standard  bend  specimen 
with  four  different  crack 
lengths 


Fig.  13  Test  arrangement  and  loading  mechanism  for  the  one-point 
bend  test 


Fig.  14  Typical  stress-intensity  and  load  histories  obtained  in  the 
one-point  bend  tost 
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history  incorporates  dynamic  effects  com¬ 
pletely.  Therefore,  no  limits  need  to  be 
imposed  on  the  impact  velocity  and  the  test 
duration  to  fracture. 

Combined  use  of  the  one-point  bend  test 
and  the  impact  response  curve  concept 
produces  a  simple,  smoothly  varying  re¬ 
sponse  curve  For  a  given  specimen  geome¬ 
try,  the  response  curve  must  be  measured 
experimentally  only  once.  It  can  then  be  ex¬ 
tended  to  other  specimen  sizes  (of  same 
aspect  ratios),  impact  velocities,  and  materi¬ 
als  by  using  the  scaling  rules  of  Eq  26  and 
27. 

Another  advantage  of  the  one-point  bend 
test  is  that  the  test  duration,  T„,  and  the  stress 
intensification  rate  can  be  varied  over  an 
order  of  magnitude  range  (80  to  800  p.s  and 
3  x  105  to  3  x  I06  MPaVm-s'1,  or 
2.7  x  I05  to  2.7  x  106  ksiVTrT  •  s_  1 )  by 
changing  the  specimen  geometry  and  the  im¬ 
pact  velocity.  Because  "loading  occurs  by 
inertia  only,  the  one-point  bend  test  response 
curve  is  insensitive  to  the  characteristics  of 
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the  loading  machine,  provided  the  hammer 
mass  is  much  larger  than  the  specimen 
mass. 

When  measurement  of  the  stress-intensity 
history  is  direct,  the  test  procedure  to  obtain 
dynamic  fracture  toughness  using  the  one- 
point  bend  test  is  also  simple.  It  requires  only 
limited  instrumentation,  and  interpretation  of 
the  test  data  is  straightforward. 

Use  of  the  one-point  bend  test  currently  is 
restricted  to  small-scale  yielding  conditions. 
Furthermore,  the  maximum  stress-intensity 
amplitude  that  can  be  obtained  using  conve¬ 
nient  specimen  sizes  and  a  conventional  pen¬ 
dulum  impact  tester  is  rather  low.  particu¬ 
larly  for  lighter  alloys  such  as  aluminum  and 
titanium  alloys.  For  these  materials,  devices 
capable  of  higher  impact  velocities  are  re¬ 
quired.  This  problem  can  also  be  circum¬ 
vented  by  attaching  ballast  plates  to  the  ends 
of  the  specimen. 

Test  Setup  and 
Instrumentation 

Figure  17  illustrates  the  typical  test  setup 
for  performing  a  one-point  bend  test  on  a 
pendulum  impact  tester.  The  specimen 
holder  used  in  a  Charpy  or  lzod  test  is 
replaced  by  a  simple  frame  that  supports  the 
specimen,  while  allowing  it  to  move  freely  in 
the  horizontal  plane.  Depending  on  the  de¬ 
sign  of  the  original  pendulum  and  hammer, 
the  impact  tester  may  require  retrofitting  with 
a  new  hammer  and  striker  that  will  not  inter¬ 
fere  with  the  specimen  edges  or  the  support 
frame.  The  sinker  should  be  instrumented 
with  two  foil  gages— one  on  each  lateral 
side — to  record  the  impact  load  on  the 
hammer. 

In  the  response  curve  time-to-fracture  ap¬ 
proach,  the  time-to-fracture  is  defined  as  the 
interval  between  the  time  of  impact  and  the 
time  at  which  crack  extension  occurs.  The 
time  of  impact  can  be  determined  conve¬ 
niently  from  the  leading  edge  of  the  hammer 
load  signal,  whereas  the  time  of  crack  exten¬ 
sion  can  be  measured  by  an  uncalibrated 
strain  gage  placed  near  the  crack  tip,  or  by  a 
magnetic  pickup  .  If  the  re¬ 

sponse  curve  is  to  be  used  later  in  conjunc¬ 
tion  with  other  impact  velocities,  it  is  also 
necessary  to  measure  the  hammer  impact 
velocity. 

When  calibrated  strain  gages  are  used  to 
perform  a  direct  measurement  of  the  stress- 
intensity  history,  the  gage  element  should  be 
approximately  2  by  2  mm  (0.08  by  0  08  in.) 
or  smaller  The  center  of  the  gage  should  be 
positioned  3  to 4  mm  (0. 12  toO  16  in.)  above 
the  crack  plane,  on  a  line  normal  to  the  crack 
path  and  passing  through  the  crack  tip 

Two  amplifiers  with  a  frequency  response 
greater  than  100  kHz  are  required  to  process 
the  hammer  and  specimen  gage  signal  for 


recording.  These  signals  are  best  recorded  on 
a  digital  storage  oscilloscope  with  a  time  res¬ 
olution  of  I  ps  or  better.  If  this  recording 
instrument  has  the  capability  to  store  signals 
that  occur  before  the  trigger  time,  the  ham¬ 
mer  or  specimen  gage  signals  can  be  used  as 
•  trigger.  Otherwise,  •  trigger  source  must  be 
incorporated  into  the  test  setup.  A  laser  beam 
interrupted  by  the  hammer  in  conjunction 
with  a  photodiode  is  one  option.  Another  is  a 
simple  electrical  circuit  closed  by  the  contact 
of  hammer  and  specimen;  however,  this 
method  may  cause  transient  noise  on  the 
strain  gage  signals. 

Specimen  Size  and  Geometry.  To 

guarantee  small-scale  yielding  conditions, 
specimens  used  in  the  one-point  bend  test  are 
subject  to  the  same  size  requirements  as 
specimens  used  for  static  fracture  measure¬ 
ments  according  to  ASTM  E  399  (Ref  I).  In 
addition,  restrictions  on  the  value  of  the 
length-to-width  ratio,  L/W.  must  be  imposed 
because  of  the  natural  modes  of  vibration  of 
the  specimen.  It  has  been  shown  experimen¬ 
tally  that  for  large  values  of  LtW,  higher 
mode  frequencies  are  excited  when  the  spec¬ 
imen  is  impacted,  thus  causing  large  oscilla¬ 
tions  in  the  stress-intensity  history  (Ref  74). 
To  avoid  this  undesirable  tendency,  values  of 
LiW  shoulJ  not  exceed  7. 

In  principle,  there  is  no  limit  to  the  abso¬ 
lute  size  of  the  specimens  beyond  those  dic¬ 
tated  by  the  development  of  plasticity.  In 
practice,  however,  small  samples  are  more 
susceptible  to  misalignment  errors  than  large 
ones.  For  a  given  impact  velocity,  small 
samples  show  more  test-to-test  variability  in 
the  maximum  amplitude  of  the  stress-inten¬ 
sity  history  Consequently,  samples  less  than 
100  mm  (4  in.)  in  length  are  not  recom¬ 
mended. 

Determination  of  Dynamic  Tough¬ 
ness  From  Strain  Gage  Records.  As 

mentioned  previously,  there  are  two  methods 
of  determining  dynamic  toughness  from 
strain  gage  records:  the  static  calibration 
method  and  the  direct  method  using  Eq  25 
The  direct  method  is  simpler,  but  requires 
care  in  positioning  the  strain  gage  to  measure 
r  and  0.  The  gage  should  be  far  enough  from 
the  crack  tip  to  avoid  the  plastic  zone,  but 
close  enough  to  the  tip  to  remain  in  the  zone 
dominated  by  the  singular  term  of  the  elastic 
solution. 

Test  results  show  that  the  recommended 
distance  of  3  to  4  mm  (0  12  to  0  16  in.) 
meets  both  requirements,  even  for  the  small¬ 
est  recommended  specimen  size.  Further¬ 
more,  the  crack  front  obtained  from  fatigue 
precracking  invariably  exhibits  a  slight  cur¬ 
vature  Thus,  when  determining  r  and  0,  it 
must  be  determined  whether  the  origin  of  the 
coordinate  system  should  be  taken  at  the 
location  of  a  virtual  average-  crack  front  (fol¬ 
lowing  the  averaging  procedure  detailed  in 


Fig.  17  Experimental  setup  used 
to  perform  the  one-point  bend 
test 


Velocity 

measurement 


Pendulum 


ASTM  E  399),  or  at  the  actual  position  of  the 
crack  tip  on  the  side  of  the  specimen  where 
the  strain  gage  is  mounted  Although  further 
verification  is  required,  experience  shows 
that  measuring  r  and  0  from  the  actual  crack 
tip  on  the  gage  side  of  the  specimen  provides 
more  consistent  results.  This  practice  is 
therefore  recommended. 

Another  concern  is  the  determination  of 
the  point  of  crack  initiation  when  the  test  is 
performed  with  strain  gages.  For  experi 
ments  with  long  times-to  fracture  ( >50  |is), 
crack  initiation  causes  a  clear  drop  in  the 
measured  strain,  and  the  point  at  which  Ku  is 
evaluated  can  therefore  be  determined  unam¬ 
biguously  (Fig  18).  For  experiments  with 
short  times-to-fracture  (<50  |xs).  experimen¬ 
tal  results  and  simulations  indicate  that  the 
strain  recorded  by  the  strain  gage  may  con¬ 
tinue  to  increase  once  the  crack  has  begun  to 
extend  This  may  result  in  overestimates  of 
the  dynamic  fracture  toughness. 

To  determine  the  correct  point  of  crack  ini¬ 
tiation.  the  procedure  illustrated  in  Fig  19 
should  be  followed.  The  reference  stress- 
intensity  history  for  the  specimen — i  e..  the 
complete  history  untruncated  by  crack  exten¬ 
sion — is  first  scaled  to  the  appropriate  test 
conditions  The  scaled  reference  stress-inten¬ 
sity  history  is  then  superimposed  or  the 
curve  recorded  during  the  fracture  test.  The 
point  where  the  two  curves  begin  separating 
is  then  selected  as  the  point  of  crack  ini¬ 
tiation. 

If  the  reference  stress-intensity  history  is 
not  available,  a  tangent  straight  line  can  be 
fitted  to  the  fracture  test  record.  As  before 
the  point  of  separation  of  the  experimental 
record  from  the  tangent  line  determines  crack 
initiation. 

Performing  Tests  With  Ballost  at  the 
Specimen  Ends.  When  testing  low-density 
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fig.  18  Stress-intensity  and  hammer-load  histories  for  fracture  test 
with  a  long  time-to-fracture 
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materials  such  as  aluminum  or  titanium  al¬ 
loys,  the  inertial  loading  capacity  of  the  one- 
point  bend  test  may  be  insufficient  to  initiate 
crack  propagation.  To  circumvent  this  limi¬ 
tation,  ballast  may  be  affixed  to  the  specimen 
ends,  as  illustrated  in  Fig.  20.  The  stress- 
intensity  history  for  an  aluminum  specimen 
outfitted  with  ballast  plates  is  illustrated  in 
Fig.  21  along  with  the  stress-intensity  history 
for  the  same  specimen  without  ballast  plates. 
The  tests  were  performed  at  an  impact  veloc¬ 
ity  of  5.3  m/s  (17.4  ft/s). 

The  dimensions  of  the  specimen  are 
L  -  241  mm  (9  5  in.).  W  =  63.5  mm  (2.5 
in.),  B  =  12.7  mm  (0.5  in  ),  and  a  =  19 
mm  (0.75  in  ).  Two  ballast  plates  of  steel  (89 
by  25.4  by  6.4  mm,  or  3.5  by  i  .O  by  0.25 
in.)  were  attached  on  each  side  of  the  speci¬ 
men.  The  maximum  stress-intensity  ampli¬ 
tude  is  significantly  higher  for  the  specimen 
with  ballast. 


Fig.  19  Stress-intensity  history  for  a  fraeture  test  with  a  short 
time-to-fraeture  and  procedure  for  determining  Kia 


is 

M 

$ 

e» 

1 

2 

<75 


B-6 


*•  '  -  - 


a;\Vv>; 


•  s  "  -  ■  W  "  a 


\  V 

*  N  *  *  *  i 


r.-x. 


.  i . "  ^ 


V-Vn'  '/•  ‘O  ^  ,-_- 

|  276/High  Strain  Rate  Tasting 


'J*  v  ~mu  -+  :\r  m~.tr  w 


S 

£ 

i 


i", 

!\- 

■  » 

V 

fc: 


P!g-  20  Specimen  and  ballatt  plates  used  to  obtain  response 
curve  in  Fig.  21 


Pig.  21  Stress-intensity  histories 
for  an  aluminum  specimen  with 
and  without  ballast  plates 
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Appendix  C 


DETERMINATION  OF  DYNAMIC  STRESS  INTENSITY  HISTORIES  BY 
STRAIN  GAGE  MEASUREMENTS 

INTRODUCTION 

The  major  difficulty  in  measuring  the  mode  I  dynamic  initiation 
fracture  toughness  is  the  lack  of  a  direct  relationship  between  the 
stress  Intensity  history  at  the  crack  tip  and  the  history  of  the  load 
applied  away  from  the  crack  tip.  Therefore,  measurements  must  be 
performed  in  the  immediate  vicinity  of  the  crack  tip  to  obtain  the  stress 
intensity  history.  Optical  techniques,  such  as  photoelasticity *  and  the 

o 

shadow  optical  method  of  caustic  have  been  successfully  used  for  these 
measurements.  Unfortunately,  these  techniques  are  complex,  and  hence  not 
suitable  for  routine  testing. 

Simpler,  yet  reliable  techniques,  are  required  for  dynamic 
initiation  toughness  measurements.  One  such  technique  uses  strain  gages 
to  measure  strain  histories  near  the  crack  tip.  The  stress  intensity 
history  is  then  related  to  the  strain  histories  using  either  of  two 
methods:  the  singularity  solution  method  or  the  static  calibration 

method. 

In  this  appendix,  we  first  review  the  principle  of  the  strain  gage 
technique  and  discuss  its  limitations.  Then  we  use  experimental  results 
to  evaluate  the  effect  of  these  limitations  on  the  measurements  of  stress 
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intensity  histories  for  pure  mode  I  and  mixed  mode  I  and  mode  II  dynamic 
loading.  Finally,  we  discuss  the  parameters  affecting  the  strain  history 
records  once  the  crack  has  started  to  extend,  with  the  purpose  of 
defining  crack  initiation  and  to  estimate  the  crack  velocity  and  the 
propagation  toughness  on  the  basis  of  the  strain  gage  records. 

PRINCIPLE  AND  LIMITATIONS  OF  THE  STRAIN  GAGE  TECHNIQUE 

The  strain  gage  technique  for  determining  the  stress  intensity 
history  is  based  on  the  assumption  that,  up  to  the  point  of  crack 
initiation,  the  strain  at  a  given  location  in  the  vicinity  of  the  crack 
tip  is  directly  proportional  to  the  stress  intensity.  Two  methods  for 
relating  the  strain  to  the  stress  intensity  are  discussed  below:  the 
singularity  solution  method  and  the  static  calibration  method. 

Singularity  Solution  Method 

The  singularity  solution  method  relies  on  the  analytical 
elastodynamlc  solution  for  the  stress  field  around  a  crack  tip  to  relate 
the  strain  to  the  stress  intendii  iactor.  Referring  to  Figure  C-l  for 
the  definition  of  the  coordinate  systems,  the  stresses  are  given  in  terms 
of  a  series  expansion  in  the  polar  coordinates  v  and  0  centered  at  the 
crack  tips  as 
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where  K^(t)  is  the  time-dependent  mode  I  stress  intensity  factor  and 
A^(t)  are  the  time-dependent  coefficients  of  the  higher  order  terms  of 
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the  expansion.  Equations  (C-L)  and  (C-2)  have  the  same  form  as  the 
elastostatic  solution,  except  that  the  coefficients  are  now  time 
dependent.  In  particular,  the  angular  functions,  f^(6)  and  g^(9)  ,  are 
the  same  as  in  the  static  case.  Assuming  plane  stress  conditions  at  the 
specimen  surface,  the  strain  measured  by  the  strain  gage,  positioned  as 
shown  in  Figure  C-l,  can  be  expressed  as 

e  =  -i—  (a  -  v  a  )  (C-3) 

y  E  '  y  x  v  ' 

Inserting  (C-l)  and  (C-2)  into  (C-3)  yields 
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Very  near  the  crack,  tip,  the  first  (singular)  term  in  equation  (C-4) 
dominates  over  the  other  terms,  so  the  strain  Is  approximately  given  by 
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Knowing  the  coordinates  r  and  0  of  the  center  of  the  gage  with  respect  to 
the  initial  crack  tip,  equation  (C-5)  can  be  used  to  calculate  the  stress 
intensity  history  from  a  single  strain  history  record  taken  during  the 
dynamic  fracture  test.  In  dynamic  fracture  tests,  the  strain  gage  is 
usually  located  as  shown  in  Figure  C-l  with  the  gage  element  oriented  in 
the  y-direction.  The  advantages  of  this  configuration  will  be  discussed 
below. 

The  singularity  solution  method  introduces  four  sources  of  error  in 
the  measurement  of  the  scress  intensity  history.  First,  the  singular 
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term  in  equation  (C-4)  dominates  only  very  close  to  the  crack  tip.  For 
practical  reasons,  however,  the  strain  gage  has  to  be  positioned  a  finite 
distance  from  the  crack  tip,  if  only  to  avoid  the  plastic  zone.  As  the 
distance  from  the  crack  tip  is  increased,  the  contribution  of  the  higher 
order  terms  to  the  strain  may  become  important.  Thus  the  use  of  equation 
(C-5)  may  Introduce  a  nonnegligible  error  in  the  estimate  of  the  stress 
intensity  history. 

A  second  source  of  error  is  introduced  by  the  finite  size  of  the 
strain  gage  active  grid.  This  causes  the  strain  gage  to  measure  a  strain 
averaged  over  an  area  rather  than  the  strain  at  a  point,  as  assumed  when 
using  equation  (C-5).  In  the  presence  of  a  strong  strain  gradient,  such 
as  near  the  crack  tip,  the  average  strain  measured  by  the  gage  could  be 
significantly  different  from  the  strain  at  the  nominal  point  of 
measurement  (the  center  of  the  gage)  and  this  difference  could  introduce 
an  error  in  the  estimate  of  the  stress  intensity. 

A  third  type  of  error  arises  if  the  position  of  the  gage  with 
respect  to  the  crack  tip  is  not  measured  precisely.  In  particular,  when 
the  crack  front  is  slightly  curved,  there  is  a  question  of  whether  the 
polar  coordinate  system  should  be  defined  by  the  average  crack  length  or 
by  the  length  of  the  crack  on  the  face  where  the  strain  gage  is 
mounted.  Also,  small  angular  misalignments  of  the  gage  element  with 
respect  to  the  y-d.trection  may  cause  some  error  when  using  equation  (C-5) 
to  obtain  the  stress  intensity  factor. 

Estimates  of  the  importance  of  these  four  error  sources  are 
discussed  in  a  later  section. 


Static  Calibration  Method 

In  the  static  calibration  method,  first  proposed  by  Loss  ,  the 
proportionality  factor  between  strain  and  stress  intensity  is  determined 
experimentally.  The  specimens,  with  the  strain  gage  mounted  at  the  crack 
tip,  are  loaded  quasi-statically  using  a  test  configuration  for  which  the 
stress  Intensity  calibration  is  available  in  handbooks.  During  this 
calibration  test,  the  applied  load  and  the  strain  gage  output  are 
recorded  simultaneously.  The  proportionality  factor  between  strain  and 
stress  intensity  factor  is  calculated  from  these  two  records  using  the 
handbook  calibration  function. 

It  is  then  assumed  that  the  relationship  between  strain  gage  output 
and  stress  intensity  factor  obtained  under  quasi-static  conditions  is 
also  valid  under  dynamic  conditions  and  can  be  used  to  relate  the  strain 
history  to  the  stress  intensity  history.  This  assumption  implies  that 
the  time-dependent  coefficients  of  the  higher  order  terms  in  equation  C-4 
all  vary  in  proportion  to  the  stress  intensity  history,  and  that  the 
proportionality  factors  are  the  same  as  for  the  static  case.  However, 
this  may  not  be  true,  particularly  in  experiments  where  stress  waves  play 
an  important  role  and  interact  both  with  the  boundaries  and  the  crack 
tip. 

Thus,  similarly  to  the  singularity  solution  method,  the  static 
calibration  method  could  Introduce  an  error  related  to  the  way  higher 
order  terms  are  taken  into  account.  The  static  calibration  method, 
however,  avoids  the  potential  error  associated  with  the  strain 
integration  over  a  finite  area  and  with  the  positioning  of  the  strain 
gage  because  these  effects  are  automatically  accounted  for  in  the 
calibration  procedure.  From  this  point  of  view,  the  static  calibration 
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method  appears  to  be  more  reliable  than  the  singularity  solution, 
although  it  has  the  disadvantage  of  being  more  complicated  to  implement. 

USE  OF  THE  STRAIN  GAGE  METHOD  IN  MIXED-MODE  DYNAMIC  FRACTURE 

The  singularity  solution  method  can  be  applied  to  mixed-mode  dynamic 
fracture  experiments.  In  the  case  of  mixed  mode  I  and  II  loading,  the 
stress  field  is  given  by 
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where  K.^(t)  the  dynamic  mode  II  stress  Intensity  factor. 

As  in  the  pure  mode  I  case,  the  strain  at  any  location  and  in  any 
direction  can  be  expressed  In  terms  of  equation  C-7.  In  the  Immediate 
vicinity  of  the  crack  tip,  the  singular  terms  will  dominate  and  the 
strain,  say,  in  the  y-direction,  can  be  written  as 
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Therefore,  given  two  independent  strain  history  records,  the  mode  I  and 
mode  II  stress  intensity  factors  can  be  determined  at  each  instant  by 
solving  two  algebraic  equations  in  the  form  of  equation  C-8.  If 


additional  strain  gages  are  introduced,  higher  order  coefficients  can 
also  be  determined  using  the  strain  gage  method. 

This  multiple  strain  gage  approach  has  been  used  to  measure  the  mode 
I  and  mode  II  stress  intensity  histories  in  mixed  mode  fracture 
experiments  and  to  evaluate  the  first  higher  order  term  in  the  series 
expansion,  A^. 

ESTIMATION  OF  ERRORS  IN  THE  DETERMINATION  OF  STRESS  INTENSITY  HISTORIES 
USING  STRAIN  GAGE  MEASUREMENTS 

Comparison  of  Singularity  Solution  and  Static  Calibration  Methods 

A  first  approach  to  evaluate  the  reliability  of  the  strain  gage 
method  of  measuring  the  dynamic  stress  intensity  factor  is  to  compare  the 
estimates  obtained  with  the  singularity  solution  and  with  the  static 
calibration  methods  for  the  same  experiment.  In  doing  so,  we  assume  that 
the  static  calibration  method  will  provide  a  more  accurate  estimate  of 
the  stress  intensity  history  than  the  singularity  solution  method  for  the 
reasons  discussed  in  the  preceding  section.  Hence,  the  static 
calibration  method  will  provide  a  reference  with  which  to  compare  the 
singularity  solution  results. 

We  have  compared  estimates  of  the  stress  intensity  based  on  the 
singularity  solution  method  and  estimates  based  on  the  static  calibration 
method  for  several  specimen  geometries  impacted  in  1PBT  experiments. 

Some  specimens  had  a  fatigue  precrack,  others  contained  a  blunt  notch 
with  a  notch  root  radius  of  approximately  0.5  mm.  All  the  strain 
measurements  were  performed  with  the  same  type  of  strain  gage,  with  a 
gage  element  length  of  1.57  mm  and  width  of  1.57  mm.  The  gages  were 
mounted  as  shown  in  Figure  C-l,  3  to  4  mm  above  the  crack  tip.  In  one 
case  two  gages  were  used,  one  mounted  1.6  mm  and  the  other  5.6  mm  from 
the  fatigue  crack  tip. 

The  results  of  the  comparisons,  summarized  in  Table  C-l,  show  that 
the  difference  between  estimates  from  one  data  reduction  method  and  the 
other  never  exceeds  16%,  If  we  consider  only  the  data  obtained  from 
fatigue-precracked  experiments,  the  differences  are  even  smaller  and  do 
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not  exceed  5%.  In  particular,  very  good  agreement  was  obtained  for  the 
test  where  two  independent  strain  histories  were  recorded  (specimen  KIcj- 
16-2).  For  the  gage  closest  to  the  crack  tip,  the  maximum  stress 
intensity  value  estimated  with  the  singularity  solution  method  agreed 
closely  with  the  value  estimated  from  the  static  calibration  method  (1.9% 
difference).  A  comparison  between  the  stress  intensity  values  obtained 
from  gage  1  and  2  using  the  static  calibration  method  gives  an  estimate 
of  the  error  associated  with  this  method.  This  error  is  about  6%  and 
probably  corresponds  to  the  random  error  associated  with  the  static 
calibration  procedure  itself.  From  these  data,  we  conclude  that,  for 
sharp  cracks,  the  singularity  solution  method  gives  results  that  are  as 
reliable  as  the  static  calibration  method. 

Larger  differences  are  observed  between  the  two  methods  when 
specimens  with  blunt  notches  are  tested.  No  particular  trend  in  the 
differences  could  be  established  with  specimen  size  and  geometry. 

Because  the  notch  root  radius  used  was  relatively  large  compared  with  the 
distance  from  the  gage  center  to  the  notch  tip,  we  believe  that  the 
observed  differences  are  due  to  the  influence  of  the  notch  on  the  near 
tip  strain  field.  As  a  result  of  the  finite  root  radius  of  the  notch, 
the  strains  at  the  gage  location  are  different  from  the  strains  given  by 
the  elastic  solution  for  the  sharp  crack,  and  an  error  is  introduced  when 
this  solution  Is  used  to  reduce  the  experimental  data.  The  variability 
in  the  magnitude  of  this  error  may  be  due  to  the  differences  in  geometry 
and  surface  finish  from  notch  to  notch,  because  these  parameters  were  not 
closely  controlled.  Nevertheless,  the  data  of  Table  C-I  demonstrate  that 
the  singularity  solution  method  provides  acceptable  estimates  of  the 
stress  intensity  history,  even  for  blunt  notches. 


Estimate  of  the  Gage  Integration  Error 


To  assess  in  more  detail  the  error  associated  with  the  singularity 
solution  method,  we  estimated  the  error  due  to  the  strain  integration  by 
the  strain  gage  over  a  surface,  using  an  approach  proposed  by  B'dhme  and 
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Kalthoff.  They  estimated  the  strain  measured  by  the  strain  gage  by 
integrating  the  singularity  solution  over  the  area  of  the  strain  gage, 
that  is 
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where  i  is  the  length  of  the  active  strain  gage  element  and  b  is  its 
width. 

Dividing  equation  C-9  by  equation  C-5  gives  an  estimate  of  the  ratio 
of  measured  strain  to  true  strain. 
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Equation  C-10  has  been  evaluated  numerically  for  the  values  of  i  - 
1.57  mm  and  b  ■  1.57  mm,  corresponding  to  the  dimensions  of  the  gages 
used  in  our  experiments.  The  results  are  plotted  in  Figure  C-2  as  a 
function  of  gage  distance  from  the  crack  tip,  r.  It  can  be  seen  that  the 
strain  integration  causes  an  underestimate  of  the  strain  and  that  the 
error  drops  off  rapidly  with  r.  For  r  equal  to  3  mm  or  greater,  the 
error  is  less  than  2%.  Therefore,  even  with  relatively  large  strain 
gages,  the  gage  integration  error  is  small. 
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Effect  of  First  Higher  Order  Term  in  the  Elastic  Solution 

Bohme  and  Kalthoff^  estimated  the  error  introduced  in  the 
singularity  solution  method  by  neglecting  the  influence  of  the  first  two 
higher  order  terms.  They  found  that,  for  a  three-point-bend  specimen 
with  aspect  ratio  S/W  »  4,  and  a  relative  crack  length  a/W  between  0,2 
and  0.7,  the  error  was  always  smaller  than  10%  for  relative  gage 
positions  r/W  smaller  than  0.005.  Their  estimate,  however,  is  based  on 
the  values  of  the  higher  order  coefficients  for  the  static  case. 
Therefore,  there  is  some  question  whether  it  can  be  extended  to  the 
dynamic  case,  where  wave  reflections  may  play  a  significant  role. 


To  investigate  this  question,  we  performed  pure  mode  I  experiments 
using  the  1PBT  and  specimens  Instrumented  with  two  strain  gages, 
positioned  at  9  =  45  and  90  degrees.  Using  the  procedure  outlined  above, 
we  calculated  the  histories  of  the  mode  I  stress  intensity  and  of  the 
coefficient  of  the  first  higher  order  term  A^  (the  r-independent  term)  in 
the  series  expansion  solution  from  the  measured  strain  histories.  The 
stress  intensity  history  resulting  from  using  only  the  singular  term  in 
the  solution  was  also  calculated  based  on  the  strain  measurement  of  a 
single  gage,  the  gage  positioned  at  90°  from  the  crack  plane.  By 
comparing  the  stress  intensity  histories  obtained  from  single  gage  and 
two  gage  strain  measurements,  we  obtain  an  estimate  of  the  error 
introduced  by  neglecting  higher  order  terms  when  evaluating  strain  gage 
data. 


Figure  C-3  shows  the  results  of  an  experiment  in  which  both  K-j-(t) 
and  A-^(t)  were  evaluated  for  a  specimen  with  a  length  of  267  mm,  a  width 
of  63  mm,  and  a  crack  length  to  width  ratio  of  0.6.  The  A-]  history  is 
negative,  indicating  a  compressive  stress  acting  parallel  to  the  crack 
plane,  which  is  consistent  with  the  hammer  load  history.  Similar  results 
were  obtained  for  a  specimen  with  a  crack  length  to  width  ratio  of  0.3. 
Comparison  of  the  stress  intensity  histories  obtained  with  and  without 
taking  into  account  the  term  shows  that  the  difference  between  the  two 
curves  is  small:  the  maximum  amplitudes  differ  by  no  more  than  2%  for  the 
specimen  with  a/W  =  0.6  and  no  more  than  10%  for  the  specimen  with  a/W  = 


0.3.  We  think  that  the  10%  difference  is  due  to  the  larger  value  of  the 
notch  root  radius  (approximately  0.75  mm  compared  with  0.125  mm  for  the 
other  specimen).  Such  an  influence  of  the  root  radius  was  reported 
earlier  in  this  section  in  connection  with  the  data  in  Table  C-l. 

We  also  evaluated  the  effect  of  and  Kj  and  Kjj  in  mixed-mode 
experiments,  using  a  double-crack  1PBT  technique  described  in  Appendix  F 
and  two  or  three  gages  to  measure  strains  near  the  crack  tip.  For  the 
mode  I  experiments,  we  found  that  neglecting  A^  for  the  strain  data 
reduction  has  little  influence  on  the  Kj  history  and  introduces 
differences  of  about  2%.  The  effect  on  Kj^  appears  more  significant.  We 
calculated  differences  up  to  35%  in  the  maximum  amplitudes  as  well  as 
differences  in  the  overall  shape  of  the  KII  histories.  This  observation 
is  consistent  with  photoelastic  results  reported  by  Kobayashi  et  al.^ 
Therefore,  in  the  mixed-mode  experiments  reported  in  Appendix  F,  we 
always  used  three  strain  gages  to  determine,  KI.  KI1*  and  Ai  . 

Our  assessment  of  the  influence  of  higher  order  terms  on  the  strain 
gage  technique  for  measuring  stress  intensity  histories  is  only 
approximate  because  we  have  performed  it  for  a  limited  number  of  specimen 
sizes  and  geometries.  Furthermore,  the  accuracy  of  the.  multiple  gage 
approach  itself  has  not  been  assessed  against  other  techniques  or  against 
numerical  simulations.  Nevertheless,  we  conclude  that,  in  the  dynamic 
case  and  for  the  size  and  position  of  the  strain  gages  that  we  used, 
neglecting  the  first  higher  order  term  in  calculating  the  mode  I  stress 
intensity  history  introduces  an  error  of  less  than  5%.  This  result  is 
consistent  with  the  conclusions  based  on  the  comparison  of  the  static 
calibration  method  and  the  singularity  solution  method. 


A.  S.  Kobayashi  and  M.  Ramulu,  ’’Dynamic  Mixed-Mode  Fracture,"  in  Mixed- 
Mode  Crack  Propagation,  G.  C.  Sih  and  P.  S.  Theocaris,  Eds.,  Si jthoff  and 
Noordhoff,  Rockville,  Maryland,  1981),  pp.  163-172. 
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Our  conclusion  can  be  rationalized  by  considering  the  elastic 
solution  for  9  =  900.  Taking  only  the  singular  term  and  the  first  higher 
order  term,  the  strain  in  the  y-direction  at  a  distance  r  =*  2.5  mm  is 

e  =  -i-  (7.61  K  -  0.28  KTT  -  0.3  A,)  (C-ll) 

y  E  I  II  I 

where  a  value  of  0.3  has  been  assummed  for  Poisson’s  ratio.  Equation 
(C— 11)  indicates  that,  if  the  strain  gage  is  positioned  immediately  above 
the  crack  tip,  the  y-strain  will  be  controlled  predominantly  by  the  mode 
I  stress  intensity  factor,  unless  the  mode  II  stress  intensity  factor  and 
the  coefficient  of  the  first  higher  order  term  are  large.  Thus,  the 
position  of  the  measuring  strain  gage  in  Figure  C-l  represents  an  optimum 
position. 

To  verify  the  basic  assumption  underlying  the  static  calibration 
method  namely,  that  the  higher  order  coefficients  vary  proportionally  to 
Kj(t),  we  consider  the  ratio  A^(t)/Kj(t).  This  ratio,  obtained  from  the 
data  in  Figure  C-3,  is  plotted  in  Figure  C-4.  Clearly,  the  ratio  A^/K-j 
is  not  constant  during  the  experiment,  so  the  assumption  underlying  the 
static  calibration  method  is  violated.  However,  from  equation  (C-ll),  we 
estimate  that  the  variation  around  a  mean  value  of  the  ratio  of  -1.5  will 
cause  errors  of  6%  in  the  value  of  Kj.  What  remains  to  be  shown  is 
whether  the  mean  value  of  A^/k^-  is  consistent  with  the  value  of  the  ratio 
In  the  static  case. 


We  have  not  directly  assessed  the  gage  position  error  (coordinates 
and  angular  position  with  respect  to  the  y-axls).  Nevertheless,  the  good 
agreement  between  the  results  of  the  scatic  calibration  method  (which 
does  not  introduce  a  gage  position  error)  and  the  results  of  the 
singularity  solution  method  indicates  that  this  error  is  only  a  few 
percent.  When  the  fatigue  crack  front  of  the  fracture  specimen  displays 
a  significant  curvature  or  is  not  absolutely  perpendicular  to  the 
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specimen  faces,  experience  has  shown  that  referring  the  gage  position  to 
the  actual  crack  tip  position  (on  the  specimen  side  where  the  gage  is 
mounted),  rather  than  to  an  average  crack  tip  position,  yields  more 
consistent  data  with  less  scatter. 
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Summary  of  Results 

The  evaluations  of  the  two  strain  gage  methods  of  measuring  dynamic 
stress  intensity  histories  discussed  above  show  that  neglecting  the 
effect  of  the  higher  order  terms  Introduces  an  error  in  both  the 
singularity  solution  method  and  the  static  calibration  method.  However, 
this  error  is  small  (~  5%)  for  the  mode  I  stress  intensity  factor.  For 
the  singularity  solution  method,  the  combined  effect  of  the  integration 
error  and  the  gage  position  error  is  also  only  a  few  percent,  as  assessed 
by  comparison  with  the  static  calibration  method  results.  Therefore,  the 
strain  gage  method  provides  an  attractive  alternative  to  the  more 
involved  stress  optic  methods  for  determining  dynamic  initiation 
toughness  values.  Furthermore,  at  a  comparable  degree  of  accuracy,  the 
singularity  solution  method  i9  simpler  to  apply  than  the  static 
calibration  method  and  is  therefore  more  advantageous  to  use. 
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FACTORS  AFFECTING  THE  STRAIN  HISTORY  AFTER  INITIATION  AND  DETECTION  OF 
THE  POINT  OF  CRACK  INITIATION  ON  THE  STRAIN  RECORDS 

An  important  step  in  estimating  the  dynamic  fracture  toughness  with 
the  1PBT  is  determining  the  instant  of  crack  initiation  from  the  near 
crack  tip  strain  history  record.  The  basis  for  determining  the  instant 
of  crack  initiation  is  that  crack  extension  will  cause  a  clear  change  in 
the  slope  of  the  strain  history  record.  As  the  crack  moves  away  from  the 
strain  gage  position,  the  strain  will  decrease  rapidly  and  this  drop  can 
be  used  to  pinpoint  the  instant  of  crack  initiation.  Once  the  crack 
starts  to  extend,  however,  the  following  factors  can  influence  the  strain 
history: 

•  The  gage  position  with  respect  to  the  crack  tip. 


C-18 


^T  i*  ^y^7^TTT7  'TJ  57WT  •'.  ^CTTVl  T.'O'.Y.WMV  r  * . 


.',  V  v  ^  «  r.»  .»  .» 


k  r*  ’».*  •'.*  v 


•  The  crack  speed  and  the  resistance  of  the  material  to  crack 
extension. 

•  The  applied  stress  intensity  rate. 

•  Plasticity  effects,  residual  strains,  and  the  formation  of  shear 
lips . 

’  Nonsingular  terms  In  the  series  expansion  solution. 

In  general,  the  effect  of  all  these  factors  will  be  combined  in  a  complex 
way  to  yield  the  experimental  strain  history  record. 

We  have  investigated  how  these  factors  influence  the  strain  history 
record,  using  an  approximate  model  of  the  1PBT  to  simulate  the  strain 
history  after  the  crack  has  started  to  extend.  A  detailed  description  of 
the  model  is  given  in  Appendix  D.  The  model  allows  us  to  load  the  crack 
at  different  rates  and  to  vary  rhe  behavior  of  the  dynamic  propagation 
toughness.  The  near  crack  tip  strain  history  is  calculated  using  the 
singular  term  of  the  elastic  series  expansion  for  a  stationary  crack. 

Figure  C-5  shows  the  results  of  simulations  of  the  strain  at  the  tip 

of  a  crack  in  a  1PB  specimen  63.5  mm  wide  and  241  mm  long  with  a  crack 

length  of  19  mm.  The  simulated  impact  velocity  was  50  m  s~^. 

Simulations  were  performed  to  calculate  the  strain  at  two  positions  near 

the  crack  tip  (x  =  -1  mm,  y  =  5.9  mm,  and  x  =  1  mm,  y  **  5.9  mm)  assuming 

three  different  behaviors  for  the  propagation  toughness  In  the 

first  case,  was  assumed  constant  and  equal  to  the  initiation 

1/2 

toughness  A  value  of  58  MPa  ra  '  was  chosen  for  KI(j,  representing 

the  initiation  toughness  of  4340  steel  heat  treated  to  HRC50t  In  the 
second  case,  the  propagation  toughness  varied  with  velocity  as 


Kid  =  58  x  106  +  395  l  2 


(C-12) 


In  the  third  case,  the  propagation  toughness  was  assumed  to  rise  with 
crack  extension  over  the  first  millimeter  of  crack  growth  as 


K  =  58  x  10  (1  +  Aa  250),  Aa  <  1  mm 


(C-13) 


=  72.5  x  10  ,  Aa  >  1  mm 
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Beyond  1  mm,  was  assumed  constant  equal  to  72.5  MPa  The 

behaviors  represented  by  equations  C-12  and  C-13  were  chosen  arbitrarily 
and  are  not  intended  to  be  those  of  4340  steel. 

For  the  first  case  of  a  constant  propagation  toughness,  the 
simulation  shows  that  the  strain  can  still  increase  slightly  after  crack 
initiation.  If  the  strain  is  measured  1  mm  ahead  of  the  crack  tip,  the 
overshoot  will  be  somewhat  larger  and  last  longer  than  when  the  strain  is 
measured  1  mm  behind  the  crack  tip.  This  behavior  becomes  more 
pronounced  for  the  propagation  toughness  assumptions  in  equations  C-12 
and  C-13.  There  is  a  significant  overshoot  of  the  strain  measured  1  mm 
ahead  of  the  crack  tip,  whereas  very  little  or  no  overshoot  is  observed 
for  the  strain  measured  1  tun  behind  the  crack  tip.  More  important,  there 
is  a  considerable  delay  in  the  drop  in  strains  when  we  assume  an 
increasing  propagation  toughness  with  crack  extension  or  crack  velocity, 
compared  with  the  case  of  a  constant 

Figure  C-6  shows  a  set  of  simulations  similar  to  those  plotteu  in 
Figure  C-5,  but  for  a  specimen  63  mm  wide  and  356  mm  long  with  a  crack 
length  of  19  mm.  The  simulated  impact  velocity  was  3.5  m  s-1,  resulting 
in  a  much  smaller  stress  intensity  rate  than  in  Figure  C-5»  Therefore,  a 
comparison  of  Figure  C-5  and  C-6  illustrates  the  effect  of  changes  in  the 
stress  Intensity  rate.  In  the  simulations  of  Figure  C-6,  the  applied 
stress  intensity  history  was  such  that  crack  arrest  was  predicted  in  all 
the  simulations  after  an  amount  of  crack  extension  varying  depending  on 
the  propagation  toughness  assumption.  For  the  case  of  the  constant 
propagation  toughness  and  the  case  of  the  toughness  varying  according  to 
equation  C-13,  the  strain  drop  appears  sharper  than  for  the  simulation  at 
higher  rate.  However,  this  difference  is  essentially  due  to  the  change 
in  the  time  scale  of  the  plots  associated  with  the  different  stress 
intensity  rates.  For  the  case  of  the  propagation  toughness  rising  with 
the  amount  of  crack  extension,  very  little  crack  extension  occurred,  and 
the  simulated  strain  history  lies  very  close  to  the  strain  history  for 
the  specimen  with  a  stationary  crack.  In  such  a  case,  the  instant  of 
crack  initiation  would  be  difficult  to  establish  unambiguously. 


C-21 


o 


o 

O 

o 

o 

o 

O 

o 

o 

Tt 

00 

C*J 

CD 

(N 

T— 

^ — 

v: 

o 

< 

a: 

_i 

< 

Z 


> 

cc 

O 

H- 

CO 


< 
QC 
►  - 
CO 


CO 

z 

o 


3 

CO 

00 

< 

CO 

CO 

LL) 

z 

X 

o 

3 

o 


e> 

< 

o 

o 


°  I 


CD 

1 

<J> 

LU 

cc 

3 

eo 

LL 


.vy\„NV>a 

k  - ■  vj 
‘.'.v-AV'J 


.• 

■V-.  r 

.«.v;.cv 


■V>A*sj 


•  -  ^  -j 


•9' 


y.Vu'v 

^-■4 


V.V.V.v 

'  *.V.v  v 

■-  v  '.\v 

V  V.  >/* 

'•>  \y,\\’ 


r«* "c  ^  i  *■ 
.  A  V.< 
>  -  V  V 

•  Vv.v. 


>V'V  .> 


*  *  * 
-"a 


'  ■  *  V  »  *  I 


*9^*1 


-  v; 


V  - 

■■.  A  -r.  ■ 

•  .V  M  "  M  - 


VVV"’ 


-  V  >  , 


m.: 


(uiejjsojDjuj)  Nivyis 


L-WL-K  v.>  V.~*»  \T* i  kH»  v.’X  L">  L'-»  L“V  L"V  L-*’  V*  ^  w  u  *  ’.  *  u  ■  i>  -  j. 


Features  similar  to  those  predicted  by  the  strain  simulations  In 
Figures  C-5  and  C-6  have  been  observed  In  the  experimental  strain  history 
records.  Figures  C-7  and  C-8  show  examples  of  strain  records  for 
experiments  with  the  same  specimen  geometries  as  the  specimens  simulated 
in  Figures  C-5  and  C-6.  In  Figures  C-7  and  C-8,  each  of  the  curves  was 
obtained  from  identical  but  separate  tests  with  different  gage  positions 
with  respect  to  the  crack  tip.  The  experimental  strain  records  have  also 
been  normalized  by  the  strain  at  initiation.  In  these  records,  the 
strain  history  after  initiation  is  significantly  affected  by  the  gage 
position.  Moreover,  some  of  the  experimental  records  also  display  a 
strain  overshoot  and  a  delayed  drop  in  the  strain,  compared  with  the 
strain  simulations  for  the  constant  propagation  toughness. 

Several  conclusions  can  be  drawn  from  the  strain  record 
simulations.  First,  because  of  the  potential  for  a  strain  overshoot 
after  crack  initiation,  we  must  be  careful  in  determining  the  instant  of 
crack  initiation  from  the  strain  history  records  to  avoid  overestimating 
the  dynamic  initiation  toughness.  In  particular,  the  instant  of  maximum 
strain  may  not  correspond  to  the  point  of  crack  initiation.  Rather,  the 
point  of  separation  of  the  strain  record  from  the  reference  curve  for  the 
specimen  without  crack  extension  is  a  better  indication  of  the  instant  of 
crack  initiation.  If  the  reference  curve  is  not  available,  it  can  be 
approximated  by  a  tangent  to  the  strain  record  in  the  region  immediately 
preceding  the  point  of  inflection. 

The  second  conclusion  drawn  from  the  strain  record  simulations  is 
that  the  toughness  characteristics  of  the  material  tested  strongly  affect 
the  strain  history  after  crack  initiation.  However,  the  simulations  also 
Indicate  that,  given  a  particular  gage  location,  the  amount  of  crack 
extension  required  to  drop  the  strain  to  zero  depends  only  mildly  on  the 
propagation  toughness  characteristics  and  hence  can  be  assumed  constant 
in  first  approximation.  These  observations  indicate  that,  as  suggested 
by  Bonme  and  Kalthoff,^  information  about  the  dynamic  propagation 
toughness  and  the  initial  crack  velocity  can  be  obtained  from  the 
analysis  of  the  strain  records. 


».  sy. 


Another  factor  that  may  affect  the  strain  history  after  crack 
initiation  I3  plasticity  in  the  crack  tip  region.  In  particular,  the 
formation  of  small  shear  lips  at  the  free  surfaces  of  the  specimen  may 
affect  the  unloading  of  the  crack  tip  strain  gage.  Indeed,  when  shear 
lips  are  present,  they  are  invariably  accompanied  by  a  certain  degree  of 
crack  tunneling  so  that  the  tip  of  the  crack  in  the  shear  lip  lags  behind 
the  tip  of  the  crack  at  the  center  of  the  specimen.  Therefore,  shear 
lips  represent  unbroken  material  ligaments  that  may  continue  loading  the 
strain  gage,  although  most  of  the  crack  front  has  already  passed  the  gage 
position. 

To  investigate  this  possibility,  we  performed  1PB  experiments  with 
small  4340  steel  specimens.  The  material  was  heat  treated  to  a  hardness 
HRC  50.  The  specimens  were  32  mm  wide  and  133  mm  long,  13  ram  thick,  and 
had  a  crack  length  of  13  mm.  Some  of  the  specimens  were  provided  with 
small  side-grooves  (0.5  mm  wide,  1  mm  deep)  to  suppress  the  formation  of 
shear  lips.  In  the  rest  of  the  specimens,  which  had  smooth  side 
surfaces,  1-mm-wide  shear  lips  were  formed  upon  fracturing.  All  the 
specimens  were  instrumented  with  a  strain  gage  positioned  3  to  3.5  mm 
from  the  crack  tip  and  as  Suown  in  Figure  C-l.  The  static  calibration 
method  was  used  to  obtain  the  stress  intensity  history  from  the  strain 
records  • 

Figure  C-9  compares  the  strain  history  records  for  side-grooved  and 
ungrooved  specimens.  For  the  records  shown,  the  gage  position  was  about 
0.3  mm  behind  the  crack  tip.  The  strain  drops  sooner  in  the  side-grooved 
than  in  the  ungrooved  specimens,  suggesting  that  the  formation  of  even 
small  shear  lips  may  have  a  significant  influence  on  the  strain  history 
after  crack  initiation.  This  influence  may  be  direct  or  indirect.  In 
the  first  case,  the  shear  lips  acts  as  a  small  tensile  ligament  behind 
the  crack  front  that  continues  loading  the  strain  gage.  According  to 
this  explanation,  the  shear  lip  would  affect  the  strain  history  but  not 
the  effective  propagation  toughness  and  the  crack  velocity.  In  the 
second  case,  the  formation  of  the  shear  lips  raises  the  apparent 
toughness  of  the  material  and  slows  the  crack  extension.  Therefore,  the 
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FIGURE  C-9  NORMALIZED  EXPERIMENTAL  STRAIN  HISTORY  NEAR  THE 
INITIAL  CRACK  TIP  FOR  SIDE-GROOVED  AND  UNGROOVED 
SPECIMENS 


strain  record  for  the  ungrooved  specimen  in  Figure  C-9  would  drop  off 
later  than  the  strain  record  for  the  side-grooved  specimen  because  there 
is  a  crack  velocity  difference  in  the  two  testa.  Such  a  crack  velocity 
difference  has  indeed  been  observed  in  experiments  in  which  the  average 
velocity  during  the  first  9  mm  of  extension  was  measured  with  strain 
gages.  Therefore,  we  conclude  that  the  shear  lips  have  an  influence  on 
the  strain  records,  but  that  this  influence  reflects  an  actual  rise  in 
the  effective  propagation  toughness.  Plasticity  at  the  crack  tip  may 
also  affect  the  strain  history  by  Inducing  residual  strains  which  then 
prevent  a  complete  unloading  of  the  strain  gage.  Residual  strains  will 
occur  particularly  if  the  gage  is  mounted  too  close  to  the  crack  tip. 

Higher  order  terms  in  the  elasto-dynamic  series  expansion  solution 
and  transient  conditions  set  up  by  stress  waves  could  also  affect  the 
strain  history  after  crack  initiation.  These  effects  can  be  studied  only 
by  detailed  photoelastic  experiments  or  by  numerical  simulations  of  the 
complete  specimen. 

CONCLUSIONS 

Our  analysis  of  the  two  strain  gage  methods  of  determining  stress 
intensity  factor  histories  in  dynamic  fracture  tests  indicates  that  both 
methods  should  give  stress  Intensity  estimates  within  10%  or  better  of 
the  true  value.  The  static  calibration  method  and  the  singularity 
solution  method  give  very  consistent  results  that  are  within  5%  of  each 
other,  provided  the  strain  gage  is  positioned  carefully.  Therefore,  the 
singularity  solution  method,  which  is  simpler  to  implement,  can  be  used 
with  confidence.  An  angle  of  90  degrees  from  the  crack  plane  and  a 
distance  of  3  to  4  ram  from  the  crack  tip  represent  an  optimum  position 
for  the  strain  gage,  for  specimens  with  length-over-width-ratios  greater 
than  2.6  and  widths  down  to  25  mm. 

Experimental  estimates  of  the  first  higher  order  term  in  the  series 
expansion  using  multiple  strain  gages  indicate  that,  with  the  optimum 
gage  position,  the  influence  of  this  term  on  the  mode  I  stress  intensity 
is  less  than  5%  and  can  be  neglected.  However,  for  the  mode  II  stress 
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intensity  history,  errors  of  about  35%  have  been  observed  for  Kj^  by 
neglecting  the  effect  of  the  first  higher  order  term. 

Simulations  of  strain  histories  near  the  crack  tip  after  initiation 
show  that,  depending  on  the  precise  gage  positon  and  the  material 
toughness,  the.  strain  can  continue  to  rise  significantly  after 
initiation.  Therefore,  the  instant  of  crack  initiation  does  not  always 
correspond  to  the  maximum  strain  and  must  be  determined  by  comparing  the 
stress  intensity  history  for  the  fracture  curve  with  the  reference  stress 
intensity  history  (without  a  crack  extension).  Finally,  the  strain 
history  simulations  suggest  that  analysis  of  the  strain  history  after 
initiation  can  provide  estimates  of  the  crack  velocity  and  information 
about  the  material  propagation  toughness. 
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ve  rrom  experimental  recoras  mat  tne  stress  intensity  nistory  ior 
PBT  has  a  smooth  approximately  sinusoidal  shape,  controlled  mainly 
e  first  mode  of  oscillation  of  the  specimen.  Because  later 


intens 


specimens  with  stationary  cracks  of  length  equal  to  the  current  length  of 
the  extending  crack.  In  estimating  the  stress  intensity,  we  do  not 
consider  any  interaction  of  stress  waves  reflected  from  the  specimen 


boundaries  with  the  crack  tip.  We  do,  however,  consider  material  inertia 
effects  at  the  crack  tip;  they  are  modeled  using  a  linear  approximation 
to  the  expression  derived  by  Freund ^  for  a  crack  extending  in  an  infinite 
plate. 

The  proposed  model  is  implemented  on  a  desk  top  computer  and 
calculates  stress  Intensity  and  crack  extension  histories 
incrementally.  A  listing  of  the  computer  program  is  given  at  the  end  of 
this  Appendix.  Given  the  initial  impact  velocity,  two  approaches  can  be 
taken:  either  the  crack  velocity  history  is  known  and  the  resulting 
stress  intensity  history  is  calculated  or  the  material  toughness  behavior 
is  known  and  the  resulting  crack  velocity  history  is  calculated.  The 
following  sections  give  a  step-by-step  account  of  the  model  development. 


PREDICTION  OF  THE  STRESS  INTENSITY  HISTORY  FOR  SPECIMENS  WITH  STATIONARY 
CRACK  LENGTH  GREATER  THAN  TIE  INITIAL  CRACK  LENGTH 


We  start  by  assuming  that,  in  a  manner  similar  to  the  static  case 

§ 

for  the  three-point-bend  test,  the  stress  intensity  history  K^for  the 
1PBT  with  a  stationary  crack  can  be  expressed  in  terms  of  the  crack  mouth 
opening  as 


,s  _  imp 

I=V? 

imp 


6  (a/wi  f(£_)  Ijl.  n^y). 

max  U/W'  UT  '  , -  V(a/w) 

o  /  a 


(D-l) 


where  F(a/W)  and  V(a/W)  are  the  static  compliance  functions  relating  the 
stress  intensity  and  the  crack  mouth  opening  displacement,  respectively, 
to  the  applied  load.  The  term  f(t/TQ)  is  a  shape  function  representing 
the  shape  of  the  stress  intensity  history  for  the  specimen  with  the 
initial  crack  length.  It  is  chosen  to  match  thp  experimentally  recorded 
stress  intensity  history  as  closely  as  possible,  within  practical  limits, 
with  a  set  of  simple  functions.  The  shape  function  has  a  maximum,  and 
the  time  at  which  this  maximum  occurs  is  denoted  t  .  As  noted  above, 
the  experimental  stress  intensity  histories  for  specimens  of  varying 
crack  length  have  a  common  initial  portion  (Figure  D-l).  The  time  at 
which  the  stress  intensity  history  for  the  specimen  with  the  initial 
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crack  length  starts  diverging  from  the  other  stress  intensity  histories 
is  denoted  t*. 

The  term  ^max(a/w)  in  equation  (D-l)  is  an  amplitude  parameter  that 
increases  with  crack  length.  The  value  of  the  amplitude  parameter  for 
the  specimen  with  initial  crack  length  is  adjusted  to  match  the 
experimental  K™ax  and  is  an  input  to  the  model.  v°mp  is  the  impact 
velocity  for  which  the  reference  stress  Intensity  history  for  the 
specimen  with  initial  crack  length  has  been  obtained.  The  product 

^max  f(t/TQ)  represents  the  crack  mouth  opening  history 


6(imE/W)or  more  precisely, 


6<a/w.t)  5  ji  -i!E-  ^  (./«)  f(t/T0) 
imp 

A  2 

where  E'  t3  equal  to  E  in  plane  stress  and  to  E/(l  -  v  )  in  plane 
strain. 


(D-2) 


From  experimental  observation,  we  have  found  that  a  shape  function 


of  the  form 


f(s-)  -  4-  [l  +  sJo  |2-  (t  -  T^)  o  <  t  <5/8  T 


f(|-)  -  co,|-(t  -ItJ  ?To<t<Tc 
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(D-3) 


gives  a  suitable  representation  of  the  stress  intensity  history  for  a 
wide  variety  of  specimen  geometries.  In  equation  (D-3),  tmax  has  been 
taken  as  (5/8)To.  Also  based  on  experimental  observations,  we  assume 
that  t*  is  given  in  terras  of  TQ  as  (7/16)TQ.  In  deriving  the  stress 
intensity  histories  for  specimens  with  longer  crack  lengths,  we  will 
constrain  all  the  curves  to  pass  through  this  point.  Further,  we  assume 
that  the  same  shape  function  applies  for  every  crack  length,  provided  the 
duration  TQ  is  made  a  function  of  crack  length. 

Equation  (D-l),  together  with  the  experimental  observation  that  the 
initial  portion  of  the  stress  intensity  history  is  independent  of  crack 
length,  implies  that: 
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The  crack  mouth  opening  history  must  have  the  same  shape 
function  as  the  stress  intensity  history. 


At  a  given  time,  the  crack  mouth  opening  must  increase  with 
increasing  crack  length,  as  illustrated  in  Figure  D-2. 


These  points  will  be  discussed  further  when  comparing  the  model 
predictions  with  experimental  results. 


From  the  foregoing  discussion,  it  follows  that  predicting  the  stress 
intensity  history  for  specimens  with  stationary  crack  lengths  longer  than 
the  initial  crack  length  simply  amounts  to  predicting  the  maximum 
^mp,  a/W),  or  equivalently  the  amplitude  parameter 


amplitude  K^aX(V 


6  (V.  ,  a/W),  and  the  duration  T  (a/W). 

max  imp  ° 


The  condition  that  at  time  t*  all  stress  intensity  histories  must 
have  the  same  value  provides  the  first  equation  to  calculate  these  two 
unknowns : 


KS(t  ,  a1/W)  -  Ka  (t  ,  a/W) 


<0-4 ) 


A  second  condition  is  derived  from  momentum  considerations.  We  consider 

the  specimen  as  made  up  of  two  rigid  halves  connected  elastically  at  the 

crack  plane,  as  shown  in  Figure  D-3.  Upon  impact,  each  half  of  the 

specimen  is  imparted  a  translational  velocity,  Vg(t),  and  a  rotational 

velocity  w  (t).  Under  these  conditions,  conservation  of  linear  momentum 
s 

requires  that 

t 

mV  (t)  =  J  F(t" ,a/W)dt"  (D-5) 

0 

where  F(t',a/W)  is  the  force  applied  to  the  specimen  by  the  impacting 
hammer  and  ra  is  the  mass  of  the  specimen.  Further,  conservation  of 
angular  momentum  requires  that 

2  t  t 

yy-w8(t)  =  /  F(t",a/W)  ~y~  dt'  -  /  M(t',a/W)dt'  (D-6) 

0  0 
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where  M(t',a/W)  is  the  resisting  moment  on  the  uncracked  ligament  and  L 
is  the  specimen  length. 

At  time  tmax,  the  whole  specimen  has  the  same  uniform  velocity 

V„(t_„)  and  w  (t  )  ■  0.  Thus,  equations  (D-5)  and  (D-6)  become 
s'  max'  s  max  ' 

t 

J  maX  F(t',a/W)dt"  -  mV  (t  )  (D-7) 

b  max 

0 

and 

1  7  t 

j  raax  F(t',a/W)-4p- dt'  -  ^  ®ax  M(t',a/W)df'  (D-8) 

0 

respectively.  If  we  assume  that  Vg  (tmax)  is  a  function  of  the  initial 
impact  velocity,  only,  and  is  independent  of  the  crack  length  of  the 
specimen,  then  equation  D-8  can  be  rewritten 
t 

max 

/  M(t",a/W)  dt'  •  constant  (D-9) 

0 

The  assumption  that  Vg(tmax)  is  independent  of  the  specimen  crack 
length  implies,  through  equation  (D-7),  that  the  total  momentum  applied 
to  the  specimen  up  to  the  point  of  maximum  stress  intensity  is  constant 
and  also  independent  of  crack  length.  This  result  will  be  compared  with 
experimental  observations  when  discussing  the  model. 

Finally,  assuming  that  the  moment  on  the  crack  plane  ligament  and 
the  stress  intensity  factor  can  be  directly  related  as  in  the  quasi- 
static  case,  that  is, 

K®  (t,a/W)  ~  M(t,a/W)F(a/W)  /“xa  (D-10) 

equation  (D-9)  reduces  to 
t  Ks 

J  raax  i  (f*  ,a/W)  d). 

0  F(a/W)/Ta 


constant 


(D-ll ) 


In  equations  (D-10)  and  (D-ll),  whe  have  used  the  9tatlc  compliance 
function  F(a/W),  which  Is  consistent  with  the  assumption  underlying 
equations  (D-l). 


The  constant  on  the  right-hand  side  of  equation  (D-ll)  can  be 

evaluated  by  Integrating  the  stress  Intensity  history  of  the  specimen 

with  the  Initial  crack  length.  Thus,  introducing  t  *  •a-  T 

max  8  o 


5/8T  (a/W)  K®(t',a/W) 

J  °  — - —  dt' 

0  F(a/W)  /la 
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(D-12) 

Equations  (D-4)  and  (D-12),  combined  with  equations  (D-l)  and  (D-3), 
provide  the  two  equations  with  which  to  calculate  l-)max(a/w)  and  T  (a/W). 

Having  expressed  6  (a/W)  and  T  (a/W),  we  have  completed  modeling 

max  o 

of  the  stress  intensity  history  for  specimens  with  stationary  crack 
lengths  larger  than  the  Initial  crack  length. 


EVALUATION  OF  THE  STRESS  INTENSITY  HISTORY  AT  THE  RUNNING  CRACK  TIP 

Once  the  material  dynamic  initiation  fracture  toughness  has  been 
exceeded  In  a  1PBT  experiment,  the  crack  will  start  to  propagate  at  a 
speed  a(t).  The  energy  to  drive  the  propagating  crack  is  supplied  by  the 
inertial  bending  of  the  specimen.  Therefore,  to  evaluate  the  stress 
intensity  at  the  extending  crack  tip,  K^(t,a,a),  we  extend  the  approach 
taken  in  evaluating  the  stress  intensity  history  at  the  stationary  crack 
tip. 

Again,  we  assume  that  the  stress  intensity  can  be  expressed  in  terms 
of  the  crack  mouth  opening  displacement  as 


Kj(t ,a,a) 


=  6 

max 


(a/W)  f( 


TQ(a/W) 


F(a/Wj_  A 
1  -  V(a/W)  ^ 


(D-13) 


approximate  way  for  material  inertial  effects  at  the  crack  tip*  Based  on 
Freund's^  analysis  of  a  crack  extending  in  an  infinite  plate  at  constant 
velocity,  and  following  an  approach  proposed  by  Homma  et  al.,  we  take 

A(a/CR)  -  1  -  (4-  (D-14) 


where  CR  is  the  Rayleigh  wave  speed 


FIGURE  D-4  SCHEMATIC  OF  CRACK  MOUTH  OPENING  BEHAVIOR 
ASSUMED  IN  1PBT  MODEL  FOR  SPECIMEN  WITH 
PROPAGATING  CRACK 


length  and  then  Introduced  In  equation  (D-13)  to  obtain  (t,a,a).  The 
process  Is  repeated  for  the  next  time  steps. 


In  simulations  where  the  crack  velocity  history  is  not  prescribed, 

p  •  e 

the  crack  Is  allowed  to  extend  as  long  as  K^(t,a,a)  exceeds  KjD(Aa,a). 


However,  K^(t,a,a)  cannot  increase  arbitrarily.  Depending  on  the  impact 


velocity  and  the  material  toughness,  inertial  bending  of  the  specimen  may 
stop  before  the  crack  reaches  the  specimen  boundary.  Referring  to  the 
simplified  model  of  Figure  D-3,  this  is  the  case  when  the  applied  force 
F(t,a/W)  and  the  moment  on  the  remaining  ligament  M(t,a/W)  have  brought 
the  specimen  to  an  instantaneous  uniform  velocity.  Consistent  with  the 
approach  taken  in  the  preceding  section,  we  assume  that  the  point  of 
maximal  inertial  bending  for  the  specimen  with  an  extending  crack  is 
reached  when  the  following  condition  is  met: 


a,' 


'  dt'’  +  / 
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(D-16) 


At  each  time  step  in  the  simulation,  the  left-hand  side  of  equation 
(D-16)  is  evaluated,  taking  into  account  that  now  the  crack  length  is 
also  a  function  of  time.  When  equation  (D-16)  is  satisfied,  the  crack 


mouth  opening  6  (a/W)f(t/T  (a/W))  is  made  to  decrease  using  the  cosine 
max  o 


expression  in  (D-3)  for  the  shape  function  f(t/TQ). 


If  the  condition  In  (D-16)  is  not  met  before  the  crack  has  extended 


to  an  a/W  value  of  0.7,  the  calculation  is  stopped.  Experimental 


3  a 

observations'*  and  numerical  results4  indicate  that,  once  the  crack 


extends  beyond  an  a/W  value  of  0.7,  the  effect  of  the  compressive  zone 
under  the  impact  point  will  significantly  affect  the  propagation 
behavior.  The  effect  of  the  compressive  zone  is  not  considered  in  our 
simple  approach. 
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SIMULATION  OF  STRAINS  IN  THE  VICINITY  OF  THE  STATIONARY  AND  THE 
PROPAGATING  CRACK 

One  of  the  objectives  in  developing  the  model  of  the  1PBT  was  to 
simulate  the  effect  of  crack,  extension  on  the  strain  records  obtained 
with  strain  gages  placed  near  the  initial  crack  tip.  Improved 
understanding  of  the  factors  affecting  these  records  allows  a  better 
determination  of  the  instant  of  crack  initiation.  A  careful  analysis  of 
the  drop  in  strain  after  crack  initiation  may  also  provide  information 
about  the  resistance  of  the  material  to  crack  extension,  as  suggested  by 
Bourne  and  Kalthoff.*’ 

Knowing  the  stress  intensity  at  the  tip  of  the  extending  crack,  we 
evaluate  the  strain  history  at  two  locations:  near  the  initial  crack  tip 
and  at  a  small  distance  from  the  crack  plane  on  the  path  of  the  extending 
crack.  These  positions  and  their  coordinates  are  shown  in  Figure  D-5. 


Simulation  of  Strain  History  Near  the  Initial  Crack  Tip 

The  strain  in  the  direction  y,  normal  to  the  crack  plane,  Figure  D- 
5,  is  evaluated  at  a  position  representative  of  the  strain  gage  position 
in  the  1PBT.  The  position  is  defined  by  the  polar  coordinates  rQ  and 
0^  with  the  coordinate  system  centered  at  the  crack  tip  (Figure  D-5). 

The  strain  is  then  calculated  using  the  singular  term  of  the  elastic 
solution  for  the  fields  at  the  stationary  crack  tip  to  relate  the  strain 
to  the  dynamic  stress  intensity  factor  evaluated  with  the  1PBT,  that  is. 


1  KI(t)  r  3 

e  =  — g - [l  —  v  +  (1  +  v)sin9sin  y  9j 

^  /2itr 


( D— 17 ) 


Up  to  the  time  of  crack  initiation,  t£nit,  the  polar  coordinates  r 
and  6  are  held  constant.  After  crack  initiation,  they  are  updated  to 
account  for  crack  extension,  with 


[r  cos9  -  a(t  -  t.  .„)]  +  (r  sinB  ) 

1  o  o  '  init  ‘  o  o 


(D-1S) 


D-14 


where 


plate  wave  speed  (plane  stress) 
shear  wave  speed 


Strain  histories  evaluated  at  several  positions  along  the  crack  plane 
using  different  assumptions  for  the  propagation  toughness  behavior 
can  be  compared  with  experimental  strain  gage  measurements  and 
conclusions  can  be  drawn  regarding  the  actual  material  propagation 
toughness. 

VERIFICATIONS  OF  1PBT  MODEL 

To  evaluate  the  validity  of  the  1PBT  model  and  of  its  underlying 
assumptions  we  have  performed  several  comparisons  with  experimental 
results,  which  are  presented  below.  Some  of  the  experimental  results 
also  provide  additional  insight  in  the  1PBT  specimen  dynamics. 


Comparison  of  Hammer  Load  History  for  Specimens  of  Different  Crack 
Lengths  " 

In  deriving  the  1PBT  model,  we  also  assumed  that  the  total  momentum 
transferred  up  to  the  maximum  In  the  stress  intensity  curve  is 
independent  of  the  specimen  crack  length.  Thus  the  time  integral  of  the 
hammer  load  histories  up  to  the  maximum  should  be  the  same.  To  verify 
this  point,  we  measured  the  hammer  load  for  two  sets  of  specimens  with 
different  crack  lengths.  The  first  two  specimens  were  229  mm  long  and  89 
ram  wide  and  had  a/W  ratios  of  0.3  and  0,6.  The  hammer  load  histories  for 
these  two  specimens  are  shown  in  Figure  D-7(a).  The  other  two  specimens 
were  267  mm  long  and  63  mm  wide  with  a/W  ratios  of  0,4  and  0.6.  The 
hammer  load  histories  and  the  stress  intensity  histories  for  these  two 
specimens  are  plotted  in  Figure  D-7(b).  As  for  the  displacement 
assumption,  it  Is  seen  that  the  area  under  the  hammer  load  curves  Is  not 
strictly  the  same  for  specimens  of  two  crack  lengths,  but  that  the 
differences  are  relatively  small.  Therefore,  the  assumption  made  In  the 
model  is,  in  first  approximation,  acceptable. 

Comparison  of  Experimental  and  Predicted  Str  as  Intensity  Curves 
for  Specimens  with  Different  Stationary  Crar  Lengths 

The  procedure  described  above  for  predicting  the  stress  intensity 
history  of  specimens  of  different  stationary  crack  lengths  was  used  to 
estimate  the  stress  intensity  history  of  small  specimens  (107  mm  by  25.4 
mm)  with  a/W  ratios  of  0.3,  0.4,  0.5,  and  0.6.  The  stress  Intensity 
history  for  the  specimen  with  the  a/W  =0.3  was  ^sed  as  input  for  the 
estimate. 

Figure  D-8  compares  the  experimental  results  with  the  results  of  the 
simulation.  The  predicted  maximum  amplitudes  and  the  time  at  which  these 
amplitudes  are  reached  agree  well  with  the  experimental  results.  The 
total  duration  of  the  stress  Intensity  oscillations  do  not  agree  as  well 
because  at  later  times  the  oscillation  mode  in  the  specimens  changes, 
which  is  not  considered  in  the  model.  The  good  agreement  between  model 
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COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  STRESS  INTENSITY  HISTORIES 
FOR  SPECIMENS  OF  DIFFERENT  CRACK  LENGTHS 


and  experiment  up  to  the  point  of  maximum  amplitude  gives  some  confidence 
that  the  model  simulates  the  essentials  of  the  specimen  dynamics. 

Comparison  of  Experimental  and  Predicted  Results  for  a  Fracturing  1PB 
Specimen 

A  definitive  verification  of  the  simple  1PB  model  can  only  be  made 
by  comparing  the  model  prediction  for  a  running  crack  with  experimental 
measurements  of  the  stress  intensity  history  or  with  numerical 
simulations  of  the  test  using  finite  element  or  finite  difference 
codes.  We  have  found  in  the  literature  the  results  of  a  single  1PB 
experiment  for  which  the  stress  intensity  history  and  the  crack  velocity 
were  measured  using  the  method  of  caustics  for  the  initial  part  of  crack 
extension.®  A  finite  difference  simulation  of  the  complete  crack 
extension  through  the  specimen  is  also  available.^  The  experiment  was 
performed  with  an  ARALDITE  B  specimen.  For  the  simulation,  the 
propagation  toughness  versus  crack  speed  curve  measured  in  separate 
experiments  was  used  to  govern  the  crack  extension,  and  the  predicted 
stress  Intensity  and  crack  extension  histories  were  compared  with  the 
experimetal  results.  They  were  found  in  excellent  agreement. 

To  test  the  1PB  model,  we  used  the  same  approach,  combining  the 
experimental  toughness  versus  velocity  curve  with  the  applied  stress 
intensity  estimated  with  the  model  to  predict  the  stress  intensity  and 
crack  extension  histories.  In  addition  to  the  toughnes  versus  crack 
velocity  curve,  the  IPB  model  requires  knowledge  of  the  stress  intensity 
curve  for  the  specimen  with  a  stationary  crack  of  length  equal  to  the 
initial  crack  length  In  the  propagation  experiment.  This  curve  was 
estimated  from  the  corresponding  curve  for  a  steel  specimen,  and  the 
initial  slope  was  adjusted  to  match  the  slope  of  the  ARALDITE  specimen 
measured  in  the  experiment. 

Figure  D-9  compares  the  results  of  references  4  and  8  and  of  the 
stress  intensity  and  crack  velocity  histories  predicted  by  the  1PB  model. 
During  the  initial  portion  of  crack  extension,  the  agreement  is 
remarkably  good,  given  the  rather  crude  nature  of  the  1P3  model.  At 
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later  times,  the  1PB  model  overpredicts  both  the  stress  intensity 
amplitude  and  the  crack  velocity;  no  maximum  had  been  reached  at  the  time 
the  calculation  was  terminated  when  a/W  exceeded  0.7.  This  simulation 
suggests  that  the  simple  model  of  the  1PBT  represents  the  phenomenon  of 
the  early  part  of  crack  extension  correctly.  At  later  times  or  for 
larger  crack  extension,  the  model  does  not  perform  as  well  and  should  be 
used  with  caution. 

USES  AND  LIMITS  OF  THE  1PBT  MODEL 

The  1PBT  model  was  developed  to  simulate  the  effect  of  crack 
extension  on  the  strain  history  near  the  initial  crack  tip  and  to  obtain 
engineering  estimates  of  the  propagation  toughness  without  having  to 
resort  to  expensive  code  calculations.  The  model  is  based  on  a  series  of 
simplifying  assumptions  and  on  the  phenomenological  understanding  of  the 
experiment.  It  is  based  on  quasi-static  considerations  and  does  not  take 
into  account  stress  wave/crack  tip  interactions,  which  are  known  to  be 
important  at  high  crack  velocities. 

Experimental  verification  of  the  assumptions  made  in  the  model  show 
that,  although  these  assumptions  are  violated,  they  represent  a 
reasonable  first  approximation.  Further  a  direct  comparison  with 
experimental  results  and  a  full-blown  computation  show  that  the  model 
approximates  well  the  stress  Intensity  history  for  specimens  with 
different  stationary  crack  lengths  and  gives  a  good  approximation  of  the 
initial  phase  of  crack  extension.  Therefore,  the  1PBT  model  is  well 
suited  for  the  applications  it  was  designed  for.  It  has  proved  very 
useful  in  understanding  what  affects  the  strain  history  near  the  initial 
crack  tip  and  in  evaluating  the  transition  from  the  stationary  to  the 
fast  extending  crack  (Appendices  C  and  G).  It  can  also  be  used  to 
estimate  what  specimen  geometry  and  impact  conditions  should  be  chosen  to 
achieve  certain  crack  extension  conditions  in  a  given  material. 

Beyond  these  applications  great  care  should  be  used.  The  results  in 
Figure  D-9  show  that  the  model  is  not  suited  to  making  predictions  for 
crack  extension  to  large  a/W  values  and  that  it  provides  only  an  estimate 


vV.’/';. 

sy./V's 

« v- 

vVN*- 


V ■ 

>  /  / ' 
m 


*•  -v-o 
/■  '■ 
•-V-V-V 

wav; 

‘«V\V 

Al.i’ 


-v 


•  %* 
fT  M-  .* 

'J-  -A 

m 

iV'// 


V  Wf\V 

tfSfcv' 


D-25 


of  the  stress  Intensity  and  crack,  velocity  histories.  Therefore,  for 
longer  crack  extensions  and  for  crack  arrest,  the  model  should  at  best  be 
used  to  make  qualitative  comparison  between  various  propagation 
behaviors. 


1PBT  MODEL:  LISTING  OF  COMPUTER  PROGRAMS 


10  !  PROGRAM  TO  SIMULATE  THE  PROCESS  OF  CRACK  INITIATION  AND  CRACK  PROPAGATION 

!  UNDER  DYNAMIC  LOADING  CONDITIONS  IN  THE  ONE-POINT-BEND  SPECIMEN 
20  ! 

30  !  THE  PROGRAM  CONTAINS  A  BEST  APPROXIMATION  OF  THE  EXPERIMENTALLY  MEASURED 

!  STRESS  INTENSITY  HISTORY  FOR  THE  DYNAMIC  LOADING  OF  THE  STATIONARY  CRAG 
AO  !  AND  USES  FREUND'S  EXPRESSION  FOR  THE  STRESS  INTENSITY  FACTOR  AT  THE  TIP 
!  OF  THE  RUNNING  CRACK 
SO  1 

60  !  DURING  PROPAGATION.  THE  CRACK  SPEED  IS  EVALUATE  USING  HOMMA'S  ET  AL . 

1  PROCEDURE 
70  ! 

80  !  SEVERAL  MATERIAL  TOUGHNESS  BEHAVIORS  ARE  INVESTIGATED  AND  CAN  BE  CHOSEN 

!  FROM  A  MENU 
90  ! 

100  !  DEFINITION  OF  MATERIAL  PROPERTIES:  Young  IS  YOUNG'S  MODULUS.  Mu  IS 
!  POISSON'S  RATIO,  Rho  IS  THE  MASS  DENSITY  AND  Cr  IS  THE  RAYLEIGH 
110!  WAVE  VELOCITY 
120  ! 

130  Young*2. 1E+1 1 

1  AO  Mu*  .  3 

ISO  Rho  =  7800 

160  Cr*  < .87+1 .12*Mu)/<1+Mu)«SQR<Young/2/(1+Mu>/Rho> 

1 70  Cp=SQR<  Young/Rho*( 1 -Hu)/< 1 -2«Mu> ) 

180  Cs=SQR( Young/2/< 1 +Mu) /Rho) 

190  ! 

200  OPTION  BASE  1 

210  DIM  Ve lo< 20 ) 

220  PRINTER  IS  1 

230  INPUT  "ENTER  SCALE  FACTOR  FOR  VELOCITY  -  ".Scvel 

240  INPUT  "SCALE  FACTOR  FOR  TIME,  <1  OR  10)". Sc 

2S0  INPUT  "SCALE  FACTOR  FOR  K.  (1  OR  10)". Sk 

2G0  INPUT  "SCALE  FACTOR  FOR  STRAIN  6000/Scs" .Scs 

270  INPUT  "PLOT  FRAME  7  <Y/N)".Cont$ 

280  INPUT  "PR1N1  DAI  A  '•>  <Y  =  1.  N  -  0>",Pnnvar 

290  PRINT  "KCONT  *  7  <1,  THEN  K<t)-CURVE  IS  PLOTTED;  OTHERWISE  Ep&iy(t>  IS 
PLOTTED  >" 

300  INPUT  Kcont 

310  IF  KcontOI  THEN 

320  PRINT  "WHERE  DO  YOU  WANT  TO  CALCULATE  STRAIN  ?  < INTIAL  CRACK  :  0 

along  crack  path  :  r 

330  INPUT  Strplot 

340  END  IF 

3S0  INPUT  "NUMBER  OF  VELOCITY  AVERAGING  POINTS  N  7",N 
360  ! 

370  * 

380  Step'30  !  CONTROLS  TIME  STEP  IN  THE  CALCULATION 
390  ! 

400  ! 

410  !  GEOMETRIC  PARAMETERS  FOR  SPECIMEN 
420  ! 

430  AO*. 033  •  INITIAL  CRACK  LENGTH 

440  W-8.89E-2  !  WIDTH 

4S0  ! 

460  T0=2.05E-4  !  STRESS  INTENSITY  HISTORY  DURATION 

470  ! 

480  T 1 =  T  0 

490  ImpfacG  .537E+7  !  STRESS  INTENSITY  AMPLITUDE  PARAMETER 
S00  X0=-8,E-4  !  X-COORDINATE  OF  POINT  WHERE  STRAIN  IS  EVALUATED 

510!  YO* . 005  !  Y-COORDINATE  OF  POINT  WHERE  STRAIN  IS  EVALUATED 

520  ! 

530  ! 

540  ! 
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550  VtO'Impfac-Scvel  ?  COMBINED  STRESS  INTENSITY  AMPLITUDE  PARAMETER 
5S0  ! 

570  ! 

580  ! 

590  ! 

600  PRINTER  IS  70S 
610  ! 

620  ! 

630  !  COMPUTATION  OF  POLAR  COORDINATES 
600  X-XO 
650  Y-YO 

660  R-S0R<X'2+Y' 2) 

670  IF  X  =  0  THEN 

680  Theta=PI/2 

690  Thetal =Theta/PI« 1 80 

.30  GOTO  790 

710  END  IF 

720  IF  X>0  THEN 

730  Theta-ATN<Y/X> 

700  Thetal »Theta/PI«180 
750  GOTO  790 
760  END  IF 
770  Theta-PI+ATN<Y/X> 

780  Thetal = Theta/P  I “1  80 
790  GOSUB  2430 
800  ! 

810  !  COMPUTATION  OF  STATIC  COMPLIANCE  AND  STRESS  INTENSITY  FUNCTIONS  FOR 
SPECIMEN  WITH  INITIAL  CRACK  LENGTH 

820  1 

830  FauO-l .09-1 . 735-A0/W+8 . 2«< AO/W)  2-14, 18*<A0/W>'3+14.57*<A0/W>  4 
840  Vau0  = . 76-2 . 28*A0/W+3 ,87*<A0/W)'2-2.O4*(A0/W>'3+.66/(1-A0/W>? 

850  KO=VtO*FawO/VauG/(AO' .5) * ( 1 +SIN< PI/5 > ) 

860  ! 

870  Impulse0=5/ 1 6»T 1  •'VtO/AO/VauO  !  REFERENCE  IMPULSE  PARAMETER 
880  * 

890  1  Kid  IS  THE  TOUGHNESS  OF  THE  MATERIAL  AT  ANY  TIME 
900  * 

910  *  Kidd  IS  THE  INITIATION  TOUGHNESS  OF  THE  MATERIAL 
920  ! 

930  K i dd*5 . 7E+7 
940  ! 

950  PRINTER  IS  1 
960  ! 

970  1  INITIALIZATION  OF  VARIOUS  VARIABLES 
980  ! 

990  A=A0 
1000  TWO 
1010  TOO-TO 
1020  T-0 

1030  T  t - 0 
1040  Vt-VtO 
1050  Ve 1 =0 
1060  Kid=K ldd 
1070  FOR  1=1  TO  N 
1080  Ve lo< I ) -0 
1090  NEXT  I 
1100  Skip'O 
1  1  1  0  Su i tch  =  0 
1120  lmpulse=0 
1130! 

1140! 
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1750  NEXT  J 

1760  Velo< 1 )=Vel 

1  770  Ve  1  =  <  1 -Kid /Kit )*1 . 00*Cr 

1780  IF  Vel<0  THEN  Vel-0 

1790  Vv-Vel 

1800  FOR  J-1  TO  N-1 

1810  Vv*Vv+Velo( J) 

1820  NEXT  J 
1830  Vel-Vv/N 
1840! 

1850' 

1860!  HERE  WE  UPDATE  THE  CRACK  LENGTH  AND  THE  COORDINATES  OF  THE  STRAIN  GAGE 
1870! 

1880  A-A+Vel*1 .E-6/Step 
1890  IF  A/W> . 7  THEN  GOTO  2470 
1900  X=X-Vel«1 .E-6/Step 
1910  R=SQR<X'2+Y  2) 

1820  IF  X=0  THEN 

1930  Theta-PI/2 

1940  The tal = The ta/P I* 1 80 

1950  GOTO  2030 

1960  END  IF 

1970  IF  X>0  THEN 

1980  Theta=ATN(Y/X> 

1990  GOTO  2020 

2000  END  IF 

2010  Theta=PI-‘ATN(Y/X) 

2020  The tal =Theta/PI* 1 80 
2030! 

2040!  HERE  WE  CALCULATE  THE  STRAIN  MEASURED  BY  THE  STRAIN  GAGE  AT  THE  INI  I 1AL 
CRACK  TIP 

2050! 

2060  IF  Kcont=1  THEN  GOTO  2370 
2070  IF  Strplot-1  THEN  GOTO  2140 

2080  H« 17<S0R<SIN<  Theta) >*C0S<Theta/2>*< I +SINtTheta/2>*SIN(3«Theta/2> -Mu« <  1  - 
SIN<  Theta/2)«SIN<3»Theta/2> ) ) ) 

2090  Epsjy*K/Young/H7S0R(2*PI«Y0) 

2100  GOTO  2330 
2110! 

2120!  HERE  WE  CACULATE  THE  STRAIN  ALONG  THE  CRACK  PATH 
2130! 

2140  IF  Ve  1  =0  THEN 

2150  Epsiy-0 

2160  GOTO  2330 

2170  END  IF 

2180  L 1 ■ < 1  - < Ve 1 / Cp >  2)  ‘  ,5 

2190  L2-(1-<Vel/Cs)"2>' .5 

2200  Cl -< 1 +L2‘2) /<  4*L 1 «L2-< 1 +L2'2) ~2) 

2210  C2= 1 +2»L1 ~2-L2"2 

2220  C3=<4*L1*L2)/<1+L2'2> 

2230  C4  =  -< 1 +L2"2) 

2240  R1 =SQR<X2+L 1 ~2-Y02> 

2250  R2-SQR1 X"2+L2' 2»Y0'2> 

2260  Cost  1 cX/R  1 
2270  Cost2“X7R2 
2? 80  Rr“S0R(X'2+Y0'2) 

2290  Theta=ACS(X/Rr  > 

2300  R  z  1  =SQR(C,ost1  +  1  )/SQR(2»R1 ) 

2310  Rz2-SQR<Cost2+1 >/SQR<2»R2> 

2320  Epsiy=C1«(Rz1*(C4-Mu»C2)+Rz2«< l+Mu>«C3)/Young/SaR(2»PI)«K 

2330  IF  Epsiy<0  AND  KcontOI  THEN 
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2890*  PRINT  "PD;PA", I-Sc/Step," 

2900  PRINT  "PD ;PA" , I«Sc/Step 

2910  GOTO  29A0 

2920  IF  Epsly<0  THEN  Epsiy-0 

2930  PRINT  "PD ;PA" , I«Sc/S tep 

29A0  RETURN 

2950  PRINTER  IS  1 

2960  END 
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CRACK  INITIATION  AND  EXTENSION  IN  STEEL  FOR  SHORT  LOADING  TIMES 

J.H.  Giovanola 

SRI  International,  333  Ravensuood  Avenue,  Menlo  Park,  California  HI 023,  U.S.A. 


Beauine'  -  On  a  etudle'  lea  conditions  d1  instability  d'une  fissure  dans  un 
acler  a  haute  resistance  pour  des  temps  de  charge  courts  (aussl  courts  que 
25  pa)  en  utlllsant  un  essal  dynamlque  de  flexion  avec  raise  en  charge  par 
lnertle  (flexion  un  point).  On  trouve  que  pour  le  du&alne  de  vltessea  de 
charge  considere,  un  crltere  de  rupture  incorperant  le  temps  n'eat  pas 
neceasalre.  Par  contre  la  resistance  a  la  ptopagation  peut  croitre  de  faqon 
marquee  avec  l'avance  de  la  fisBure,  en  ralsor.  de  la  formation  de  l>vres  de 
ciaalllement  merac  faibles. 

Abstract  -  The  one-point  bend  impact  test  haB  been  used  to  Investigate  crack 
Instability  conditions  for  short  loading  times  (aa  short  as  25  pa)  in  high 
strength  steel.  It  was  found  that,  in  the  range  of  loading  rates 
investigated,  no  tlme-modif led  fracture  criterion  need  be  invoked.  The 
resistance  to  crack  extension  may,  however,  increase  significantly  with 
crack  growth  due  to  the  formation  of  even  small  shear  lips. 


1  -  INTRODUCTION 

In  recent  years,  considerable  attention  has  been  given  tt  the  prediction  of  crack 
instability  under  l\igh  loading  rates.  It  is  now  well  understood  that,  because  of 
inertia  and  stress  wave  effects,  the  stress  intensity  history  for  dynamically  loaded 
cracks  is  often  significantly  different  from  the  stress  intensity  history 
anticipated  from  quasi-static  considerations  /1-5/. 

Further,  Beveral  experimental  results  suggest  that,  for  very  rapid  loading  (a  few 
microseconds  or  leas),  the  classical  quasi-static  linear  elastic  fracture  mechanics 
criterion  for  crack  instability  must  be  modified  to  account  for  microatructural  rate 
effects  /5-8 /.  On  the  basis  of  experiments  for  very  short  duration  loads  (less  than 
2.2  ps),  Kalthoff  and  Shockey  / 5/  postulated  that  cracks  would  become  unstable  only 
after  the  dynamic  fracture  toughness  has  been  exceeded  for  a  minimum  time-  Baaed  on 
fracture  tests  with  pulse  loads  of  several  durations  (18  to  80  ps),  Homma  et  si.  /6/ 
have  estimated  that,  for  AISI  4340  steel  in  the  Rc  50  condition,  the  minimum  time 
for  fracture  should  be  about  7  ps.  There  is,  however,  some  uncertainty  regarding 
this  value,  because  the  fracture  toughness  reported  in  reference  6  is  only  about  one 
half  th*  value  reported  in  the  literature  for  similar  loading  rates.  More  recently, 
Ravi-Chandar  and  Knauss  /?/  found  that,  when  the  time  required  to  load  a  crack  to 
instability  in  Homalite  100  is  gradually  decreased  below  50  ps,  the  stress  intensity 
value  at  which  the  crack  becomes  unstable  increases  rapidly.  For  a  loading  time  to 
fracture  of  10  ps,  the  stress  Intensity  factor  at  Instability  is  about  twice  the 
value  when  the  loading  time  is  50  pa  or  more.  Kalthoff  /S/  has  rationalited  this 
behavior  by  suggesting  that  a  so-called  Incubation  time  is  necessary  to  nucleate, 
grow,  and  coalesce  microdefects  ahead  of  the  macrocrack.  Conceptually,  the 
incubation  time  is  similar  to  the  minimum  time  postulated  by  Kalthoff  and  Shockey 
/57. 

This  work  attempted  to  establish  whether  fracture  Incubation  time  has  a  significant 
influence  on  crack  instability  conditions  in  high  strength  steel  at  conventional 
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dynamic  loading  rates,  that  Is,  for  fracture  times  down  to  25  ns.  In  evaluating  the 
experimental  crack  Initiation  results,  attention  was  alao  given  to  the  conditions 
controlling  the  first  increments  of  crack  extension. 

II  -  EXPERIMENTAL  PROCEDURES 

Investigation  of  the  crack  instability  conditions  for  high  loading  rates  requires 
detailed  knowledge  of  the  dynamic  stress  intensity  history  experienced  by  the  crack 
tip.  To  obtain  well-characterized  stress  Intensity  histories  In  which  the  loading 
rate  could  be  easily  varied,  we  uaed  a  novel  experimental  technique,  the  one-point 
bend  impact  teat  (1PBT)  /9-11/* 

The  1PBT  uses  a  single-edge-cracked  specimen  and  the  same  testing  arrangement  as  a 
conventional  three-point  bend  impact  test,  except  that  the  end  supports  are 
removed.  When  the  hammer  strikes  Che  specimen,  the  center  portion  of  the  specimen 
la  accelerated  away  from  the  hammer,  whereas  the  end  portions  of  the  specimen  lag 
behind  because  of  inertia.  This  causes  the  specimen  to  bend  and  to  load  the  crack 
tip.  The  resulting  stress  Intensity  histories  do  not  show  the  undesirable 
oscillations  that  often  characterize  the  three-point  bend  Impact  test,  but  rather 
are  smoothly  varying,  as  shown  In  Figure  1  for  two  specimens  with  blunt  notches  and 
with  different  length-over-wldth  ratios.  The  maximum  stress  intensity  amplitude  and 
the  duration  of  loading  are  readily  adjusted  by  changing  the  Impact  velocity  and  the 
specimen  geometry,  respectively.  Varying  these  parameters  in  combination  also 
allows  the  loading  rate  to  be  varied. 


TIME,  ms 

JA-2777-37 


FIGURE  1  STRESS  INTENSITY  HISTORY  FOR  SPECIMEN 
USED  IN  DYNAMIC  FRACTURE  TESTS  I  Impact  Velocity 
Vjfnp  e  5.3  m/s.  Specimen  Dimensions  Thickness  B  =  13  mm. 
Width  Wc  64  mm.  Crack  Length  a  *  19  mm;  Length  L  ■  241  mm 
(full  line  curve)  and  356  mm  (dashed  line  curve)} 


The  atreaa  Intensity  histories  in  Figure  1  w^»re  obtained  by  measuring  the  crack  tip 
strain  intensification  with  a  calibrated  strain  gage  placed  near  the  crack  tip.  The 
measured  strain  was  then  related  to  the  applied  stress  intensity  using  a  calibration 
factor,  determined  in  an  independent  quaal-statlc  three-point  bend  calibration 
experiment  .12/.  This  method  was  also  used  for  all  the  dynamic  fracture  experiments 
reported  in  this  paper. 
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To  investigate  whether  Incubation  tine  affects  the  dynamic  crack  instability 
condition,  we  performed  a  series  of  dynamic  fracture  experiments  at  three  stress 
intensity  rates  (3.3  »  10s,  5  x  105  ,  and  3  x  106  MPa  m1  11  a"1),  using  the  1PBT.  The 
material  tested  was  AISI  4340  aircraft  quality  steel,  quenched  and  tempered  to  a 
Rockwell  hardness  of  50.  The  specimens  used  for  the  tests  at  the  low  and 
Intermediate  stress  intensity  rates  were  356  mm  long,  64  era  wide,  and  13  mm  thick. 

The  ratio  of  crack  length  to  specimen  width  was  0.3  for  the  low  rate  tests  and  0.4 
for  the  intermediate  rate  teats.  The  lower  loading  rates  were  obtained  by  impacting 
the  specimen  at  3.5  and  5.3  m /•  with  a  striker  mounted  on  a  Charpy  pendulum 
machine.  The  specimens  used  for  the  testa  at  the  high  stress  Intensity  rate  were 
241  am  long,  64  mm  wide,  and  13  mm  thick.  The  ratio  of  crack  length  to  specimen 
width  was  0.3.  The  specimens  were  impacted  aL  50  m/a  using  a  compressed  air 
launcher  to  accelerate  a  steel  striker  plate.  The  reference  stress  Intensity 
history  for  these  specimens  is  shown  in  Figure  1.  All  specimens  were  fatigue- 
precracked  according  to  ASTM  standard  E399,  instrumented  with  a  strain  gage  near  the 
crack  tip,  and  calibrated. 

To  complement  the  dynamic  fracture  testa,  the  quasi-static  fracture  toughness  was 
measured  in  triplicate  following  the  standard  ASTM  E399  procedure. 

The  effect  of  Incubation  time  was  tested  as  follows.  In  the  tests  at  the  low  stress 
intensity  rate,  the  specimens  were  impacted  at  a  velocity  sufficient  to  produce  a 
maximum  stress  intensity  slightly  above  the  anticipated  critical  value.  The  loading 
times  to  fracture  were  about  100  pa,  well  above  the  anticipated  incubation  time. 

Thus,  these  testa  provided  a  measure  of  the  dynamic  fracture  toughness,  K1(J.  The 
high  rate  tests  were  designed  to  determine  an  eventual  incubation  time.  In  these 
tests  the  specimens  were  grossly  overloaded  by  increasing  the  impact  velocity.  This 
resulted  in  times  to  fracture  of  20  to  25  pa.  Thus  if  fracture  incubation  time  is 
important  and  of  the  order  of  several  microseconds,  we  expected  fracture  to  occur  at 
significantly  higher  etresa  lntenaity  values  than  the  value  measured  In  the  low  rate 
tests.  The  difference  between  the  time  at  which  the  K1(1  value  was  reached  and  the 
time  of  instability  would  provide  an  estimate  of  the  incubation  time.  The 
intermediate  rate  teats  provided  additional  data  to  monitor  the  influence  of  the 
loading  rate  on  the  value  of  the  stress  intensity  factor  at  crack  initiation. 

111  ~  EXPERIMENTAL  RESULTS 

Figure  2  shows  typical  stress  intensity  histories  measured  at  the  low  and  high 
loading  rateB.  The  point  of  crack  initiation  is  also  indicated.  Since  the  stress 
intensity  factor  is  calculated  from  a  strain  gage  measurement,  the  curves  in  Figure 
2  represent  the  true  ntress  intensity  history  only  up  to  the  point  of  crack 
initiation.  Beyond  thla  point,  the  stress  intensity  at  the  running  crack  tip  is  no 
longer  simply  proportional  to  the  recorded  strain.  Therefore,  the  corresponding 
strain  scale  la  also  given  in  Figure  2. 

At  the  low  and  intermediate  loading  rates,  ersek  initiation  caused  a  clear  drop  in 
the  measured  strain;  therefore,  the  point  at  which  the  stress  intensity  for 
initiation  la  evaluated  can  be  easily  determined.  At  the  high  loading  rate,  the 
strain  history  haH  a  smoother  appearance.  A  simulation  of  the  effect  of  crack 
growth  on  the  strain  record  (discussed  in  the  next  section)  indicates  that, 
depending  on  the  strain  gage  position,  the  loading  rate,  and  the  crack  growth 
behavior,  a  strain  overshoot  above  the  point  of  crack  initiation  is  possible.  Thus 
the  point  of  initiation  was  chosen  as  the  point  of  separation  of  the  recorded  strain 
history  from  a  straight  line  fitted  to  the  main  portion  of  this  curve,  see  Figure 
2(b). 

At  the  low  loading  rate,  the  crack  started  to  propagate  and  then  arrested  after 
tunnelling  2  to  5  mm.  The  arrest  point  can  be  clearly  seen  in  Figure  2(a),  At  the 
intermediate  and  high  rates,  the  specimens  fractured  completely.  In  the  low  and 
intermediate  rate  teats,  the  Jagged  appearance  of  the  strain  records  immediately 
after  intlatlon  suggests  that  crack  extension  occurred  by  a  aeries  of  run-arrest 
events,  Figure  2(a). 
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FIGURE  3  SUMMARY  OF  FRACTURE  TEST  RESULTS 
FOR  FOUR  STRESS  INTENSITY  RATES 


crack  tip-  Tt  is  Implemented  on  a  microcomputer  and  calculates  the  stress  Intensity 
at  regular  small  time  Intervals.  During  each  time  increment,  the  crack  velocity  is 
calculated  using  the  current  stress  intensity  and  the  current  toughness  value 
according  to  a  procedure  proposed  by  Homma  et  al.  / 13/.  The  current  toughness  is  an 
adjustable  variable  of  the  model. 

The  strain  records  obtained  in  the  1PBT  can  be  simulated  by  calculating  the  stress 
intensity  history  with  the  model  and  by  relating  stress  Intensity  and  strain  with 
the  elastic  singularity  solution  /14/.  Ve  studied  the  effect  on  the  records  of  gage 
position  with  respect  to  the  crack  tip  as  well  as  the  Influence  of  different 
initiation  and  propagation  fracture  toughness  assumptions. 

We  examined  two  postulated  toughness  behaviors.  In  one  case,  we  assumed  the 
material  toughness  was  constant  equal  to  the  initiation  toughness  Kjj 

In  the  other  case,  a  toughness  rising  with  crack  extension  was  assumed,  i.e., 

K  »  58  for  Aa  <  0.5  (1) 

K®®*  -  58  -*  14(Aa  -  0.5)  for  0.5  <  Aa  <  1.5 

Kmat  •  72  for  An  >1*5 

mat 

where  is  in  MPa  m1  and  Aa  Is  in  mm. 

Strain  simulations  for  the  low  rate  and  high  rate  experiments  «r^  shown  in  Figure  4 
for  the  two  assumed  toughness  behaviors,  using  the  classical  fracture  criterion  for 
initiation  and  extension 


K  ■  K  . 
mat  appl 

where  K^ppj  is  the  applied  stress  intensity. 


(2) 
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(a)  Low  Rate  Test 


TIME  (ms) 

(b)  High  Rate  Test 

JA-2777-39 

FIGURE  4  SIMULATION  OF  EXPERIMENTAL  STRAIN  RECORDS 
FOR  TWO  TOUGHNESS  ASSUMPTIONS 


The  experimental  strain  records  were  simulated  well  by  the  model,  provided 
expression  (2)  was  used  for  the  toughness,  A  constant  toughness  predicts  a  very 
rapid  strain  drop  at  the  crack  tip  and  an  excessive  arrested  crack  length. 
Expression  (2)  yields  better  agreement  with  the  experimental  strain  profiles;  the 
arrest  crack  length  Is  somewhat  underestimated.  An  optimum  simulation  of  both 
strain  profile  and  arrest  crack  length  was  achieved  when  slightly  different 
resistance  curve  expressions  were  used  for  the  low  and  the  high  rate  experiments. 
Note  that,  as  mentioned  before,  the  strain  In  Figure  4(b)  continues  to  rise  beyond 
the  point  of  crack  Initiation,  Irrespective  of  the  thoughness  assumption. 

V  -  DISCUSSION 

The  results  of  the  dynamic  fracture  tests  show  that,  within  the  range  of  loading 
rates  studied,  the  stress  Intensity  factor  at  Initiation  Is  58.5  KPa  m1  ^  and 
essentially  constant.  This  value  Is  about  927.  of  the  static  toughness  and 
represents  the  dynamic  Initiation  fracture  toughness  Kj^.  This  suggests  that  a 
fracture  Incubation  time  Is  unimportant  for  these  loading  rates.  Further,  our 
results  Imply  that.  If  a  fracture  Incubation  time  exists,  It  Is  much  smaller  than 
the  times  to  frscture  of  20  to  25  ps  considered  here. 
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These  conclusions  are  consistent  with  estimates  of  the  incubation  time  based  on  the 
kinetics  of  aicrostructural  damage.  Mi.crodamage  models  /15,1ft/  based  on  a  viscous 
damage  growth  law  and  calibrated  with  high  rate  experimental  data  indicate  that  the 
characteristic  time  for  the  nucleation  and  growth  of  microvolds  and  microcracks 
should  be  of  the  order  of  a  microsecond.  If  the  incubation  time  is  interpreted  as 
the  time  it  takes  for  microvoids  or  oicrocracks  to  nucleate  and  grow  at  the  crack 
tip,  then  it  should  be  of  the  order  of  the  characteristic  time  for  these  processes. 

Our  conclusions  are  also  consistent  with  recent  results  obtained  for  a  high  strength 
steel  by  Kalthoff  /17 / •  It  was  found  that,  down  to  fracture  times  of  about  21  ps, 
the  stress  intensity  at  initiation  remains  more  or  less  constant.  However,  when  the 
time  t<  fracture  was  decreased  to  about  5  ps,  a  sharp  increase  in  the  stress 
intensity  at  initiation  was  observed. 

The  simulation  of  the  1PBT  predicted  initiation  well  using  the  classical  initiation 
criterion  Kfl  j  ■  Kjj.  However,  to  simulate  the  strain  records  after  initiation  we 
had  to  postulate  a  rising  dynamic  toughness  resistance  curve.  Such  an  assumption 
may  be  Justified,  considering  that  small  shear  lips  (roughly  1  mra  in  width)  were 
formed  during  the  crack  extension  phase  and  that  the  arrested  cracks  in  the  low 
loading  rate  experiments  had  a  thumbnail  shape.  Shear  lips  and  a  bowed  crack  front 
are  known  to  increase  the  apparent  toughness  in  static  fracture  tests.  In 
additional  experiments,  in  which  the  specimens  had  been  side ‘grooved  to  suppress  the 
formation  of  shear  lips,  the  strnin  records  dropped  much  more  abruptly  (as  predicted 
by  the  simulations  with  a  constant  toughness)  than  in  the  experiments  where  shear 
lips  formed,  although  the  value  of  the  stress  intensity  factor  at  initiation 
remained  the  same.  These  latter  results,  together  with  the  simulations  of  the 
tests,  therefore  suggest  chat  the  formation  of  even  small  shear  lips  significantly 
affects  the  strain  history  recorded  after  crack  initiation.  The  shear  lip  acts  as 
an  unbroken  ligament  that  continues  to  load  the  strain  gage  as  the  crack  extends. 

As  a  consequence,  the  shear  lip  may  also  considerably  increase  the  material 
resistance  to  dynamic  crack  extension.  This  effect  needs  to  be  assessed  more 
precisely,  in  particular  to  determine  if  it  is  sensitive  to  crack  velocity. 

VI  -  CONCLUSIONS 

W*  conclude  that,  for  fracture  Limes  in  excess  of  20  ps,  the  crack  initiation 
conditions  for  4340  steel  In  the  Rockwell  hardness  SO  (and  probably  for  other 
similar  high  strength  alloys)  are  not  influenced  by  a  fracture  incubation  time. 
Rather,  our  work  indicates  that  the  initiation  conditions  can  be  determined  solely 
by  using  the  classical  concepts  of  fracture  mechanics.  In  relatively  thin  sections, 
the  resistance  to  dynamic  crack  extension  may  rise  significantly  due  to  the 
formation  of  shear  lips  and  to  crack  tunnelling.  In  these  situations  the  concept  of 
resistance  curve  muat  be  Invoked  to  predict  crack  extension. 
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Appendix  F 


DYNAMIC  MIXED-MODE  CRACK  INITIATION 

INTRODUCTION 

In  numerous  practical  applications,  cracks  are  not  subjected  to  pure 
mode  I  loading,  but  rather  to  a  combination  of  mode  I  and  mode  II 
loading.  Mixed-mode  loading  particularly  occurs  under  dynamic  loading 
conditions  where  loading  from  reflected  waves  or  from  structural 
vibrations  can  be  superimposed  on  the  primary  loads.  In  addition,  the 
degree  of  mixed-mode  loading  strongly  influences  the  crack  propagation 
path.  Therefore,  we  need  to  understand  how  mixed-mode  loading  conditions 
affect  crack  initiation  and  extension.  Although  several  investigations 
have  been  performed  for  crack  initiation  under  mixed-mode  Btatic  loading, 
no  consistent  mixed-mode  crack  initiation  theory  is  yet  available.  For 
the  dynamic  case,  very  few,  if  any,  initiation  data  are  available  to 
evaluate  mixed-mode  crack  initiation  models. 

The  objectives  of  the  work  described  in  thiB  appendix  were  to 
develop  a  testing  technique  capable  of  producing  controlled  dynamic 
mixed-mode  loading  conditions  at  a  stationary  crack  tip  and  to  generate 
dynamic  mixed-mode  Initiation  data  in  order  to  develop  failure 
criterion.  The  significant  progress  made  toward  these  objectives  is 
described  below. 

DEVELOPMENT  OF  A  METHOD  TO  PRODUCE  CONTROLLED  DYNAMIC  MIXED-MODE  LOADING 
CONDITIONS. 

To  produce  controlled  dynamic  mixed-mode  loading  conditions,  we 
chose  to  build  on  our  experience  with  the  inertia-loaded  1PBT  specimen 
and  to  modify  it  to  introduce  a  mode  II  loading  component  at  the  crack 


tip.  Figure  F-l  shows  three  configurations  capable  of  producing  dynamic 
mixed-mode  loading.  The  configuration  in  Figure  F-l(a)  has  been 
investigated  by  Barnes  et  al.^  The  configuration  in  Figure  F-l(b)  is  a 
variation  designed  to  preserve  specimen  symmetry.  Variation  in  the 
mixed-mode  ratio  can  be  achieved  by  changing  the  inclination  of  the 
cracks  with  respect  to  the  impact  direction.  The  configurations  in 
Figure  F-l (a)  and  (b)  both  present  the  disadvantage  that  they  make  it 
difficult  to  introduce  a  controlled  sharp  fatigue  crack.  Therefore,  we 
adopted  the  double-crack  configuration  shown  in  Figure  F-l(c),  which 
allows  fatigue  precracking  in  pure  mode  I  in  an  easy  way. 

For  all  the  work  reported  here  we  chose  specimens  267  mm  long  and 
63.5  mm  wide.  The  mode  1  and  mode  II  stress  intensity  factors  were 
obtained  from  strain  measurements  near  the  crack  tip  using  three  strain 
gages  mounted  at  angles  of  45°,  90° ,  and  135°.  The  method  is  described 
in  detail  in  Appendix  C. 

Double-Crack  1PBT  Configuration 

We  investigated  the  double  crack  specimen  response  by  performing 
calibration  experiments  on  specimens  with  blunt  notches  to  avoid  crack 
initiation.  Figure  2  shows  the  mode  I  and  mode  II  stress  intensity 
histories  for  the  unsupported  double-crack  specimen  later  used  in  the 
fracture  experiments  for  impact  velocities  of  1.7,  2.3,  and  4.3  m  s-^. 

The  crack  length  is  38  mm,  and  the  distance  between  cracks,  2D,  is  102  mm 
as  shown  in  Figure  F-l(c).  Comparison  of  the  stress  intensity  histories 
indicates  that,  as  the  impact  velocity  is  increased,  the  mode  of 
vibration  that  dominates  the  specimen  response  changes.  This  change  is 
evidenced  by  the  shorter  duration  of  the  mode  I  stress  intensity  history 
at  4.3  m  s-^  than  at  2.3  m  s-^  and  by  the  oscillations  appearing  on  the 


^C.  R.  Barnes,  J.  Ahmad,  and  M.  F.  Kanninen,  "Analysis  of  Mixed-Mode 
Crack  Growth  in  a  Ship  Hull  Material,"  Proceedings  of  the  Army  Symposium 
on  Solid  Mechanics,  1982-Critical  Mechanics  Problems  in  Systems  Design, 
AMMRC-MS  82-4  (1982). 
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FIGURE  F-2  STRESS  INTENSITY  HISTORIES  FOR  1PB  DOUBLE-CRACK  SPECIMEN 
WITH  2D  =  102  mm,  a/W  =  O.G 


(a)  Vjmp  =  1 ,7  ms'1 ,  (b)  V,mp  =  2.3  ms'1 ,  (c)  V,mp  =  4.3  ms" 
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mode  II  stress  intensity  history  for  2.3  m  s-1. 
vibrational  mode  change  is  not  yet  understood.  The  later  part  of  the 
stress  intensity  histories  was  not  obtained  at  higher  veloctiy  because  of 
triggering  problems  for  the  test  at  4.3  m  s-1  and  because  the  specimen 
broke  during  the  next  calibration  test  at  12  m  s  The  stress  intensity 
calibrations  of  interest  for  the  fracture  tests  are  those  obtained  for 
impact  velocities  of  4.3  m  s~^  and  higher. 

The  calibration  tests  at  4.3  and  12  m  s-1  demonstrated  that,  once 
the  higher  vibrational  inodes  dominate  the  specimen  response,  no  further 
change  in  the  shape  of  the  stress  intensity  history  occurs.  Further,  the 
at  plitudes  of  the  stress  intensity  history  scale  proportionally  to  the 
impact  velocity  as  is  the  case  in  the  single  crack  1PBT.  Figure  F-2(c) 
shows  that,  on  impact,  K-j-j  first  goes  through  a  maximum  and  then  drops 
back  to  essentially  zero.  K^,  however,  becomes  slightly  negative  at 
first  (compressive  load  on  the  crack  tip)  and  then  has  a  sinusoidal 
history.  Thus,  Kj  is  either  negative  or  small  during  the  major  portion 
of  the  Kjj  history. 

We  tested  additional  specimens  to  investigate  how  changing  the  crack 
length  and  the  crack  spacing  affects  the  stress  intensity  histories.  For 
these  additional  calibration  experiments  no  changes  in  the  dominant 
vibrational  mode  were  observed  for  impact  velocities  up  to  4.3  m  s-^,  and 
we  presume  that  no  change  will  occur  for  higher  velocities. 

Figure  F-3  shows  the  effect  of  varying  the  distance  2D  between 
cracks  from  50  to  100  mm  keeping  the  crack  length  constant  at  19  mm. 

For  the  Kj  histories  in  Figure  F-3(a),  the  time  scales  have  been  shifted 
by  about  20  pa  with  respect  to  each  other,  and  for  the  histories  in 
Figure  F-3(b),  they  have  been  shifted  by  about  8  ps  to  account  for  the 
difference  in  arrival  time  at  the  crack  tip  of  the  flexural  and  shear 
waves,  respectively,  caused  by  che  difference  in  crack  spacing. 

The  early  part  of  the  Kj  histories  is  the  same  for  both  crack 
positions,  except  that  Kj  becomes  slightly  negative  at  first  for  the 
specimen  with  the  larger  value  of  2D.  They  start  diverging  once  the 
stress  waves  reflected  at  the  far  end  of  the  specimen  reach  the  crack 


F-5 


farthest  from  the  specimen  mid-plane.  The  history  for  this  crack  has 
a  shorter  oscillation  duration  and  a  lower  maximum  amplitude  than  the  one 
for  the  crack  closest  to  the  crack  plane. 

The  Kj-j  histories  in  Figure  F-3(b)  are  essentially  the  same  for  both 
crack  spacings  at  early  times,  with  small  differences  at  later  times  when 
the  Kj  histories  also  start  to  diverge.  In  the  two  cases,  the  mode  II 
stress  intensity  amplitudes  are  small. 

Figure  F-4  illustrates  the  influence  of  varying  the  crack  length  on 
the  stress  intensity  histories,  keeping  the  crack  spacing  constant  (50 
mm).  In  this  Figure,  the  time  origin  corresponds  to  the  time  of  impact, 
and  the  stress  intensity  histories  have  not  been  shifted  with  respect  to 
each  other.  For  reasons  not  yet  understood,  the  mode  I  stress  intensity 
for  the  shorter  crack  rises  earlier  than  the  stress  intensity  for  the 
longer  crack,  although  the  initial  slope  is  the  same.  Both  the  duration 
and  the  maximum  amplitude  of  the  Kj  history  increase  with  crack  length. 
The  shape  of  the  early  times  of  the  Kjj  histories  is  similar  for  both 
crack  lengths,  with  the  amplitude  increasing  with  crack  length.  At  later 
time  the  K-q  histories  differ,  probably  because  of  the  different  flexural 
behaviors . 

The  calibration  results  presented  above  demonstrate  that,  although 
some  features  of  the  mode  I  stress  intensity  histories  are  anticipated 
based  on  the  results  of  the  single  crack  1PBT,  the  dynamic  response  of 
the  double-crack  test  is  more  complicated  and  will  require  further 
studies  to  understand  the  coupling  between  shear  and  flexural 
deformation.  The  calibration  results  also  illustrate  an  important 
characteristic  of  dynamic  mixed-mode  loading,  namely,  that  the  mixed-mode 
ratio  Kjj/Kj  varies  continuously,  but  not  necessarily  raonotonicaly  with 
time.  This  variation  makes  the  interpretation  of  dynamic  mixed-mode 
initiation  tests  more  complicated,  because  possible  loading  history 
effects  must  be  taken  into  account.  However,  the  variation  of  the  Kjj/Kj 
ratio  with  time  offers,  in  principle,  a  simple  mean  of  varying  the  degree 
of  mixed-mode  at  initiation  by  varying  the  impact  velocity. 
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FIGURE  F-4  STRESS  INTENSITY  HISTORIES  FOR  DOUBLE-CRACK  1 PB  SPECIMENS 
WITH  DIFFERENT  CRACK  LENGTHS  (2D  =  50  mm  V  =4.3  ms'1) 
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The  calibration  result?  of  Figure  2  further  show  that  a  significant 
mcde  II  loading  component  can  be  produced  with  the  double-crack  one- 
point-bend.  With  impact  velocities  of  30  m  a-1,  maximum  amplitudes  of  54 
MPa  m^^  have  been  achieved  In  mixed-mode  fracture  experiments.  These 
values,  however,  proved  too  small  to  Initiate  a  crack  under  mainly  mode 
II  conditions,  and  there  is  a  limit  to  raising  the  impact  velocity  to 
obtain  higher  stress  intensity  values  because,  beyond  a  certain  velocity, 
a  large  plastic  indent  develops  under  the  hammer.  Moreover,  the  maximum 
K-j-j-  amplitude  occurs  when  Kj  is  either  compressive  or  zero.  In  turn,  the 
mode  II  loading  decreases  to  zero  while  the  mode  I  loading  becomes 
significant.  Thus  the  test  configuration  in  Figure  F-l(c)  produces 
loading  conditions  in  which  the  mode  I  and  the  mode  II  loading  are 
essentialy  decoupled. 

The  configuration  of  Figure  F-l(c)  must  be  modified  to  increase  the 
maximum  achievable  amplitude  and  to  obtain  larger  ratios  at 

stress  intensity  levels  close  to  the  material  mode  I  and  mode  II 
toughnesses  and  for  reasonable  impact  velocities. 

Modified  Double-Crack  1PB  Configuration 

To  achieve  higher  mode  II  stress  intensity  amplitudes,  at  least 
three  solutions  are  possible.  The  first  solution  is  to  increase  the  crack 
length,  keeping  the  specimen  width  constant.  A  second  solution  is  to  use 
a  wider  hammer  p?ate  to  minimize  the  plastic  indentation  and  permit 
higher  impact  velocities.  We  also  anticipate  that  this  solution  will 
Induce  a  sharper  acceleration  of  a  larger  portion  of  the  center  part  of 
the  specimen  and  hence  a  higher  shear  force  on  the  uncracked  ligamer.t  and 
a  higher  K-q  amplitude.  The  third  solution  Is  to  increase  the  specimen 
Inertia  between  the  crack  and  the  far  end  uf  the  specimen  either  by 
Increasing  the  thickness  of  th-2  specimen  in  this  area  or  by  adding 
ballast  plates.  These  solutions  have  net  been  tested  in  this 
Investigation. 


To  achieve  higher  mixed-mode  ratios  at  higher  stress  intensity 
levels,  at  least  two  solutions  are  feasible-  The  first  solution  is  to 


control  the  amount  of  mode  I  loading  by  statically  prestressing  the 
specimen.  This  solution  presents  the  additional  advantage  of 
eliminating  the  compressive  mode  I  loading  component,  which  complicates 
interpretation  of  the  fracture  experiments. 

The  second  solution  is  to  introduce  a  delayed  shear  loading  at  the 
time  in  the  loading  history  when  the  mode  1  loading  becomes  important. 
This  delay  can  be  achieved  by  using  unrestrained  end  supports  that  are 
free  to  move  with  the  specimens  once  they  are  impacted,  as  shown  in 
Figure  F-5(a).  At  early  rimes  after  the  hammer  impact,  the  double-crack 
specimen  does  not  interact  with  the  supports  and  the  stress  intensity 
histories  are  similar  to  the  histories  shown  in  Figure  F-2.  Later, 
however,  the  specimen  comes  in  contact  with  the  unrestrained  supports, 
which  transmit  renewed  inertial  shear  loading  on  the  uncracked  ligament. 
Although  we  have  not  investigated  this  point,  the  amplitude  of  the  second 
shear  loading  could  be  adjusted  to  a  certain  extent  by  varying  the 
support  position  with  respect  to  the  crack  position 

To  perform  our  mixed-mode  fracture  experiments,  we  used  a 
combination  of  static  prestressing  and  additional  inertial  shear  loading 
by  unrestrained  supports.  The  experimental  arrangemunt  is  shown  in 
Figure  F-5(b).  The  specimen  is  loaded  statically  in  three-point-bend  by 
means  of  a  support  plate  and  a  yoke.  The  amount  of  prestressing  can  be 
adjusted  with  the  screw  positioned  opposite  the  impact  point.  The  whole 
prestressed  specimen  assembly  is  supported  in  the  air  launcher  apparatus 
used  for  the  experiments  in  a  way  that  permits  free  movement  in  the 
direction  parallel  to  the  impact  direction.  The  effect  of  the 
prestressing  and  of  the  prestreasing  fixture  on  the  stress  intensity 
histories  is  illustrasted  in  the  results  of  the  fracture  experiments 
presented  below. 

DYNAMIC  MIXED-MODE  EXPERIMENTS 

We  used  the  experimental  arrangement  described  in  the  previous 
paragraph  to  perform  dynamic  mixed-mode  experiments.  One  of  the  two 
notches  of  the  specimen  was  fatigue-precracked ,  whereas  the  other  had  a 
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FIGURE  F-5  ARRANGEMENTS  TO  CONTROL  MIXED-MODE  RATIO 
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(b)  Combined  prestress  and  unrestrained  support  arrangemen 


root  radius  of  approximately  0.125  mm.  The  total  length  of  the  sharp 

crack  was  the  game  as  the  length  of  the  blunt  notch,  38  mm.  The  material 

in  these  tests  was  the  same  as  in  the  rest  of  the  program,  aircraft 

quality  4340  steel,  quenched  and  tempered  to  HRC  50.  The  amount  of 

prestress  measured  in  terms  of  the  mode  I  stress-intensity  factor  was 
S£d  t 

varied  from  K-  ■  0  to  a  level  equal  to  the  maximum  amplitude  applied 
during  fatigue  precracking,  i.e.,  60%  of  the  dynamic  mode  1  fracture 
toughness,  58  MPa 

Figure  F-6  shows  the  stress  intensity  histories  for  three  tests  with 
different  levels  of  prestress  (specimens  MMF-2,  MMF-6,  and  MMF-8) .  These 
histories  were  calculated  from  the  records  of  three  strain  gages  placed 
near  the  crack  tip  (see  Appendix  C).  The  point  of  crack  initiation  is 
defined,  by  analogy  with  the  pure  mode  I  dynamic  fracture  initiation,  as 
the  point  where  the  mode  I  stress  Intensity  history  departs  significantly 
from  the  tangent  to  the  rising  portion  of  the  curve.  The  impact 
velocities  for  each  test  are  indicated  in  the  figure.  Tests  performed 
under  essentially  the  same  conditions  demonstrated  that  the  fracture 
initiation  conditions  as  well  as  the  crack  extension  path  are 
reproducible. 

The  results  of  the  dynamic  mixed-mode  initiation  experiments  are 
summarized  in  Table  F-l  and  Figure  F— 7.  Figure  F-7  plots  Kjj,  normalized 
by  the  dynamic  initiation  toughness  KId,  as  a  function  of  Kj  also 
normalized  by  Kjd»  Solid  circles  represent  conditions  under  which  crack 
initiation  occurred;  open  circles  represent  mixed-mode  loading  states 
that  were  traversed  during  the  dynamic  experiments  without  causing  the 
crack  to  extend.  Three  loading  paths,  corresponding  to  the  stress 
intensity  histories  in  Figure  F-6,  are  indicated  by  the  dashed  lines. 

The  open  triangles  represent  initial  prestress  conditions  used  in  the 
test3,  and  the  single  solid  triangle  indicates  a  value  of  the  static  mode 
II  initiation  toughness  for  4340  steel  (HRC  50),  reported  by  Shah.^ 


^R.  C.  Shah,  "Fracture  under  Combined  Modes  in  4340  Steel,"  ASTM-STP  560, 
pp.  29-52  (1974). 
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Table  F-l 

RESULTS  OF  MIXED-MODE  FRACTURE  EXPERIMENTS 
acture  At 


Test  No. 


Configuration 


Kt  [MPa  b^^““ktt  IMP8  ml/2] 


Kr/K 


„Max 


Lilli 


— LL» 


F--13 


L  -t-T  If.'w-W  w* ->-'.•/ ••.aV.VjV.'  1 


La  jfj  ■£_ 


Angle  of 
Crack  Extension 


MMF-l 

No  Prestreae 

54.3 

11.0 

0.2 

44.3 

49” 

MMF-2 

No  Preatresa 

56 

0 

0 

49.0 

10” 

MMF-3 

No  Preatresa 

54 

3 

0.06 

34.5 

3°  -  9° 

MMF-4 

No  Prestress 

63.4 

32.4 

0.51 

41.4 

35° 

MMF-5 

No  Preatresa 

63.5 

27.6 

0.43 

35.9 

35” 

MMF-6 

No  Prestresa 

53.8 

18 

0.34 

25.5 

18”  ♦  5“ 

MMF-7 

No  PrestreBs 

98 

31 

0.32 

53.8 

23°  -  33” 

MMF-8 

No  Preatresa 

82.7 

31.7 

0.38 

52  .4 

41“  +  59” 
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The  shaded  band  In  Figure  F-7  is  a  tentative  lower  bound  mixed-mode 
failure  curve  that  will  be  discussed  below. 

Figure  F-7  presents  dynamic  initiation  data  for  mixed-mode  ratios 
Kn/K-j  up  to  0.5.  Most  of  the  initiation  data  points  fall  around  a 
vertical  line  -  1,  except  for  the  two  points  on  the  far  right  of 

the  figure.  This  results  suggests  that,  under  relatively  mild  mixed-mode 
conditions,  dynamic  crack  initiation  is  essentially  controlled  by  the 
mode  I  stress  intensity  amplitude.  The  two  failures  that  occurred  at  Kj 
values  significantly  larger  than  the  dynamic  mode  I  fracture  toughness 
(specimens  MMF-7  and  MMF-8)  represent  specimens  that  were  tested  at 
higher  impact  velocity  than  other  specimens  with  comparable  mixed-mode 
ratios.  The  significance  of  these  two  data  points  will  be  discussed 
below. 

Figure  F-8  illustrates  the  crack  path  for  each  of  the  specimens 
tested,  except  of  specimens  MMF-1  and  MMF-7,  which  had  a  rough  and  uneven 
crack  path.  With  the  exception  of  specimen  MMF-8,  the  initial  crack 
propagation  angle  appears  to  increase  with  the  mixed-mode  ratio  at 
initiation.  The  propagation  angle  changed  rapidly  in  some  specimens, 
reflecting  the  changing  mixed-mode  condition  at  the  extending  crack  tip. 

For  most  specimens,  the  macroscopic  fracture  surface  appearance  was 
very  similar  to  the  fracture  surface  appearance  in  pure  mode  I 
experiments.  However,  in  experiment  MMF-7  and  particularly  in  experiment 
MMF-8,  the  fracture  surface  has  a  much  rougher  appearance  than  in  the 
other  specimens.  Preliminary  scanning  electron  microscope  observations 
confirm  these  differences. 

Figure  F-9  compares  the  fracture  surface  of  specimen  MMF-8,  Figure 
F-9(a)  and  (b),  with  the  fracture  surface  of  specimen  MMF-5,  Figure  F- 
9(c).  Specimen  MMF-8  has  a  narrow  region  near  the  fatigue  crack  tip  of 
relatively  smooth  appearance,  which  corresponds  to  a  region  of  failure  by 
ductile  growth  of  small  voids.  Beyond  this  narrow  band,  the  fracture 
surface  is  very  rough  and  displays  large  inclined  steps  Figure  F-9(a).  A 
view  of  this  area  at  higher  magnification.  Figure  F-9(b) ,  suggests  that 
the  Inclined  steps  may  correspond  to  grains  that  have  pulled  out  along 
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MMF-3  Ky/Kj  =  0.06 


MMF-2  Ky/K,  =  0 


MMF-6  Kn/Kj=0.34 


MMF-4  Kjj/Kj  =  0.51 


MMF-5  Kn/Kj  =  0.43 


MMF-8  Ky/K,  =  0.38 


RA-1  750-82 


FIGURE  F-e  CRACK  PATH  IN  MIXED-MODE  FRACTURE  SPECIMENS 
(Arrows  indicate  initial  crack  tip  position.) 
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grain  boundaries*  In  contrast,  such  steps  are  not  visible  on  the 
fracture  surface  of  specimen  MMF— 5,  except  for  a  few  limited  occurrences. 
Figure  F-9(c).  Note,  however,  that  the  difference  in  fracture  surface 
appearance  arises  only  about  1  mm  away  from  the  initial  fatigue  crack 
front;  therefore,  it  can  be  argued  that  the  fracture  surface  appearance 
reflects  differences  in  crack  propagation  behavior  and  not  difference  in 
crack  initiation  behavior. 

DISCUSSION  OF  DYNAMIC  MIXED-MODE  FRACTURE  EXPERIMENTS  AND  CONCLUSIONS 

The  dynamic  mixed-mode  fracture  data  presented  here  represent  the 
first  set  of  results  obtained  with  a  new  experimental  technique; 
therefore,  they  must  be  analyzed  with  caution  before  conclusions  can  be 
drawn.  One  of  the  limitations  of  the  results  is  the  small  number  of 
specimens  tested  and  the  restricted  range  of  mixed-mode  ratios 
investigated.  A  more  serious  limitation  is  that  the  strain  gage 
technique  for  measuring  the  Kj  and  K1X  histories  at  the  stationary  crack 
tip  has  not  been  assessed  against  another  experimental  technique  or 
against  numerical  simulations.  Although  the  stress  intensity  histories 
measured  with  the  strain  gage  technique  qualitatively  agree  with 
histories  anticipated  from  simple  specimen  dynamic  considerations,  an 
assessment  against  other  methods  is  necessary  to  gain  full  confidence  in 
the  strain  gage  technique. 

Another  limitation  of  the  approach  taken  in  this  investigation  is 
the  detection  of  crack  initiation  under  mixed-mode  conditions.  As  for 
the  pure  mode  I  experiments,  we  have  determined  the  instant  of  crack 
initiation  based  on  the  strain  history  measured  by  the  gage  mounted  90° 
from  the  crack  plane  and  immediately  above  the  crack  tip.  A  sharp  drop 
in  the  strain  signal  recorded  by  this  gage  has  been  interpreted  as  an 
indication  of  crack  initiation.  However,  in  contrast  to  the  pure  mode  I 
case,  we  have  not  yet  investigated  how  crack  extension  under  mixed-mode 
conditions  affects  the  strain  signals  after  initiation.  For  small  mixed- 
mode  ratios,  the  initial  angle  of  crack  extension  is  small,  and  crack 
extension  will  probably  affect  the  strain  histories  in  a  manner  similar 


to  the  pure  mode  I  case.  As  the  degree  of  mixed-mode  and  the  initial 
angle  of  crack  extension  increase,  the  crack  front  will  tend  to  move 
toward,  rather  than  away  from,  the  strain  gages.  Immediately  after 
initiation,  this  may  at  first  cause  a  rise  in  the  strain  measured  by  the 
90°  gage.  As  a  consequence,  the  time  of  crack  initiation  and  the  value 
of  Kj  would  be  overpredicted.  This  argument  seems  borne  out  by  our 
experimental  data,  which  show  that  the  K-j-  values  at  initiation  increase 
with  the  initial  angle  of  crack  extension.  In  particular,  the  very  high 
values  of  Kj  measured  for  tests  MMF-7  and  MMF-8  maybe  due  to  this 
extension  effect.  Therefore,  we  strongly  suspect  that  the  data  points  on 
the  far  right  in  Figure  F-7  should  be  moved  toward  the  line  %/Kj^  -  1, 
as  indicated  by  the  small  arrows  in  the  figure. 

A  final  comment  on  the  mixed-mode  experiments  concerns  the  number 
of  parameters  that  must  to  be  considered  when  interpreting  the  fracture 
results.  At  least  three  factors  may  influence  the  dynamic  mixed-mode 
initiation  results:  the  loading  history,  the  mixed-mode  ratio  at 
initiation,  and  the  loading  rate.  With  the  test  method  we  used  (as 
probably  with  any  other  dynamic  mixed-mode  test  method) ,  it  is  not 
possible  to  vary  these  three  factors  independently.  This  limitation  may 
make  it  difficult  to  make  an  unambiguous  interpretation  of  the  results. 

Bearing  in  mind  the  above  limitations,  we  can  draw  the  following 
conclusions  from  the  present  dynamic  mixed-mode  fracture  investigation. 

(1)  We  have  demonstrated  the  applicability  of  a  simple  experimental 
technique  for  obtaining  dynamic  mixed-mode  crack  initiation 
data.  The  technique  permits  us  to  initiate  cracks  with  dynamic, 
mixed-mode  ratios  Kjj/K-j-  varying  between  0  and  0.5.  The 
dependence  of  the  stress  intensity  history  on  crack  length  and 
position  and  on  impact  velocity  has  been  established.  From 
these  results,  we  have  proposed  modifications  to  achieve  mixed- 
mode  ratios  at  crack  initiation  larger  than  0.5. 

(2)  The  dynamic  mixed-mode  crack  initiation  data  generated  for  the 
range  Kj-j/Kj  =  0  to  0.5  show  that,  for  this  range  of  mixed  mode 
ratios,  crack  initiation  is  essentially  controlled  by  the  mode  I 
stress  intensity  and  that  it  occurs  when  K1  =  KId’  the  pure  mode 
I  initiation  toughness. 


(3)  Based  on  the  points  of  the  Kj-Kjj  plane  traversed  during  the 

fracture  experiments  for  which  no  failure  occurred,  we  can  draw 
a  tentative  lower  bound  mixed-mode  failure  curve  for  4340  steel 


HRC  50  (Figure  F-7) .  The  pure  mode  II  dynamic  toughness  K^Id  *s 
at  least  equal  to  the  pure  mode  I  initation  toughness  hi* 

Kild  >  58  MPa  m1'^.  Assuming  that  KIId  has  the  same  rate 
sensitivity  as  Kj^,  which  is  very  low,  and  taking  for  the  static 
mode  II  fracture  toughness  the  value  reported  by  Shah,^  88  MPa 
nr-'  ,  we  estimate  that  KIId  should  be  about  85  MPa  m  '  •  This 
value  is  consistent  with  the  proposed  lower  bound  failure  curve. 


(4)  We  have  observed  a  significant  difference  in  the  fracture 

surface  appearance  both  macroscopically  and  microscopically,  in 
the  specimens  with  the  largest  overall  crack  propagation  angle. 
These  specimens  are  also  the  specimens  that  were  tested  at  the 
highest  impact  velocities.  This  result  suggests  that  the 
micromechanism  of  dynamic  mixed-mode  crack  extension  is  affected 
either  by  the  degree  of  mixed-mode  or  by  the  crack  velocity  or 
by  both  factors. 


For  future  studies ,  we  recommend  that  the  dynamics  of  the  double¬ 
crack  specimen  configuration  be  further  investigated  to  answer  some  of 
the  questions  in  this  study.  The  methods  proposed  to  increase  the 
achievable  mode  II  component  of  loading  should  be  tried,  and  a  procedure 
for  a  more  reliable  detection  of  crack  initiation  should  be  developed. 
Simply  placing  the  strain  gages  on  the  exterior  face  of  the  crack  may 
prove  adequate;  with  this  arrangement,  the  extending  crack  will  move  away 
from  the  gages  instead  of  toward  them.  More  mixed-mode  initiation  data 
should  then  be  generated  with  the  improved  experimental  proceedure  to 
augment  and  verify  the  failure  curve  obtained  in  this  study. 


F-21 


Appendix  G 


INVESTIGATION  OF  THE  INITIAL  PHASE  OF  DYNAMIC  CRACK 
EXTENSION  FOR  IMPACT  LOADING  CONDITIONS 


INTRODUCTION 

The  relationship  between  dynamic  initiation  and  propagation 
toughnesses  is  an  important  aspect  of  dynamic  fracture  from  the 
fundamental  as  well  as  the  practical  point  of  view.  Initiation  and 
propagation  are  not  two  separate  processes  but  two  phases  of  one 
continuous  process  controlled  by  the  same  microdamage  mechanisms. 
Therefore,  from  a  fundamental  point  of  view,  a  complete  understanding  of 
dynamic  fracture  will  be  achieved  only  if  a  unique  model  can  be  developed 
that  consistently  explains  experimental  observations  for  both  dynamic 
crack  initiation  and  propagation.  In  addition,  interpretation  of  dynamic 
initiation  behavior  at  very  high  loading  rates  may  also  hinge  on 
understanding  the  initial  phase  of  crack  extension.  From  a  practical 
point  of  view,  an  understanding  of  the  relationship  between  dynamic 
initiation  and  propagation  toughnesses  may  open  the  possibility  for 
obtaining  estimates  of  the  propagation  toughness  from  values  of  the 
dynamic  initiation  toughnesses,  which  are  significantly  easier  to 
measure. 

Few  studies  have  focused  on  the  initial  part  of  dynamic  crack 
extension  immediately  following  initiation. ^  Most  dynamic  crack 
propagation  studies  have  been  conducted  on  cracks  initiated  at  blunt 
notches  under  quasi-static  loading  conditions.  In  many  practical  cases, 


■^K.  Ravi-Chandar  and  W.  G.  Knauss,  "An  Experimental  Investigation  into 
Dynamic  Fracture:  I.  Crack  Initiation  and  Arrest,"  Int.  J.  Fract.,  Vol. 
25,  pp.  247-262  (1984). 
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iwever,  running  cracks  will  initiate  from  sharp  precursor  fatigue  cracks 
id  the  loading  conditions  will  be  dynamic.  This  latter  situation 
ffers  significantly  from  the  situation  of  a  crack  initiated  from  a 
itch  under  quasi-static  loading.  When  the  crack  initiates  at  a  blunt 
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MATERIAL  INVESTIGATION  AND  EXPERIMENTAL  PROCEDURE 

The  material  used  in  this  investigation  was  AISI  aircraft  quality 
4340  steel  quenched  and  tempered  to  HRC50.  Two  batches  of  material  were 
studied:  9.5-mm-thick  plates  and  12.7-mm-thick  plates. 

For  each  batch  of  material,  the  static  and  the  dynamic  fracture 
toughnesses  were  measured,  the  latter  by  means  of  the  1PBT.  The  results 
are  summarized  in  Table  G-l.  The  point  of  dynamic  crack  initiation  was 
determined  from  the  near  crack  tip  strain  records  by  means  of  the 
tangent-to-the-loading-curve  method  discussed  in  Appendices  B  and  C.  The 
results  for  each  of  the  two  thicknesses  of  material  will  be  discussed 
separately.  More  data  were  obtained  for  the  9.5-mm-thick  material;  the 
data  from  the  12.7-mm-thick  material  were  used  as  a  complement  to  assess 
the  effect  of  the  specimen  thickness  on  the  results. 

Table  G-l 


STATIC  AND  DYNAMIC  FRACTURE  TOUGHNESS 


Specimen  Thickness  (ram) 

Ktj  (MPa  m1/2) 

K1C  (MPa  m1/2) 

9.5 

63.8  ±  6 

59.9  ±  5 

12.7 

58.5  ±  5 

63.7 

Three  methods  were  used  to  estimate  the  dynamic  propagation 
toughness  and  to  relate  it  to  the  initiation  toughness:  crack  .lump 
experiments,  strain  record  analysis,  and  direct  strain  gage 
measurements.  In  the  first  method  we  performed  a  series  of  experiments 
in  which  the  impact  velocity  was  gradually  increased  to  produce 
controlled  amounts  of  crack  extension.  These  experiments  are  termed 
crack  jump  experiments  (CJ).  By  simulating  these  experiments  with  the  1 
PBT  model  described  in  Appendix  D,  we  assessed  the  crack  extension 
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behavior  over  short  propagation  distances  after  initiation  The 

specimens  used  in  the  CJ  experiments  were  98.5  mm  long,  38.1  mm  wide,  9.5 

ram  thick.,  and  had  a  crack  length  over  width  ratio  of  0.3.  The  stress 

Intensity  history  for  thus  specimen  geometry  and  a  nonpropagating  crack 

is  shown  in  Figure  G-l  for  an  impact  velocity  of  5.3  m  s-^.  It  is  seen 

that  a  stress  intensity  pulse  shorter  than  90  ps  can  be  obtained;  this 

short  pulse  combined  with  varying  the  impact  velocity  (and  hence  the 
Max 

amplitude  )  permits  a  precise  control  of  the  crack  jump  distance.  In 
the  CJ  experiments,  some  of  the  specimens  were  subjected  to  a  single 
impact,  each  specimen  at  a  different  impact  velocity,  to  produce  a 
different  amount  of  crack  extension.  In  these  experiments,  crack 
extension  was  detected  by  inspecting  the  near  crack  tip  strain  records. 

As  the  crack  jump  distance  decreases,  it  becomes  increasingly 
difficult  to  detect  crack  initiation  from  the  strain  records.  To  attempt 
to  circumvent  this  difficulty,  we  subjected  the  remainder  of  the 
specimens  to  multiple  impacts,  starting  at  low  velocity  and  gradually 
increasing  it  to  bracket  as  closely  as  possible  crack  initiation.  After 
each  impact,  the  specimens  were  inspected  for  signs  of  crack  initiation 
by  a  compliance  measurement  method  or  by  replicating  the  crack  tip  region 
with  acetate  tape  and  observing  the  replica  with  a  light  microscope. 
Although  neither  technique  gave  satisfactory  results,  we  were  able  to 
produce  crack  extension  increments  as  small  as  50  pm  in  the  multiple 
impact  experiments. 

In  addition  to  the  small  98.5-mm  specimens,  one  larger  specimen  (203 
by  98  by  9.5  mm,  a/W  =  0.35)  was  also  tested  to  determine  the  influence 
of  the  specimen  size  and  consequently  the  stress  intensity  pulse  duration 
on  che  crack  jump  behavior. 

Ten  CJ  experiments  were  performed.  After  the  tests,  the  crack  jump 
specimens  were  subjected  to  low  amplitude  fatigue  loading  to  mark  the 
position  of  the  arrested  crack  front,  and  they  were  broken  open.  Crack 
extension  was  then  measured  on  the  fracture  surface  using  scanning 
electron  or  optical  microscopy. 


the  stress  intensity  at  the  propagating  crack  tip.  However,  if  the  crack 
jump  distance  is  much  smaller  than  the  specimen  width,  the  stress 
intensity  history  for  the  specimen  with  the  stationary  crack  provides  an 
estimate  of  the  stress  intensity  history  experienced  by  the  extending 
crack. 


Figure  G-2  shows  that,  initially,  a  large  increase  in  *  ant* 

hence  in  impact  velocity  is  required  to  produce  a  small  increase  in  crack 

jump  distance.  Beyond  a  Kapp^  value  of  about  80  MPa  the  trend  is 

IMax 

reversed  and  small  increases  of  Kapp  above  80  MPa  m1''2  cause 

■Max 

comparatively  large  crack  jumps.  The  smallest  crack  extension  detected 

in  the  Cl  experiments  was  50  pm,  for  a  KapP^  of  58  MPa  m^2.  From  this 

IMax 

observation,  we  conclude  that  crack  initiation  must  occur  around  57  MPa 


Note  that  data  points  for  specimens  that  have  undergone  multiple 
impacts  fall  on  the  same  curve  as  specimens  impacted  only  once. 

Figure  G-2  Is  reminiscent  of  the  resistance  curves  of  static 
fracture  mechanics,  which  plot  the  Increase  in  toughness  as  a  function  of 
crack  extension.  However,  it  could  be  argued  that  the  shape  of  the  curve 
in  Figure  G-2  Is  simply  a  result  of  the  loading  history  at  the  crack 
tip.  That  is,  for  KapJ^  values  only  slightly  over  the  value  of  the 
initiation  toughness,  the  crack  is  loaded  only  long  enough  to  propagate  a 
very  small  distance  before  It  is  unloaded  by  the  one-point-bend  specimen 
bending  back  after  reaching  its  maximum  inertial  deflection.  If  this  was 
the  case,  changing  the  stress  intensity  duration,  by  changing  the 
specimen  size,  and  keeping  Kapp^  the  same  would  change  the  resulting 
crack  jump  distance.  Figure  G-2  shows  that  this  is  not  the  case  because 
the  data  point  for  a  geometrically  similar  specimen  but  scaled  up  in  size 
by  a  factor  of  2.3  falls  on  the  same  curve  as  the  data  points  for  smaller 
specimens.  Therefore  Figure  G-2  suggests  that  the  propagation  toughness 
increases  significantly  above  the  Initiation  toughness  level  as  the  crack 
extends  a  short  distance. 

This  conclusion  is  further  supported  by  predictions  of  the  crack 
jump  distance  as  a  function  of  the  Kapp'''  obtained  with  the  simple  1PBT 


model,  assuming  two  behaviors  for  the  propagation  toughness.  In  one 
case,  the  propagation  toughness  was  assumed  constant  equal  to  the  Initial 
toughness,  KjD  »  -  57  MPa  m^2.  In  the  other  case,  was  assumed 

to  depend  both  on  the  crack,  velccity  a  and  the  amount  of  crack  extension 

Aa  and  to  have  the  form 

KID  "  57  1C>6  +  526  Aa>  +  f<a)  0  <  Aa  <  0.5  mm 

(G-l) 

KID  "  72  1()6  (1  +  89  <Aa  "  0,5  10"3^  +  0,5  ™i  <  Aa  <  2mm 

K  =  81.6  106  +  f(a)  Aa  >  2  mm 

where 

f(£)  =■  40703  a  -  275  a2  +  0.2094  a3  +  0.00155  a4 

Bearing  in  mind  the  limitations  of  the  simple  1PBT  model,  comparison  of 
the  experimental  data  and  the  simulations  shows  that  the  experimental 
results  fall  much  closer  to  the  predictions  based  on  equation  (G-l)  which 
assumes  an  increase  in  toughness  with  both  crack  length  and  crack 
velocity  than  to  the  predictions  based  on  the  constant  propagation 
toughness  assumption,  although  for  longer  crack  jumps,  the  crack  jump 
distance  is  underpredicted  when  using  equation  (G-l). 

RESULTS  OF  STRAIN  RECORD  ANALYSIS 

The  normalized  3train  gage  records  for  several  1PB  experiments  are 
plotted  in  Figures  G-3  and  G-4.  Figure  G-3  shows  the  records  for  three 
crack  jump  experiments  in  which  the  crack  arrested  after  less  than  2  mm 
extension.  The  records  in  Figure  G-3(a)  and  G-3(b)  were  obtained  from 
small  specimens  (Specimens  CJ7  and  GJ5,  respectively),  whereas  the  data 
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in  Figure  G-3(c)  were  obtained  from  a  large  specimen  (Specimen  HI). 

Figure  G-4  shows  the  stress  intensity  records  for  1PB  experiments  in 
which  the  specimens  were  completely  broken.  Figure  G-4(a)  shows  the 
record  for  a  small  specimen  (C6)  and  Figures  G-4(b)  shows  the  record  for 
a  large  specimen  (H7)  used  to  determine  the  dynamic  initiation  toughness 
data  in  Table  G-l.  Figure  G-3  also  plots  the  strain  history  that  would 
have  resulted  had  the  crack  not  propagated.  The  horizontal  dashed  line 
represents  the  strain  level  at  which  the  crack  initiates  (corresponding 
to  a  stress  intensity  of  57  MPa  m^^). 

Simulations  of  the  strain  records  for  the  experiments  illustrated  in 
Figures  G-3  and  G-4  are  presented  in  Figures  G-5  and  G-6,  respectively. 
These  simulations  were  performed  with  the  simple  1PBT  model  assuming  the 
same  two  toughness  behaviors  as  for  the  simulations  in  Figure  G-2.  For 
each  simulation,  the  strain  was  calculated  at  the  position  of  the  strain 
gage  in  the  corresponding  experiment.  On  some  simulated  strain  records, 
we  have  indicated  the  amount  of  crack  extension  associated  with  given 
amounts  of  9train  drop.  This  value  can  be  used  to  estimate  the  crack 
velocity  from  the  experimental  strain  records. 

In  the  test  from  which  the  strain  history  in  Figure  G-3(a)  was 
obtained,  the  crack  extended  by  0.5  mm,  which  is  probably  the  smallest 
crack  growth  Increment  that  can  be  detected  based  on  the  strain 
records.  It  causes  a  very  mild  flattening  of  the  top  part  of  the  strain 
history  curve.  Further,  the  strain  level  continues  to  rise  significantly 
above  the  crack  initiation  level.  This  rise  is  also  observed  in  Figures 
G-3(b)  and  G-3(c) ,  where  the  instant  of  crack  arrest  can  be  established 
as  a  point  of  inflection  in  the  decreasing  part  of  the  strain  records. 
Comparison  of  the  time  of  arrest  and  the  time  at  which  the  dashed  curves 
are  maximum  suggests  that  arrest  occurred  when  the  maximum  in  applied 
stress  intensity  was  reached. 

Comparison  of  the  experimental  records  and  the  strain  simulations 
based  on  Equation  (G-l)  shows  good  agreement  in  the  general  shape  of  the 
strain  histories  (Figures  G-3  and  G-5  and  Figures  G-4  and  G-6).  In 
particular,  the  duration  between  the  time  of  crack  initiation  (i.e.,  time 
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FIGURE  G-5  SIMULATION  OF  NORMALIZED  STRAIN  HISTORIES  FOR  CJ 
EXPERIMENTS  SHOWN  IN  FIGURE  G-3 
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FIGURE  G-6  SIMULATION  OF  NORMALIZED  ST  RAIN  HISTORIES  FOR  FRACTURE 
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when  Kj  -  57  MPa  m1^)  and  the  time  of  crack  arrest  or  the  time  at  which 
the  strain  becomes  zero  Is  well  predicted.  The  only  exception  Is  Figure 
G-6(b),  which  underestimates  the  time  required  for  the  strain  to  drop  to 
zero  in  Figure  G-4(b). 

The  good  qualitative  agreement  between  experimental  data  and 
simulations  shows  that  the  Interpretation  of  the  strain  records  is 
correct,  particularly  in  Identifying  the  inflection  In  the  unloading 
portion  as  the  point  of  crack  arrest.  Moreover,  the  fact  that 
simulations  using  a  constant  propagation  toughness  equal  to  the 
initiation  toughness  do  not  represent  the  experimental  strain  records 
suggest  that  that  assumption  Is  Incorrect.  Rather  an  increasing 
propagation  toughness  must  be  assumed. 

We  can  estimate  the  average  crack  velocity  during  the  first 
millimeters  of  crack  extension  by  comparing  the  results  of  simulation  and 
experiments.  Indeed,  as  discussed  in  Appendix  C  and  illustrated  in 
Figure  G-6,  for  a  given  strain  measurement  location,  the  distance  that 
the  crack  must  extend  to  relax  the  strain  to  zero  is  essentially 
constant.  For  a  given  experiment,  this  distance  can  be  estimated  from 
the  strain  simulations.  Then,  by  measuring  the  time  from  initiation  to 
the  point  where  the  experimental  strain  drops  to  zero,  we  can  estimate 
the  average  crack  velocity. 

The  average  crack  velocities  estimated  for  several  tests  are  given 

in  Table  G-2,  together  with  average  velocities  for  the  experiments  in 

Figure  G~3(b)  and  G-3(c).  The  latter  data  were  obtained  by  dividing  the 

crack  jump  distance  measured  after  the  test  by  the  time  between 

initiation  and  arrest  measured  on  the  experimental  records.  Table  G-2 

also  gives  the  loading  rate  to  crack  initiation  and  for  each 

Ihax 

test.  The  results  in  Table  G-2  show  that  the  initial  average  velocities 
in  all  the  tests  are  low,  not  exceeding  250  m  s-^.  A  more  detailed  look 
at  the  strain  history  records  and  simulations  indicates  that  crack 
extension  is  even  slower  during  the  first  0.5  mm  cr  so,  and  then 
accelerates  to  velocity  about  double  the  average  velocity. 


Table  G-2 


AVERAGE  CRACK  VELOCITIES  FOR  FIRST  MILLIMETERS  OF  CRACK  EXTENSION 
Specimen  K^(MPa  s  a  (m  8  ^  a^(m  8  *) 


HI 

106 

80 

34 

38 

H7 

1.9  x  106 

220 

130 

190 

CJ5 

2.5  x  106 

81 

76 

109 

CJ6 

2.85  a  106 

99 

192 

253 

CJ1 

5.7  x  106 

280 

247 

400 

• 

Estimates  of  the  crack  velocity,  aj 

extension  and  up  to  the  point  where  the 

given  in  the  last  column  of  Table  G-2. 

,  after  the  first  0.5 

strain  has  dropped  to 

The  very  low  average 

mm  of 

zero  are 

crack 

velocities  for  testa  HI  and  CJ5  indicate  that,  in  these  experiments,  the 
crack  was  extending  under  dynamically  stable  conditions;  that  is,  the 
driving  force  applied  to  the  crack  tip  by  the  bending  specimen  was  always 
equal  to  the  material  resistance  and  had  to  continuously  increase  to 
maintain  crack  extension.  Once  the  driving  force  reached  its  raaximum; 
the  crack  arrested. 

RESULTS  OF  DIRECT  KID  AND  a  MEASUREMENTS 

One  specimen  (Specimen  H7)  instrumented  with  four  strain  gages 
placed  along  the  crack  path  was  tested  in  1PB.  The  strain  records  are 
shown  in  Figure  G-7 ,  where  the  strain  history  from  four  gages  have  been 


referred  to  the  same  initial  time  and  superimposed  on  the  same  time 
scale.  Strains  in  the  vicinity  of  the  running  crack  tip  are  much  higher 
than  the  strain  near  the  initial  crack  tip,  although  the  gages  were 
placed  further  away.  Residual  strain  is  a) so  significant  for  the  gages 
placed  along  the  crack  path.  Note  that,  according  to  the  elastodynamic 
solution  for  an  extending  crack,  the  maximum  in  the  strain  history  does 
not  occur  when  the  running  crack  tip  is  immediately  underneath  the  gage; 
rather,  the  strain  is  maximum  while  the  crack  is  still  approaching  the 
gage  and  the  line  connecting  the  crack  tip  to  the  center  of  the  gage 
makes  an  angle  of  70  degrees  with  the  crack  plane.  This  distinction  is 
important  for  correct  evaluation  of  the  crack  velocity  immediately  after 
Initiation^ 

Figure  G-7  also  shows  the  results  of  a  simulation  of  the  experiment, 
using  the  propagation  toughness  behavior  given  by  equation  (G-l). 

Although  the  position  of  the  simulated  peaks  in  strain  do  not  match 
perfectly  with  the  experimental  data  and  their  amplitude  is  smaller,  the 
agreement  between  simulation  and  experiment  is  qualitatively  good.  A 
simulation  with  a  constant  propagation  toughness  equal  to  the  initiation 
toughness  is  unable  to  reproduce  any  of  the  experimental  strain 
features.  This  again  suggests  that  the  propagation  toughness  increases 
rapidly  over  the  first  millimeters  of  cracx  extension. 

From  the  strain  gage  data,  we  determined  the  average  crack  velocity 
between  two  strain  gages  and  the  dynamic  stress  intensity  factor  at  the 
propagating  crack  tip.  The  velocity  results  are  plotted  in  Figure  G-8  as 
a  function  of  the  crack  extension  Aa  .  The  symbols  in  the  figure 

represent  the  average  crack  velocities,  and  the  horizontal  lines 

represent  the  interval  over  which  they  were  measured.  The  strain  gage 

measurements  yield  an  average  velocity  over  the  first  3  to  4  mm  of  120  m 

s-*,  which  is  in  very  good  agreement  with  the  estimate  of  the  crack 
velocity  basea  on  the  drop-off  of  the  strain  near  the  initial  crack  tip 
(see  Table  G-2,  Specimen  H7).  This  result  demonstrates  that  the  latter 
procedure  is  legitimate. 
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The  serrated  strain  history  during  unloading  suggests  that  crack 
extension  occurred  by  a  series  of  small  jumps  and  arrests.  Average  crack 
speeds  can  be  calculated  from  the  arrested  crack  length,  measured  on  the 
fracture  surface,  and  from  the  time  from  initiation  to  final  crack 
arrest,  measured  on  the  strain  history  records.  The  average  velocities 
are  between  60  and  110  ms1.  Because  of  the  reloading  time  between 
incremental  jumps,  the  actual  crack  velocities  are  higher  than  the 
average  velocities.  Corrected  crack  velocity  values  are  between  110  and 
130  ms1. 

Figure  G-12(a)  shows  the  normalized  strain  history  record  for  a 

5 

specimen  loaded  at  an  intermediate  stress  intensity  rate  of  5.5  x  10 

MPa  s’1.  Kappl  in  these  tests  was  approximately  140  MPa  m1^2  and 

IMax 

the  specimens  fractured  completely  during  the  tests.  Figure  G-12(b) 
shows  strain  simulations  for  the  case  corresponding  to  Figure  G-12(a)  and 
for  two  toughness  assumptions.  In  one  case  the  propagation  toughness  is 
assumed  constant  equal  to  the  initiation  toughness;  in  the  other  case  it 
varies  according  to  equation  (G-l). 

The  serrated  appearance  of  the  unloading  portion  of  the  strain 
history  again  suggests  that  the  crack  extended  by  a  series  of  jumps  and 
arrests.  Comparison  of  the  experimental  records  with  the  strain 
simulation  shows  that  the  simulation  with  a  constant  toughness  assumption 
predicts  a  much  more  rapid  strain  drop  than  is  observed  in  reality.  A 
simulation  using  equation  (G-l)  predicts  a  strain  history  that  resembles 
the  experimental  one  more  closely,  although  the  jumps  and  arrests  are  not 
reproduced.  In  particular,  the  simulation  predicts  an  average  crack 
velocity  during  the  unloading  portion  of  the  strain  history  of  about  120 
m  s"1.  This  compares  well  with  the  velocities  calculated  from  the 
experimental  records  assuming  (based  on  the  indication  of  the  simulation) 
that  the  crack  has  extended  by  approximately  9  mm  when  the  strain  has 
dropped  to  zero.  Because  of  the  run-arrest  nature  of  the  crack 
extension,  these  estimates  represent  lower  bounds  for  the  crack 
velocity.  Corrected  crack  velocities  are  about  170  m  s-1 . 
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FIGURE  G-12  NORMALIZED  STRAIN  HISTORY  FOR  1 
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Figure  G-13(a)  shows  a  normalized  strain  history  record  for  an 
experiment  at  high  loading  rate  (3  x  10^  MPa  s-^).  This  strain 

history  is  much  smoother  than  the  records  In  Figures  G-10  and  G-12, 
suggesting  that  crack  extension  was  continuous,  Also,  the  strain  does 
not  drop  significantly  during  the  first  10  ps  after  initiation. 
Simulations  of  the  high  loading  rate  experiments  in  Figure  G-i3(b) 
indicate  that,  for  a  constant  propagation  toughness  equal  to  the 
initiation  toughness,  the  strain  should  drop  much  more  rapidly.  A 
simulation  using  equation  (G-l)  predicts  a  strain  history  resembling  more 
closely  the  experimental  one.  From  the  comparison  of  the  simulation  and 
the  experimental  records,  we  estimate  that  the  crack  extended  first  at  a 
very  slow  speed  of  about  100  ra  s  1  and  then  rapidly  accelerated  to  about 
500  m  s"1. 

The  observation  of  crack  tunneling  and  the  small  but  not  negligible 
shear  lips  on  the  fracture  surfaces  raises  the  possibility  that  the  crack 
front  geometry  and  the  shear  lips  may  affect  the  strain  history  and  the 
overall  crack  extension  behavior.  To  investigate  this  possibility,  we 
performed  experiments  on  smooth  and  side-grooved  specimens  of  Identical 
shape  and  under  the  same  loading  conditions.  The  specimens  were 
instrumented  with  a  strain  gage  at  the  fatigue  precrack  tip  and  with  an 
additional  strain  gage  positioned  9  mm  away  from  the  crack  tip  along  the 
prospective  crack  path.  The  second  strain  gage  determined  the  average 
crack  velocity  immediately  after  crack  initiation.  A  static  calibration 
of  the  strain  gages  at  the  fatigue  crack  tip  was  performed  to  relate  the 
strain  history  to  the  stress  intensity  history.  For  the  side-grooved 
specimens,  a  thickness  correction  factor  was  introduced  to  account  for 
the  effect  of  the  side  grooves. 

Figure  G-14  shows  strain  histories  near  the  fatigue  crack  tip  for  a 
smooth  and  for  a  side-grooved  specimen.  After  crack  initiation,  the 
curves  show  a  definite  difference  in  shape.  For  the  smooth  specimen,  the 
strain  first  decreases  very  gradually  before  it  rapidly  drops  off.  For 
the  side-grooved  specimen,  the  strain  drops  off  immediately ,  suggesting  a 
higher  Initial  crack  velocity  when  shear  lips  are  suppressed. 
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FIGURE  G-13  NORMALIZED  STRAIN  HISTORY  FOR  1 2.7-mm-TH ICK  SPECIMEN 
LOADED  AT  HIGH  RATE 


Measurements  of  the  average  crack  velocity  with  strain  gages  supports 
this  interpretation  with  an  average  crack  velocity  of  360  m  in  the 
side-grooved  specimen  and  270  m  s-^  in  the  the  smooth  specimen.  However, 
the  rate  of  the  rapid  dropoff  in  the  smooth  and  side-grooved  specimens  is 
the  same,  indicating  similar  crack  velocities  during  this  part  of  the 
crack  extension.  Therefore,  the  measured  difference  in  the  average  crack 
velocity  appears  essentially  due  to  the  very  slow  initial  crack  velocity 
in  the  smooth  specimen.  In  turn,  this  indicates  that  the  resistance  to 
crack  propagation  increases  significantly  above  the  initiation  toughness 
level  as  the  crack  extends.  Strain  history  simulations  for  the 
experiments  of  Figure  G-13  show  that  a  significant  increase  in 
propagation  toughness  above  the  initiation  toughness  must  be  postulated 
to  duplicate  the  experimental  results  for  the  smooth  specimen. 

In  summary,  the  experiments  with  12.7-mm  material  have  provided  the 

following  results.  The  crack  velocity  in  the  1PB  tests  increases  with 

the  loading  rate.  At  lower  loading  rates  <3.2  x  105  and  5.5  x  106  MPa 

m  ' "  s  A),  crack  extension  is  discontinuous  and  proceeds  by  a  series  of 

small  crack  jumps  and  arrests.  Estimates  of  the  average  crack  speed 

range  from  130  m  s-1  at  a  loading  rate  of  3.2  x  105  MPa  s'-1  to  170  m 

at  a  loading  rate  of  5.5  x  10^  M?e.  s-^  to  500  m  s~*  for  the 

highest  loading  rate  of  3  x  106  MPa  a-1.  Simulations  of  the  strain 

histories  indicate  that  the  propagation  toughness  must  be  increased 
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significantly  (20  MPa  tn  'or  more)  above  the  initiation  level  during  the 
first  millimeters  of  crack  growth  to  reproduce  the  observed  crack 
speeds.  Finally,  the  experiments  with  smooth  and  side-grooved  specimens 
demonstrate  that  the  formation  of  even  small  shear  lips  markedly  aifects 
the  crack  propagation  behavior. 

DISCUSSION 

All  the  experimental  results  discussed  above  together  with  the 
strain  history  simulations  show  that,  immediately  after  initiation,  the 
crack  extension  resistance  increases  above  the  crack  initiation  toughness 
level.  The  crack  jump  experiments  indicate  that  an  increase  in  toughness 
of  over  20  MPa  m^2  occurs  during  the  first:  2  mm  of  crack  extension. 


Direct  measurements  with  strain  gages  yield  a  propagation  toughness  of 
140  MPa  at  a  crack,  speed  of  about  400  m  s-^  (Figure  G-9),  which  is 

over  twice  the  value  of  the  initiation  toughness.  These  observed 
increases  in  propagation  toughness  may  be  a  function  of  at  least  two 
parameters:  the  amount  of  crack  extension  Aa  and  the  crack  velocity. 

Most  of  the  dynamic  crack  propagation  work  has  focused  on  the 

2  3 

relationship  between  crack  velocity  and  propagation  toughness  KjD«  ’ 
Although  it  is  generally  admitted  that  K-^  13  a  unique,  increasing 
function  of  the  crack  velocity,  it  is  unlikely  that  the  increase  in 
toughness  observed  in  the  crack  .lump  experiments  is  due  to  crack  velocity 
effects.  First,  the  veloctles  measured  in  the  crack  jump  experiments  are 
low,  whereas  significant  increases  in  for  high  strength  steels  are 

generally  observed  above  500  m  s~ ^  only.  Second,  specimens  having 
different  crack  velocities  displayed  the  same  behavior  In  the  crack  jump 
experiments  (specimens  CJ5  and  HI  in  Table  G-2). 

Therefore,  we  conclude  that  the  Initial  increase  in  the  propagation 
toughness  above  the  initiation  value  in  the  9.8-mm-thick  material  is  a 
function  of  the  amount  of  crack  extension  Aa.  In  other  words,  there  is  a 
resistance  curve  effect,  in  the  dynamic  case  as  is  often  observed  in 
static  fracture  experiments.  The  crack  jump  experiments  and  their 
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simulations  indicate  that  Xjq  increases  from  the  value  of  57  MPa  m  at 
initiation  to  81  MPa  after  2  mm  of  extension.  They  further  suggest 

that,  beyond  this  steady  state  value,  Kj-p  does  not  depend  significantly 
on  the  amount  of  crack  extension.  We  presume  that,  as  the  crack  extends 
and  the  crack  speed  increases,  the  dependence,  on  the  velocity  becomes 
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more  important  and  contributes  to  raise  the  propagation  toughness 
further.  Evidence  for  this  behavior  is  provided  by  the  direct 
measurements  of  the  propagation  toughness  with  9train  gages  9hown  in 
Figure  G-9  which  indicate  values  of  140  MPa  nr  ^. 

Therefore,  the  experiments  with  the  9.8-mm-thick  specimens  indicate 
that  the  propagation  toughness  depends  both  on  the  amount  of  crack 
extension  and  on  the  crack  velocity.  Equation  (G-l),  which  was  used  for 
the  strain  simulations,  is  a  crude  attempt  to  represent  and  quantify  thiB 
dependence.  It  was  chosen,  based  on  the  experimental  results,  to 
represent  all  the  data  with  a  single  toughness  behavior.  In  particular, 
the  dependence  on  the  amount  of  crack  extension  was  determined  by  fitting 
linear  segments  to  the  versus  Aa  curve  of  Figure  G-2.  Due  to  the 

limitations  of  the  1FBT  model,  complete  agreement  with  all  the 
experimental  results  cannot  be  expected,  particularly  for  the  larger 
crack  propagation  distances.  However,  we  believe  that  the  simulations 
are  adequate  for  predicting  the  early  stage  of  crack  extension  and  for 
qualitatively  assessing  the  overall  crack  propagation  behavior. 

The  crack  propagation  results  for  the  9.8-mm-thick  specimens  are  in 

contrast  with  several  published  results  for  very  similar  materials,  which 

indicate  that  the  propagation  toughness  i9  a  unique  monotonic  function  of 

the  crack  velocity,  which  is  essentially  equal  to  the  initiation 

—1  7  8 

toughness  up  to  crack  velocities  of  500  m  s  or  more.  *  However,  other 
publications  also  report  values  of  the  propagation  toughness  in  4340 
steel  that  are  significantly  larger  than  the  initiation  toughness  at 
crack  velocities  as  low  as  80  m 

The  results  of  the  experiments  on  12. 7-mm-thick,  smooth  and  side- 
grooved  specimens  provide  the  basis  for  a  tentative  explanation  of  the 

^G.  T.  Hahn,  R.  G.  Hoagland,  A.  R.  Rcsenfield,  and  R.  Sejnoha,  "Rapid 
Crack  Propagation  in  a  High  Strength  Steel,”  Met.  Trans.  Vol.  5,  pp.  475- 
482  (1974). 

^S.  J.  Hudak,  R«  J.  Dexter,  J.  H.  Fitzgerald,  and  M.  F.  Kanninen,  "The 
Influence  of  Specimen  Boundary  Conditions  on  the  Fracture  Toughness  of 
Running  Cracks,"  Eng.  Fract.  Mech.,  Vol.  23,  pp.  201-213  (1986). 


Therefore,  our  investigation  shows  that  nart  of  the  difference 
between  the  dynamic  initiation  and  propagation  toughness  is  due  to  a 
geometrical  and  a  structural  effect  associated  with  changes  in  constraint 
along  a  propagating  crack  front  and  to  the  formation  of  small  shear 
lips.  This  explanation  is  consistent  with  the  observations  of  Hahn,  et. 
al.4 

Practical  conclusions  can  be  drawn  from  the  results  reported  here. 
First,  unless  the  thickness  of  the  fractured  specimen  or  part  is  large, 
there  will  be  a  resistance  curve  effect,  even  in  the  dynamic  case,  due  to 
changes  in  the  crack  front  geometry  during  the  early  part  of  crack 
extension.  Because  of  this  resistance  curve  effect,  a  rising  applied 
stress  intensity  is  required  to  drive  the  crack  at  a  significant  speed. 

In  other  words,  there  is  the  possibilty  for  dynamic  stable  crack 
extension  over  very  short  distances.  Second,  our  results  show  that  the 
thicknesses  sufficient  to  obtain  valid  dynamic  initiation  toughness 
values  may  not  be  sufficient  for  obtaining  thickness-independent 
propagation  toughness  data.  Third,  because  of  the  three-dimensional 
nature  of  the  crack  propagation  process,  caution  should  be  applied  in 
interpreting  experimental  data  obtained  from  surface  measurements. 
Finally,  the  results  presented  here  have  demonstrated  that  the  strain 
history  record  obtained  from  a  strain  gage  placed  near  the  initial  crack 
tip  can  be  succesfully  analyzed  to  obtain  crack  velocity  estimates  and  at 
least  a  qualitative  assessment  of  the  propagation  toughness. 


Appendix  H 


INVESTIGATION  OF  DYNAMIC  CRACK  EXTENSION  BY  MEANS  OF  HEAT 
OF  FRACTURE  MEASUREMENTS 

This  appendix  reports  on  a  series  of  crack  extension  experiments  in 
which  the  crack  was  initiated  from  a  blunt  notch  by  quasi-static 
loading.  The  first  objective  of  these  experiments  was  to  measure  the 
dynamic  propagation  toughness  for  approximately  steady  crack  velocity 
conditions  (in  contrast  to  the  work  reported  in  Appendix  G,  where  we  were 
interested  in  the  transient  regime  immediately  following  initiation). 

The  second  objective  of  the  experiments  was  to  compare,  by  means  of 
independent  measurements,  the  energy  delivered  to  the  crack  tip  and  the 
energy  absorbed  during  crack  propagation.  The  third  objective  was  to 
evaluate  several  experimental  techniques  for  measuring  the  dynamic 
propagation  toughness:  thermal  measurement  techniques  with  thermistors  or 
thermocouples  and  a  strain  gage  technique. 

For  reasons  that  will  be  discussed  below,  a  direct  comparison 
between  energy  delivered  to  the  crack  tip  and  energy  absorbed  could  not 
be  performed  from  the  experiments.  However,  the  experiments  yielded 
propagation  toughness  values  consistent  with  the  values  reported  in 
Appendix  G  and  highlighted  several  advantages  and  limitations  of  the 
experimental  techniques  we  used. 

PRINCIPLES  OF  THE  EXPERIMENTS 

During  the  experiments  described  in  this  Appendix,  we  intended  to 
measure  the  applied  dynamic  stress  intensity,  Kj^^n(t),  by  means  of 
strain  gages  and  to  derive  the  dynamic  propagation  toughness,  KjD(a), 
from  temperature  measurements  obtained  using  thermocouples  and 
thermistors.  The  principles  of  these  two  approaches  are  summarized 
below. 


■A'Jl/Jf/VVY,rfVV  f  X.j!i  vV  .  i  j  ri  Vj  vi  ^\r«-  -j»\.  rv*  n.%  - 


where  H(P ,  k,  0^)  =  /sin  &max 


[♦  (p.lt.e  )  -  v*  (P.k.O  )], 

y  max  x  max 


0max  can  be  obtained  by  differentiating  the  expression  for  Gy  with 
respect  to  0  and  is  a  function  of  crack  velocity  and  Poisson's  ratio.  As 
tabulated  in  Table  H-l,  the  values  of  0max  over  the  range  of  interest  for 
a/C  and  v  (d/C,.  =  0  to  0.6  and  v  =  0.1  to  0.49)  vary  only  between  69° 

t  G 

and  73°.  The  function  H(P,k,0max) 
is  plotted  in  Figure  H-2. 

Thus  it  is  seen  from  equation  (H-4),  that,  if  the  crack  velocity  a 
and  the  strain  e^,  at  a  small  distance  y  from  the  crack  plane  are  known, 
say,  from  strain  gage  measurements,  the  stress  intensity  factor  can 
easily  be  determined.  The  crack  velocity  can  be  estimated  by  measuring 
the  time  Interval  between  the  peak  strains  recorded  by  two  neighboring 
gages  placed  at  a  known  distance  from  each  other  along  the  crack  path. 

By  mounting  a  series  of  strain  gages  parallel  to  the  specimen  center 
line,  we  can  measure  a  Kj^yn-d  history  at  regular  intervals  as  the  crack 
propagates.  This  approach  was  taken  in  the  present  experiments. 


The  strain  gage  method  of  measuring  K^dyn(t)  is  based  on  two 
limiting  assumptions.  First,  it  uses  the  analytical  solution  for  a  crack 
extending  at  a  constant  velocity  in  an  infinite  plate,  and  therefore,  it 
does  not  account  for  the  effect  on  the  crack  tip  field  of  stress  waves 
reflected  at  the  boundaries.  For  a  finite  geometry,  it  has  been  shown 
that  stress  waves  emanating  from  the  extending  crack  tip,  reflecting  at 
the  boundaries,  and  returning  to  the  crack  tip  can  cause  significant 
changes  in  the  instantaneous  stress  intensity  factor  and  result  in  a 
transient  change  in  the  crack  velocity.  The  second  limitation  of  the 
strain  gage  method  is  that  it  considers  only  the  singular  term  of  the 
series  expansion  for  the  crack  tip  fields.  However,  in  finite  size 
specimens,  higher  order  terms  may  also  be  important.^  An  additional 
limitation  is  the  restricted  number  of  measurement  positions,  which 
obviously  limits  the  ability  to  measure  rapid  transients  in  crack 
velocity  and  stress  intensity  factor.  Therefore,  this  technique  must  be 
used  with  caution. 
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The  dynamic  fracture  toughness  of  the  propagating  cracks  was  derived 
from  the  temperatures  measured  on  the  specimen  surface  near  the  running 
crack.  Nearly  all  the  energy  absorbed  by  a  propagating  crack  is  absorbed 
by  plastic  flow  of  material  in  a  relatively  small  (compared  with  specimen 
dimensions)  volume  in  the  vicinity  of  the  crack  tip.  Therefore,  the 
fracture  toughness  associated  with  the  propagating  crack  is  determined  by 
the  plastic  energy  expended  at  the  crack  tip.  Taylor  and  Quinney^  showed 
that  for  metals  nearly  all  the  plastic  energy  is  dissipated  as  heat. 

Heat  generated  at  a  running  crack  tip  flows  away  from  the  crack  plane 
into  both  halves  of  the  specimen,  causing  small  transient  temperature 
elevations.  If  the  temperature  history  of  material  near  a  running  crack 
can  be  measured,  the  heat  energy  and  hence  the  plastic  energy  and  the 
propagation  toughness  can  be  calculated  by  solving  the  appropriate  heat 
flow  equation. 


Thermocouples  and  thermistors  were  used  in  these  experiments  to 
measure  near-crack-tip  temperature  histories.  This  method  was  used  first 
by  Wells**  and  later  improved  and  applied  by  several  other 
investigators. ^-1° 


To  obtain  the  heat  generated  by  the  propagating  crack  Q  and  the 

dynamic  fracture  toughness  KjD,  we  used  the  measured  maximum  temperature 

values  ATmax  and  the  solution  of  the  one-diraensional  heat  flow  equation 

11-13 

for  an  instantaneous  heat  source. 

KTD 

Q  =  AT  y  /Hie'  pc  =  (H-5) 

max  E 


where  y  is  the  distance  between  the  fracture  surface  and  the  temperature 
tranducer,  e  =  2.718,  and  p,  c,  and  E  are  the  mass  density,  the  specific 
heat,  and  Young’s  modulus,  respectively. 

The  adequacy  of  the  models  used  for  interpreting  the  experiments  and 
the  measurements  error  for  the  strain  and  temperature  data  are  discussed 
along  with  the  results. 


EXPERIMENTAL  PROCEDURES 


Six  transversally  wedge-loaded  compact  tension  (CT)  specimens  were 
prepared  from  a  commercial  aircraft  quality  4340  steel  plate.  The  rough 
machined  specimens  were  heat-treated  to  HRC  50  and  then  ground  to  the 
final  dimensions  and  tolerances  shown  in  Figure  H-3.  The  loading  hole 
and  the  starter  notch  were  put  In  by  electrodischarge  machining. 

Chevrons  of  varying  acuity  were  introduced  at  the  bottom  of  the  notch 
using  a  thin  grinding  wheel. 

Each  specimen  was  instrumented  with  four  miniature  strain  gages 
evenly  spread  along  the  anticipated  crack  path  and  about  1/8  inch  from 
the  specimen  center  line  (Figure  H-3).  These  strain  gages  measured  the 
crack  velocity  as  well  as  the  applied  dynamic  stress  intensity  factor 
K^yn  at  the  running  crack  tip.  in  addition  to  the  strain  gages,  two 
thermocouples  and  four  thermistors  were  mounted  on  each  specimen.  The 
active  junction  of  the  thermocouples  consisted  of  a  0.2-mm  constantan 
wire  spot-welded  to  the  specimen  a  few  millimeters  away  from  the  crack 
path.  Another  contantan  wire,  spot  welded  in  a  corner  of  the  specimen 
away  from  the  crack,  acted  as  reference  junction.  The  thermomistors 
(Siemens  K19)  were  less  than  0.4  mm  in  diameter  and  were  held  on  the 
specimen  with  heat  conducting  paste. 

The  precise  positions  of  the  temperature  sensors  on  the  specimen 
varied  from  test  to  test,  but  they  were  more  or  less  evenly  spaced  along 
the  crack  path  at  a  distance  not  exceeding  10  mm.  In  some  cases  one 
thermocouple  and  one  thermistor  were  mounted  as  close  to  each  other  as 
possible  so  that  they  would  see  essentially  the  same  temperature 
history.  This  arrangement  made  a  direct  comparison  of  the  two  types  of 
transducers  possible. 

During  loading  of  the  specimen,  the  crack  mouth  opening  displacement 
(CMOD)  was  measured  using  a  clip  gage. 

Figure  H-4  diagrams  the  electronic  circuitry  that  was  used  to 
process  and  record  the  signals  from  the  various  transducers.  The  strain 
gage  signals  were  recorded  on  dual-channel  digital  oscilloscopes  with  a 
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1-ms  time  window.  The  temperature  histories  were  recorded  on  a  Honeywell 
eight-beam  strip  chart  recorder.  During  loading,  the  CMOD  was  also 
recorded  on  the  strip  chart  recorder. 

To  trigger  the  oscilloscopes,  we  used  a  laser  beam  directed  at  a 
spot  a  few  millimeters  ahead  of  the  specimen  starter  notch  and  reflected 
by  the  polished  specimen  surface  onto  a  photodiode.  When  the  crack 
becomes  unstable,  it  displaces  a  small  plastic  zone  around  its  moving 
tip.  As  the  plaBtic  zone  crosses  the  laser  beam,  it  deflects  it  away 
from  the  photodiode  and  the  voltage  change  at  the  photodiode  output 
triggers  the  oscilloscopes. 

After  checking  that  the  response  of  the  thermistors  was  linear  in 
the  temperature  range  of  interest,  we  calibrated  the  thermistors  against 
a  thermocouple  of  known  temperature  sensitivity.  A  similar  procedure  was 
used  to  calibrate  the  thermocouples. 

The  specimens  were  loaded  to  fracture  in  a  10-kips  screw-driven 
Instron  universal  testing  machine.  The  loading  arrangement  is  shown  in 
Figure  H-5.  Two  support  steel  plates  were  inserted  between  the  two 
halves  of  the  split-pin  and  the  base  plate  to  prevent  any  axial  loading 
of  the  specimen  and  to  suppress  friction  at  the  specimen  support  plate 
interface.  For  all  the  experiments,  the  machine  cross-head  speed  was  8 
mm  s~*.  At  this  cross-head  rate,  loading  to  fracture  took  less  than  2 
seconds • 

The  testing  machine  was  manually  operated.  Due  to  the  slow  reaction 
time  of  the  operator  and  the  inertia  in  the  loading  train,  it  took  about 
half  a  second  for  the  cross-head  to  stop  after  the  instability  and  crack 
arrest  events.  As  discussed  below  this  delay  in  interrupting  the 
specimen  loading  after  the  test  was  completed  is  significant  for 
interpretation  of  the  test  results. 

EXPERIMENTAL  RESULTS 

Two  of  the  six  experiments  that  were  performed  gave  a  complete  set 
of  strain  and  temperature  measurements  (specimens  CT3  and  CT6).  For  one 


specimen  (Cl’l),  the  strain  gage  signals  could  not  be  recorded  because  the 
sensitivity  of  the  oscilloscopes  was  improperly  set.  For  two  other 
specimens  (CT4  and  CT5),  the  strain  gage  data  were  lost  because  of 
trigger  problems.  Finally,  for  specimen  CT2,  there  was  a  first 
initiation  at  a  relatively  low  KIQ;  then  the  crack  arrested,  reinitiated 
upon  further  loading,  and  arrested  again  after  a  total  extension  of  only 
24  mm.  The  multiple  initiations  rendered  the  interpretation  of  the 
temperature  data  impossible.  Small  shear  lips  formed  in  all  six 
specimens  tested.  The  width  and  height  of  the  shear  lip  were  about  1.5 
mm  near  the  initial  notch  and  tapered  to  about  1  mm  in  the  arrest  region. 

The  Kjq's  and  arrest  crack  lengths  for  all  six  specimens  are  given 
in  Table  H-2.  Tests  CT3  through  CT6  have  similar  K^q  values  and  arrest 
crack  lengths.  Based  on  this  observation  and  on  the  work  of  Hoagland  et 
al.,^  we  concluded  that  all  four  specimens  had  about  the  same  crack 
propagation  history.  Thus  the  propagation  toughness  data  obtained  in 
tests  CT4  and  CT5  can  be  evaluated  using  the  strain  histories  and  crack 
velocities  measured  in  tests  CT3  and  CT6. 

Figure  H-6  plots  the  strain  histories  obtained  with  the  four  strain 
gages  in  tests  CT3  and  CT6.  In  each  plot,  the  strain  gage  signals  have 
been  superimposed  in  such  a  way  that  they  have  the  same  strain  baseline 
and  the  same  time  base.  The  strain  histories  for  the  two  experiments  are 
very  similar.  The  observable  differences  between  the  strain  histories  in 
tests  CT3  and  CT6  occur  because  the  row  of  evenly  spaced  strain  gages  was 
positioned  further  away  from  the  original  notch  in  test  CT6  than  in  test 
CT3. 

The  sharp  strain  peaks  are  caused  by  the  passage  of  the  crack  tip 
near  the  strain  gage.  The  signals  from  gage  4  in  both  test  have  a 
different  shape  than  the  signals  from  gages  1  to  3  because  in  both  tests 
strain  gage  4  was  located  very  close  to  the  arrested  crack  tip.  Also, 
the  strain  histories  for  gage  4  do  not  have  a  sharp  peak,  the  strain  does 
not  decrease  to  zero  after  going  through  the  maximum,  and  there  is  a 
significant  oscillation  in  the  late  time  strain  history.  This 
oscillation  indicates  that,  after  crack  arrest,  the  material  In  the 
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region  of  the  arrested  crack  tip  continues  to  be  loaded  by  the  vibrating 
specimen. 

The  propagation  time  from  one  gage  location  to  the  other  can  be 
directly  measured  in  Figure  H-6  as  the  time  difference  between  successive 
peaks.  Using  the  singularity  solution  to  estimate  the  position  of  the 
crack  tip  with  respect  to  the  gage  position  at  the  time  of  the  strain 
maximum,  we  can  estimate  the  crack  velocity.  The  strain  signals  before 
the  sharp  peaks  give  a  measure  of  the  strain  in  the  specimen  during 
quasi-static  loading  before  crack  extension,  and  the  signals  after  the 
sharp  peaks  show  the  residual  strains  after  crack  propagation.  Figure  H- 
6  shows  that  the  maximum  strains  measured  by  gages  1,  2,  and  3  were  very 
large  (in  excess  of  6000  raicrostralns)  and  that  significant  residual 
stresses  were  left  in  the  wake  of  the  crack  (in  excess  of  1000 
microstrains).  These  high  strain  values  suggest  that  the  material  at  or 
near  the  gage  location  (about  3  ram  from  the  crack  plane)  may  have 
yielded.  This  result  is  consistent  with  the  estimate  of  the  radius  of 
the  plastic  zone  of  about  2  mm  based  on  the  shear  lip  size. 

Figure  H-7  plots  the  crack  velocity  profiles  as  a  function  of  crack 
position  obtained  from  the  strain  gage  measurements  in  tests  CT3  and 
CT6.  Estimates  for  the  crack  velocity  profiles  in  tests  CT4  and  CT5  are 
also  indicated  in  the  figure.  These  estimates  were  obtained  in  the 
following  way.  Because  the  Kjq  values  and  the  arrested  crack  length  were 
similar  in  all  the  tests,  we  assumed  that  the  crack  velocity  profiles 
from  tests  CT3  and  CT6  would  also  apply  to  tests  CT4  and  CT5.  The 
assumed  initial  velocity  was  slightly  raised  or  lowered  to  account  for 
the  differences  in  K-j-q.  The  velocity  curves  of  test  CT3  (for  test  CT5) 
and  CT6  (for  test  CT4)  were  then  displaced  to  fit  through  the  points  Aa  ■ 
0,  a  -  a^nit  and  Aa  «  Aaar,  A  =  0.  The  crack  velocity  profiles  of 
Figure  H-7  were  used  in  plotting  the  stress  intensity  and  propagation 
toughness  data  measured  with  the  strain  gages,  the  thermocouples,  and  the 
thermistors • 

Figures  H-8  and  H-9  show  temperature  histories  recorded  near  the 
crack  plane  in  tests  CT3  and  CT6,  respectively.  The  dashed  vertical 
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lines  in  the  figures  indicate  the  instant  of  crack  initiation.  After 
crack  extension  and  arrest,  simple  one-dimensional  heat  flow 
considerations  show  that  the  temperature  at  the  transducer  location 
should  gradually  rise,  pass  through  a  maximum,  and  then  slowly  decrease 
as  the  specimen  reaches  thermal  equilibrium.  The  amplitude  of  the 
temperature  maximum  is  a  function  of  the  position  of  the  transducer  and 
of  the  intensity  of  the  heat  source  (equation  H-5). 

During  quasi-static  loading  of  the  specimens,  the  temperature 
transducers  register  a  temperature  change  that  is  linear  with  time 
(portion  of  the  temperature  histories  to  the  left  of  the  vertical  dashed 
lines).  This  temperature  change  is  associated  with  the  thermoelastic 
effect  in  the  specimen,  which  causes  cooling  of  the  material  at  points 
where  the  mean  stress  is  tensile  and  heating  at  points  where  the  mean 
stress  is  compressive.  If  the  loading  is  fast  enough,  this  temperature 
change  occurs  quasi-adiabatically,  and  it  is  reversible  once  the  elastic 
stress  is  relieved  by  the  passage  of  the  crack.  Thus  after  crack  arrest, 
the  temperature  is  expected  to  go  back  to  the  base-line  value  almost 
instantaneously,  before  gradually  rising  because  oii  heat  flow  from  the 
fracture  surface. 

Figures  H-8  and  H-9  indicate  that,  alt,.uugh  for  the  therrac uples 
the  recovery  of  the  thermoelastic  temperature  excursion  is  relative]” 
fast  and  complete.  Figure  H-8(a)  and  (b) ,  and  Figure  H-9(a),  this  is  not 
the  case  for  the  thermistors,  which  in  some  instances  show  considerable 
delay  and  only  partial  temperature  recov  y.  The  reasons  for  this 
behavior  are  not  entirely  clear.  It  may  be  due  to  poor  heat  conduction 
of  the  cement  used  for  mounting  the  thermistors  or  to  residual  stresses 
that  prevent  complete  material  unloading  after  the  passage  of  the  crack 
and  hence  complete  temperature  recovery. 

Another  anomalous  temperature  behavior  is  evidenced  in  Figure  H- 
9(b) ,  where  the  temperature  histories  measured  by  thermistors  TH2A  and 
TH2B  first  decrease  after  the  fracture  event,  go  through  a  minimum,  and 
then  increase  as  expected  from  heat  flow  considerations.  This  behavior 
is  characteristic  of  the  transducers  placed  near  the  arrested  crack  tip 
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region  and  ia  due  to  the  thermoelastic  effect.  After  crack  extension  and 
arrest,  the  loading  machine  was  stopped  manually  by  the  operator. 
Therefore,  a  delay  of  0.5  to  1  second  occurred  between  crack  arrest  and 
shut-off  of  the  loading  machine.  Because  high  displacement  rates  were 
used  to  drive  the  wedge  (to  minimize  heat  flow  caused  by  the 
nonequilibrium  thermoelastic  temperature  distribution  during  loading), 
this  small  delay  is  sufficient  to  permit  reloading  of  the  arrested  crack 
tip,  which  in  turn  results  in  thermoelastic  cooling  of  the  material  in 
the  crack  tip  region  that  is  under  tensile  stress.  The  thermoelastic 
cooling  is  superimposed  on  the  temperature  elevation  due  to  heat  flow 
from  the  crack  plane  and  makes  interpretation  of  the  temperature 
histories  impossible.  Therefore,  the  temperature  data  from  such  records 
have  oeen  discarded. 

The  stress  intensity  factor  K^^11  was  calculated  from  the  strain  gage 
data  using  equation  (H-3).  The  propagation  toughness  was  calculated  from 
t  '•  temperature  data  using  equation  (H-5)  and  from  the  maximum 
temperature  amplitude  ATmax,  defined  as  the  difference  between  the 
maximum  in  the  temperature  history  and  the  baseline  temperature  just 
before  the  test.  Because  the  thermoelastic  temperature  excursion  was  not 
always  recovered,  this  definition  yields  a  lower  bound  for  the 
propagation  toughness.  The  propagation  toughness  values  calculated  with 
equation  (H-5)  were  also  corrected  to  account  for  the  finite  width  of  the 
plastic  zone.  The  correction  is  necessary  because,  as  evidenced  by  the 
shear  lips  at  the  specimen  edge,  the  plastic  zone  is  about  4  mm.  Thus, 
the  distance  from  the  temperature  transducer  to  the  crack  plane  can  no 
longer  be  considered  large  compared  with  the  width  of  the  instantaneous 
heat  source.  To  correct  the  propagation  toughness  values,  we  calculated 
the  temperature  history  at  various  distances  from  the  heat  source  for  the 
case  of  an  instantaneous  line  heat  source  and  for  the  case  of  an 
instantaneous  heat  source  distributed  uniformly  over  a  width  2ry.  The 
total  heat  dissipation  per  unit  crack  area  was  assumed  the  same  in  both 
calculations.  Figure  H-10  compares  the  results  for  a  calculation 
assuming  a  heat  source  width  2ry  of  4  mm  with  the  results  for  a  line  heat 
source.  The  history  of  temperature  elevation,  normalized  by  the  heat 


source  amplitude  Q,  is  plotted  for  several  distances  from  the  center  of 
the  heat  source.  At  a  given  point  and  for  a  given  value  of  Q,  use  of  the 
line  heat  source  model  predicts  a  lower  temperature  elevation  than  the 
finite  width  heat  source.  Conversely,  use  of  a  given  temperature 
elevation  with  the  line  heat  source  model  predicts  a  higher  value  of  Q 
than  the  use  of  the  same  temperature  elevation  with  the  finite  width  heat 
source  model.  The  correction  factor  for  the  appropriate  measurement 
distance  for  each  temperature  transducer  was  calculated  from  the 
differences  in  T/Q-amplitude  in  Figure  H-10  and  applied  to  the  Q  values 
calculated  with  equation  (H-5). 

The  stress  intensity  values  calculated  from  the  strain  gage  data  and 
the  propagation  toughness  values  calculated  from  the  temperature  data  are 
plotted  in  Figure  H-ll  as  a  function  of  crack  velocity  and  in  Figure  H-12 
as  a  function  of  distance  from  the  original  crack  tip.  Also  plotted  in 
Figure  H-12  is  a  curve  giving  the  statically  calculated  stress  intensity 
factor  as  a  function  of  crack  length  for  initial  conditions 
representative  of  those  in  test  GT3  through  CT6  (crack  opening 
displacement  at  initiation  =  1.7  mm,  aQ/Wg  =>  0.402).  The  value  KID  shown 
in  Figures  H-ll  and  H-12  represents  the  dynamic  initiation  toughness 
measured  independently  with  the  one-point-bend  test  for  the  particular 
batch  of  4340  steel  used  for  the  CT  crack  extension  experiments. 

In  Figure  H-ll,  the  K^n  values  obtained  from  strain  gage 
measurements  lie  significantly  above  the  data  obtained  from  the 
temperature  measurements.  The  KjD  data  show  a  large  scatter  and  most  of 
the  data  are  grouped  in  the  velocity  range  300  to  600  m  because  of 
the  elimination  of  the  data  points  affected  by  the  post-arrest 
thermoelastic  effect.  In  Figure  H-12,  the  data  again  lie  well 

above  the  KID  data  but  also  above  the  static  curve.  This  result 
contradicts  previously  published  results,  which  show  that  during  the 
early  part  of  crack  extension  K  J  is  below  K^.  and  at  later  times 
the  opposite  situation  prevails. “  However,  the  Kjp  values  show  the 
expected  trend,  first  falling  below  the  static  curve  and  then  approaching 
or  crossing  the  static  curve  close  to  the  arrest  crack  length.  Note  that 
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the  kID  and  the  measured  In  the  CT  experiments  are  much  higher  than 

the  dynamic  initiation  toughness 

The  objective  in  using  two  types  of  temperature  sensors  for  the 
thermal  measurements  was  to  evaluate  the  relative  advantages  of 
thermocouples  and  thermistors.  Table  H-3  compares  the  results  of  pairs 
of  thermocouple-thermistors  positioned  at  the  same  location.  In  each 
case  the  general  shape  of  the  temperature  history  pairs  agree  well  (see 
Figures  H-8  and  H-9).  There  is,  however,  a  systematic  difference  in  the 
maximum  temperature  amplitude  that  results  in  a  systematic  difference  in 
the  K-j-q  values.  Thus,  the  KID  values  obtained  from  the  thermistor 
measurement  are  lower  than  the  values  obtained  from  the  thermocouple 
measurements.  The  reason  for  this  systematic  difference  has  not  been 
established.  Note  that  because  the  toughness  depends  on  the  square  root 
of  the  temperature,  the  differences  in  propagation  toughness  estimates 
are  relatively  small  compared  with  the  temperature  differences. 

DISCUSSION  AND  CONCLUSIONS 

The  propagation  toughness  measured  in  the  present  experiments  is 
much  higher  than  expected  from  published  data  for  similar  materials1*’’17 
or  from  comparisons  with  the  dynamic  initiation  toughness.  This  high 
value  of  the  propagation  toughness  is  associated  with  a  large  plastic 
zone,  which  affects  our  ability  to  measure  the  applied  stress  intensity 
factor  K^yn  using  strain  gages.  Indeed,  the  high  strains  measured 
during  the  experiments  (Figure  H-6)  and  the  large  residual  strains 
measured  after  crack  arrest  demonstrate  that  the  material  at  and  near  the 
strain  gage  locations  (3  mm  from  the  crack  plane)  was  deformed  beyond  the 
elastic  limit.  Therefore,  the  proposed  approach  of  using  the 
ela8todynamic  solution  to  obtain  the  stress  intensity  applied  at  the 
running  crack  tip  cannot  be  used,  and  the  values  calculated  with  equation 
(H-4)  are  gross  overestimates  because  they  do  not  account  for 
plasticity.  This  discrepancy  is  further  demonstrated  in  Figure  H-12 
where  the  Kj^n  values  calculated  from  the  strain  gage  measurements  are 
greater  than  the  Katat  values.  It  is  physically  impossible  for  K^n 
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to  be  greater  than  because,  during  the  early  part  of  crack 

extension,  part  of  the  energy  released  must  be  used  to  overcome  the 
inertia  of  the  specimen  and  hence  cannot  be  delivered  to  the  crack  tip. 

Although  the  strain  gage  measurements  failed  to  provide  values 

because  of  plasticity,  they  are  still  usable  to  determine  the  crack 
velocity.  The  main  limitation  on  the  crack  velocity  measurements  during 
the  CT  experiments  is  that  they  are  performed  at  a  limited  number  of 
locations  and  give  average  crack  velocity  values.  This  limitation  is 
more  important  when  using  small  specimens  because  in  these  specimens  the 
velocity  gradients  are  larger.  Further,  the  strain  gages  cannot  detect 
rapid  transients  in  the  crack  velocity  caused  by  stress  waves  reflected 
from  the  boundaries  and  impinging  on  the  crack  tip.  Again  we  expect  that 
these  transients  would  be  more  pronounced  in  small  specimens  than  in 
larger  specimens.  However,  because  the  average  crack  velocities  measured 
during  the  experiments  never  exceeded  600  m  s”^,  we  think  that  the  stress 
wave-crack  tip  interactions  may  have  been  small  in  the  CT  experiments. 

Thus  we  conclude  that  the  strain  gage  method  gives  reliable  first- 
order  estimates  of  the  crack  velocity,  but  is  not  suited  for  studying 
details  of  the  crack  tip  dynamic  unless  an  impractical  high  number  of 
gages  are  used.  With  respect  to  the  CT  crack  extension  experiments, 
tests  CT3  and  CT6  yielded  consistent  crack  velocity  data  in  the  range 
where  most  of  the  temperature  measurements  were  performed,  although 
insufficient  data  were  available  to  characterize  precisely  the 
deceleration  phase  (Figure  H-7). 

In  contrast  to  the  strain  gage  technique,  the  heat  of  fracture 
approach  to  measuring  propagation  toughness  is  independent  of 
constitutive  assumptions  for  the  material  and  can  be  used  even  in  the 
presence  of  plasticity.  In  the  case  of  significant  plastic  deformation 
at  the  crack  tip,  however,  the  finite  width  of  the  heat  source  must  be 
taken  into  account  when  reducing  the  temperature  data.  We  did  this  in 
the  present  study,  but  for  the  sake  of  simplicity,  we  assumed  a  uniform 
intensity  of  the  heat  source  over  the  plastic  zone  width.  This 
uniformity  is  not  true  in  reality  because  plastic  strains  decay  rapidly 
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with  the  distance  from  the  crack  tip;  hence  plastic  dissipation  Is  most 
intense  at  the  crack  tip  and  decreases  as  the  boundary  of  the  plastic 
zone  is  approached*  The  assumption  of  uniform  heat  source  intensity  will 
cause  an  underestimate  of  the  propagation  toughness* 

The  heat  of  fracture  approach  is  based  on  a  one-dimensional  heat 
flow  interpretation  of  the  experimental  data.  This  interpretation 
implies  that  energy  dissipation  at  the  crack  tip  is  the  same  over  large 
distances  of  crack  extension*  Assuming  that  the  propagation  toughness  is 
a  function  of  the  crack  velocity,  the  crack  velocity  must  be  constant 
through  the  region  where  the  temperature  measurements  are  performed. 

This  requirement  was  met  for  most  of  the  data  points,  except  those  taken 
near  the  arrested  crack  tip. 

The  heat  of  fracture  approach  is  a  surface  measurement  and  it 
measures  the  dissipation  in  the  plane  stress  plastic  zone.  As  a 
consequence,  it  may  overestimate  the  plane-strain  propagation  toughness, 
particularly  when  significant  shear  lips  are  formed.  This  point  can  be 
assessed  by  comparing  the  KjD  values  obtained  from  transducers  placed  at 
various  distances  from  the  crack  plane.  If  the  heat  dissipation  at  the 
surface,  is  significantly  different  from  that  in  the  interior,  there  will 
be  a  three-dimensional  heat  flow  field  after  crack  extension. 

Temperature  transducers  placed  farther  from  the  crack  plane  will  be 
affected  differently  than  transducers  placed  closer  to  the  crack  tip.  In 
our  experiments,  the  transducers  that  were  placed  3  mm  and  7  mm  from  the 
crack  plane  gave  results  that  are  within  10%  suggesting  that  there  are  no 
strong  three  dimensional  effects.  In  addition  the  potential  overestimate 
caused  by  the  surface  measurements  is  at  least  partially  offset  by  the 
underestimate  introduced  by  the  finite  width  heat  source  correction  and 
by  the  definition  of  the  maximum  temperature  elevation  caused  by  crack 
extension. 

A  final  comment  about  the  heat  of  fracture  approach  concerns  the 
type  of  transducers  most  appropriate  for  this  type  of  experiments. 
Thermistors  have  a  better  resolution  than  thermocouples,  but  they  are 
more  difficult  to  use  because  they  have  a  nonlinear  response  for  large 
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temperature  fluctuations  and  because  they  are  more  difficult  to  bond  to 
the  specimen.  This  interfacing  difficulty  appears  to  affect  their  time 
response.  Comparison  of  the  performance  of  the  thermistors  and  the 
thermocouples  in  our  experiments  indicates  a  systematic  difference  (20% 
or  less  in  the  temperature  measurements)  between  the  results  of  the  two 
transducer  types  that  could  not  be  explained.  Beside  the  difference  in 
the  absolute  amplitudes,  the  temperature  histories  agreed  very  well,  with 
the  thermocouple  displaying  a  faster  time  response.  For  this  reason  and 
because  of  their  ease  of  application,  thermocouples  should  be  used  in 
future  heat  of  fracture  experiments. 

In  summary,  the  CT  crack  extension  experiments  failed  to  provide 
K^yn  measurements  because  a  large  plastic  zone  at  the  crack  tip 
invalidated  the  assumption  underlying  the  strain  gage  approach.  However, 
reliable  estimates  of  the  propagation  toughness  were  obtained  with 
the  heat  of  fracture  approach.  First-order  estimates  of  the  crack 
velocity  during  the  steady-state  crack  extension  phase  were  obtained  from 
the  strain  gage  measurements  with  an  error  bound  of  50  m  s-^. 

The  propagation  toughness  measured  in  the  12.7-mm-thiek  CT  specimens 

varied  between  172  MPa  at  a  crack  speed  of  575  m  s-1  and  110  MPa 

m1/2  at  a  crack  speed  of  330  m  s-1  with  an  estimated  scatter  of  +_  30  MPa 
12 

m  .  These  values  are  more  than  twice  the  value  of  the  dynamic 
initiation  toughness  and  they  are  also  significantly  higher  than  the 
values  reported  by  Kobayashi  and  Dally^  and  Rosakis  et  al.^2  However, 
other  investigations  of  4340  steel,  in  a  comparable  hardness  condition 
and  for  comparable  crack  velocities,  have  also  yielded  high  propagation 
toughness  values. 1^,19  Furthermore,  the  propagation  toughness  values 
obtained  with  the  heat  of  fracture  measurements  are  consistent  with  the 
values  reported  in  Appendix  G.  As  discussed  in  Appendix  G,  and  following 
the  work  of  Hahn  et  al.,  we  think  that  the  large  increase  in  toughness 
above  the  initiation  toughness  is  caused  by  the  formation  of  shear  lips 
and  by  the  associated  curvature  of  the  extending  crack  front. 

Apparently,  the  resistance  to  crack  extension  is  much  more  sensitive  to 
the  degree  of  constraint  than  the  resistance  to  crack  initiation.  Also, 


the  toughness-crack  velocity  relationship  may  depend  on  the  initiation 
conditions  (such  as  the  intial  notch  root  radius)  because  the  size  of  the 
plastic  zone  will  affect  the  size  of  the  shear  lips  along  the  extending 
crack  path,  over  distances  that  are  significant  compared  with  the  initial 
plastic  zone  size. 
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